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over  finite  ground. 


SUMMARY 


Previous  thin-wire  antenna  programs  have  either  been  very 
specialized  or  all-encompassing.  A beginning  or  occasional 
user  does  not  need  expertise  in  programming  to  gain  insight 
into  wire  antenna  structures,  using  this  general  purpose 
user-oriented  code.  The  revisions  contained  herein  correct 
deficiencies  of  handling  the  image  problem  in  the  original 
code  and  improve  the  accuracy  of  calculations  of  structures 
over  finite  ground. 
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EXPLANATION  OF  REVISIONS  TO  ASAP 


Since  the  issuance  of  the  original  ASAP  report  in  December 
1974,  and  distribution  of  the  source  program,  feedback  from  many 
users  and  direct  assistance  from  Dr.  Robert  Bevensee  of  Lawrence 
Livermore  Labs,  University  of  California  enabled  the  preparation 
of  this  revised  version.  During  the  past  2 years,  no  additional 
discrepancies  have  arisen  and  it  is  now  felt  that  the  code  can  be 
safely  called  "revised". 

The  nature  of  the  improvements  was  as  follows: 

1.  To  correct  the  manner  in  which  structures  were 
treated  when  elevated  over  a ground  plane. 

2.  Improvements  and  corrections  in  calculating  finite 
ground  plane  effects. 

3.  Text  corrections  in  explanation  of  ground  effects. 

4.  Instructional  changes  in  the  DESCRIPTION,  FREQUENCY, 
and  CHANGE  cards. 

5.  Corrected  sample  problem  outputs. 

The  user  should  note  that  some  facilities  will  experience 
printing  errors  when  NEAR  FIELDS  are  called  for.  Statement  58  in 
the  MAIN  program  should  be  suitably  rewritten  if  that  occurs. 
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!•  INTRODUCTION 


Although  many  thin-mire  computer  programs  have  been 
developed  for  the  purpose  of  analyzing  antennas  and 
scatterers,  few  of  these  programs  have  been  directed  toward 
the  student  of  electro-magnetic  theory.  The  majority  of  the 
programs  are  directed  to  the  engineer  or  advanced  student 
for  the  purpose  of  analyzing  designed  structures  or 
verifying  experimental  data. 

The  purpose  of  the  study  is  to  develop  a computer 
program  by  modifying  an  existing  computer  code  which  can  be 
utilized  as  an  educational  method  to  develop  insight  into 
radiating  structures  by  the  beginning  student  of 
electro-magnetic  theory. 

The  modified  Ohio  State  University  Antennas-Scatterers 
Analysis  Program  (OSUHOD  or  ASAP)  is  directed  toward  the 
beginning  student  who  does  not  yet  have  the  expertise 
necessary  to  manipulate  the  input  data  for  proper  execution 
of  the  larger  more  comprehensive  analysis  program.  Even 
though  ASAP  is  small  in  core  requirements  and  is  fast  in  tun 
time,  it  is  capable  of  analyzing  structures  to  assist  the 
engineer  with  design  problems. 

Since  the  resulting  program,  ASAP,  is  primarily  directed 
toward  students,  the  program  has  been  limited  to  structures 
which  contain  less  than  50  monopoles  (segments) , no  longer 
than  one-fourth  of  a wavelength,  and  which  have  less  than  51 
nodes  (intersections  and  endpoints) . If  a ground  plane, 
either  perfect  or  finite  is  present;  the  stated  limits  above 
are  halved  due  to  the  generation  of  an  image  structure. 
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II*  ORIGINAL  PROGRAM 


A.  THEORY 

Reference  1 presents  the  electro-naqnetic  theory  for  the 
analysis  of  antennas  and  scatterers  in  an  isotropic,  linear, 
and  hoBogeneous  ambient  aediua.  The  analysis  is  performed  in 
the  frequency  domain  with  an  excitation  caused  by  either  a 
generator  or  an  incident  wave. 

In  the  analysis,  a piecewise-sinusoidal  expansion  is 
used  for  the  current  distribution.  The  matrix  equation 
Z 1 = V is  generated  by  enforcing  reaction  tests  with  a set 
of  sinusoidal  dipoles  located  in  the  interior  region  of  the 
wire.  Since  the  current  distribution  has  the  same  form  as 
the  expansion  mode,  this  formulation  is  known  as  the 
"sinusoidal  reaction  technique". 

B.  COMPUTER  PROGRAM 

Reference  2 presents  the  computer  program  corresponding 
to  the  theory  presented  in  Ref.  1. 

1.  Input  Poraat 

In  the  program,  the  input  data  must  specify  the 
frequency,  wire  radius,  wire  conductivity,  the  parameters  of 
the  exterior  medium,  coordinates  of  the  points  to  describe 
the  shape  and  size  of  the  wire  configuration,  a list  of  the 
wire  segments,  and  the  indicators  for  the  various  outputs. 
Table  1 is  the  input  data  necessary  to  analyze  a half-wave 
dipole. 


2.  Output  Format 

In  the  original  fora,  the  only  outputs  which  could 
be  requested  by  the  input  data  stream  are  the  following: 
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a.  Antenna  Problems 


(1) 

Current  Distribution  on  the  Structure. 

(2) 

Input  Impedance. 

(3) 

Radiation  Efficiency. 

(“) 

Near-Zone  Field. 

(5) 

Far-Zone  Field. 

b. 

Backscatter ing  Problems 

(1) 

Absorption  Cross  Section. 

(2) 

Scattering  Cross  Section. 

• 

(3) 

Extinction  Cross  Section. 

Scattering 

(4)  Complex  Elements  of  the  Polarization 

Matrix 

c. 

Bistatic  Scattering  Problems 

Echo 

Area. 

Table  2 is  an  example  of  the  output  data  available  for 
data  of  table  1. 

3.  LIMITATION 

Although  the  program  can  analyze  a structure  with  up  to 
50  segments,  55  points  and  60  dipoles  modes;  it  can  not 
analyze  a structure  in  the  presence  of  a finite  ground 
plane. 
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AN  EXAMPLE  OF  THE  OUTPUT  DATA  FOR  THE  ORIGINAL  PROGRAM 


Ill*  MODIFIED  COHPUTEE  PROGRAM 


A.  Input  Format 

As  illustrated  in  table  1 the  format  for  the  input  data 
cards  is  not  self  explanatory.  This  format  can  be  determined 
by  referring  to  the  FORMAT  statements  of  the  program  of  Ref. 
2.  Since  the  modified  program  is  directed  toward  the 
student,  the  input  data  format  was  changed  to  allow  free 
format.  Reference  2 was  written  in  a form  which  permitted 
modifications  to  allow  flexibility  in  specifying  input  data 
for  the  analysis  program.  Appendix  B,  titled  "User's 
Manual",  discusses  the  input  data  cards  necessary  for  proper 
execution  of  an  analysis  problem.  Appendix  B is  self- 
contained  and  may  be  used  independently  of  the  remainder  of 
this  document. 

B.  Outpat  Format 

In  the  original  computer  program,  the  absence  of  labels 
encumbered  the  output  data  and  lessened  the  usefulness  of 
the  program.  To  improve  the  usefulness  of  the  modified 
version,  detailed  labels  were  added  to  the  output  data.  As 
with  the  input  data.  Ref.  2 was  written  in  a form  which 
enabled  modification  to  allow  more  specific  output  data  for 
the  analyzed  problem,  with  the  addition  of  the  polar 
plotting  package,  the  far-zone  electric  field  intensity 
polar  radiation  and  reradiation  patterns  can  be  plotted.  A 
sample  problem  can  be  found  on  page  120  in  Appendix  B, 

User's  Manual. 

C.  Finite  Ground 

To  enable  the  student  or  the  engineer  to  have  an 
improved  analysis  program,  the  finite  ground  effects  were 
added  to  A3AP.  The  theory  corresponding  to  the  ground 
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effects,  which  utilize  Fresnel  reflection  coefficients,  is 
discussed  in  Appendix  A,  titled  "System  Manual".  Also 
discussed  in  Appendix  A is  the  modified  computer  program  and 
the  corresponding  theory.  The  electro-magnetic  theory  was 
developed  in  Refs.  1,  2,  and  3;  and  it  is  restated  with  its 
corresponding  computer  code  to  assist  in  the  understanding 
of  the  methods  applied.  Appendix  A is  self-contained  and 
may  be  used  independently  of  the  remainder  of  this  document. 


IV*  CONCLUSION 

The  addition  of  ground  effect  techniques  to  the  original 
program  did  not  alter  the  accuracy  or  the  computational 
capabilities  of  the  program.  The  ground  effect  techniques 
utilized  the  results  of  the  original  program  and  modified 
these  results  to  account  for  the  effects  of  the  presence  of 
the  finite  ground. 

To  verify  the  numerical  results  of  ASAP,  the  input 
impedances  of  both  a horizontal  and  a vertical  dipole  were 
compared  to  the  solutions  of  the  exact  form  of  the 
Sommer field ' s equation.  As  can  be  seen  in  table  3 the 
finite  ground  treatment  of  ASAP  agrees  favorably  with 
Sommerf ield ' s solutions.  The  ASAP  finite  ground  results  are 
also  in  excellent  agreement  with  the  previous  computer 
solutions  of  Refs.  4 and  5. 

V*  RECOMMENDATIONS 

Although  the  program  is  a general  analysis  tool  for 
students,  several  future  modifications  will  enhance  the 
program  as  a design  tool  for  engineers.  These  items  include: 
varying  the  wire  radius  on  the  structure;  incorporation  of 
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non-radiating  elements  such  as  transmission  lines;  varying 
the  wire  insulation  radius,  conductivity,  and  dielectric 
constant;  and  a geometry  generation  package  such  as  dipole 
array  or  helix.  One  aaior  cnange  that  would  both  improve  tne 
speed  and  reduce  the  core  requirement  is  that  of  symmetry. 

No  attempt  was  made  to  utilize  the  symmetry  in  the 
admittance  matrix  when  the  ground  plane  is  present.  If 
symmetry  were  applied,  the  structure  size  limit  with  the 
ground  plane  present  would  be  approximately  that  of  the 
structure  without  the  ground  plane. 
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STSTEH  MANUAL 


INTRODUCTION:  The  Antennas-scatterers  Analysis  Program 

(ASAP)  for  thin  Mire  structures  in  a homogenous  conducting 
medium  performs  a freguency  domain  analysis  of  antennas  and 
scatterers.  The  program  is  applicable  in  the  presence  of  a 
ground  either  perfect  or  finite.  This  appendix  Mill  describe 
the  computer  program  which  accomplishes  this.  Although  the 
program  was  written  for  the  IBM  360  computer  system  it  can 
be  executed  on  another  system  with  minor  modifications. 

A piecewise-sinusoidal  expansion  is  used  for  the  current 
distribution.  The  matrix  equation  ZI  = V is  generated  by 
enforcing  reaction  tests  with  a set  of  sinusoidal  dipoles 
located  in  the  interior  region  of  the  wire.  Since  the  test 
dipoles  have  the  same  current  distribution  as  the  expansion 
modes,  this  may  be  regarded  as  an  application  of  Salerhin's 
method.  Rumsey’s  reaction  concept  was  most  helpful  in  this 
development,  and  therefore  the  formualtion  is  known  as  the 
"sinusoidal  reaction  technique". 

The  main  routine  and  each  subroutine  is  discussed 
separately  in  this  appendix.  The  writeups  for  the 
subroutines  are  arranged  alphabetically  by  subroutine  name 
after  the  main  program.  Each  of  the  discussions  includes  the 
purpose  of  the  subroutine,  brief  description,  and  a listing. 
After  the  subroutine  writeups  is  a table  of  the  more  common 
symbols  used  in  this  program. 

The  input  data  and  program  limits  are  discussed  in 
detail  in  the  next  appendix  titled  "USERS  MANUAL". 
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GROUND  EFFECTS:  In  the  aodified  antenna  analysis 

# computer  program  finite  and  infinite  ground  effects  were 

added  by  using  the  reflection  coefficient  technique.  The 
method  in  which  tnis  technique  was  used  required  the 
generation  of  an  image  structure.  In  this  section  the 
reflection  technique  will  be  discussed  in  detail. 

In  order  to  apply  ground  effects  to  the  electric  field, 
the  field  for  the  image  structure  was  first  calculated  as  if 
a ground  were  not  present.  Then,  the  field  was  decomposed 
into  parallel  and  perpendicular  components.  (A  parallel 
component  is  the  component  which  is  parallel  to  the  plane  of 
incidence.  A perpendicular  component  is  one  which  is 
perpendicular  to  this  plane.  The  plane  of  incidence  is  the 
plane  containing  the  normal  to  the  reflecting  surface  and 
the  incident  ray.) 


Z 1 


Considtic  an  image  monopole  with  the  electric  field  in 
the  1 direction.  The  ray,  6,  is  a vector  which  is 
perpendicular  to  1 and  passes  thru  the  point  of  interest.  To 
apply  reflection  technique,  the  plane  of  incident  must  be 
found.  It  is  advantageous  to  define  a new  coordinate  system 
(a«/3/Z)  where  a and  p are  parallel  to  the  xy  plane  with  a in 
the  plane  of  incident  and  perpendicular. 


If  the  direction  cosines  (cosx,  cosy,  and  cos^)  are 
known,  it  can  be  shown  that  the  components  of  the  field  in 
the  a/3  (xy)  plane  have  the  following  relationship: 


, 



— — 1 

1 1 

E 

cos  p sin  t 

E 

X 

sin  ♦ -cos  (> 

E 

J 

_ _ 

- y .. 

where  ^ = arctan (cosy/cosx)  . 


Now  the  reflection  coefficients  for  the  interface  can  be 
applied  as: 


1 1 (R) 
E 


R 

N 


I 


±(R) 

E 


where  R 

I I 

Applying 


and  R will  be  defined  later  in  this  section, 
i 

the  matrix  equation  above  yeilds: 
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X 

cos  9 

sin  9 

1 1 (R) 

E 

e‘"' 
_ y _ 

sin  9 

-cos  9 

1 (R) 

E 

(the  square  natrix  is  unique,  in  that,  the  inverse  is  equal 
to  the  original  matrix) . Since  the  image  direction  is 
opposite  to  the  original  nonopole,  that  is. 


(1  X z)  . , = - (T  X z)  . , 

original  mage 

the  2 component  of  the  field,  which  is  in  the  plane  of 
incident,  is  given  by: 


E 


(R) 


z 


E . 

z 


Erom  electro- magnetic  theory  the  reflection  coefficients 
for  the  fields  in  medium  (1)  at  the  interface  with  another 
medium  (2)  are  defined  as: 

for  perpendicular 

cos  9 - 6 ' - sin*  9 

R - — ■ , ij, 

H COS  9 y & ' - sin*  9 


and  for  parallel 

6*  cos  9 - y 6 ' - sin*  9 
V &'  cos  9 ♦ y & ' - sin*  9 

where  9 is  the  angle  of  incident  as  measured  from  the  normal 
to  the  interface  and 

6*  * (6  ♦ o/j{o)/(6  ♦ a /j<o) 

2 2 I I 

where  the  subscripts  correspond  to  the  mediums  above. 
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To  determine  the  relationship  between  R , R and  R , R 

II  j_  V f 

a perfect  ground  (&  = 0,  a*x)  nas  investigated. 

r 

limit  H = -1 


limit  H = ♦! 

V 

But,  for  a perfect  ground  the  contributions  to  the  field 
from  the  image  monopole  wojld  be  equal  to  the  field  of  tne 
original  monopole  but  opposite  in  sign  due  to  the  chosen 
reference  direction. 


1 I (R) 


II  11 

: = - E 


±(B) 


therefore 


H = - a 

I I V 

d = H 

1 H 


In  summary,  the  contribution  to  the  electric  field  of  a 
monopole  over  a ground  plane  at  a given  point  is  given  by: 

(R)  (H)  (R)  (R) 

E = E cosx  ♦ E cosy  ♦ E cosz 
z y z 

where 


= R E cosz  ♦ (P  - R ) E cosx  cos*  t 

_L  I I X 


♦ (B  - R I E cosy  sin  # cos  # 

I I X 
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E cosy 


E 


(R) 

y 


R 

I I 


E cosy  cos^  t 


(R  - R ) E cosx  sin  ^ cos  ^ 

I I X 


E 


(R) 


z 


R E cosz 


where  E is  the  field  without  the  ground  plane  present  and 

- &•  cos  0 ♦ \/^ S'  - sin2  6 

R = 7 ' 

I I S'  cos  6 ♦ v &*  - sin*  0 

cos  0 - y' ~ sin*  0 

H = — ^ 

cos  0 y S ' - sin*  0 

S'  = S ~ j (a/s  6>) 
r 0 
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MAIN 


PURPOSE:  to  control  the  input,  output,  and  the  flow  of 

calculations. 

HETHOD:  The  main  prograa  controls  the  flow  of  the 

required  calculations  by  calling  only  a few  subroutines. 
These  subroutines  in  turn  call  other  subroutines  which 
actually  do  the  required  calculations.  The  order  of  the 
calling  sequence  is  diagramed  after  the  listing  for  the  main 
program. 

The  DIHENSION  statements  at  the  beginning  of  the  main 
routine  provides  the  required  storage  for  a wire  structure 
with  up  to  50  segments,  60  nodes  and  60  dipoles  without  the 
presence  of  a ground  plane.  If  a ground  plane  is  present 
one-half  of  the  reserved  storage  is  required  for  the  image, 
therefore  a wire  structure  with  up  to  25  segments  and  30 
nodes  can  be  analyzed. 

NH  denotes  the  actual  number  of  monopoles  (segments) , 

INH  is  the  corresponding  dimension,  and  the  dimension  for 
CG,  VG,  and  ZLD  is  twice  INM.  The  second  subscript  for  MD 
always  has  a dimension  of  4 to  correspond  to  the  number  of 
segments  meeting  at  a given  node. 

N denotes  the  number  of  simultaneous  linear  equations 
and  ICJ  is  the  corresponding  dimension.  The  dimension  for  C 
is  (ICJ  * ICJ  ♦ ICJ)/2. 

In  the  statements  above  statement  4,  the  initial 
conditions  and  defaults  are  established.  After  calling 
subroutine  READ  to  determine  the  input  parameters,  the  IF 
statements  output  the  parameters  to  be  used  for  the 
calculations.  In  the  DO  LOOP  ending  at  statement  7,  the 
the  input  data  of  the  structure  geometry  is  stored  in  order 
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to  recall  if  the  structure  is  to  be  moved  for  ground  plane 
calculations. 

After  the  image  structure  is  generated  and  structure 
location  is  moved,  subrou* ine  SORT  is  called  to  determine 
the  aipole  modes.  Prior  'O  calling  SCANT,  the  load  and 
generator  information  is  established. 

Subroutine  SCANT  is  then  called  to  calculate  the 
elements  of  the  impedance  matrix.  If  FEEDS  or  CENERATORS 
are  specified  by  the  input  data  stream,  subroutine  GANTI  is 
called  to  solve  for  the  current  distribution  due  to  these 
forcing  functions. 

In  the  DO  LOOP  ending  with  statement  29,  subroutine 
GNFLD  is  called  to  calculate  the  near-zone  field  for  the 
current  distribution  of  the  subroutine  GANTI. 

The  subroutine  GFFLD  is  called  for  the  far-zone  field  of 
the  current  distribution  of  the  subroutine  GANTI  in  the  DO 
LOOP  ending  at  statement  35.  The  subroutine  GFFLD  is  called 
again  in  DO  LOOPs  ending  at  statements  42  and  51,  if 
bistatic  and  backscattering  calculations  are  requested  by 
the  input  data  stream. 

CALLS  TO:  GANTI 

GFFLD 
GNFLD 
POLPRT 
READ 
SCANT 
SORT 
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oc^(N$!aN  KI60I,  r<60»t  2<60}t  xGi60>f  yatbOif  /&i6p) 
plM€N$fON  n(60f.  12(60)1  I|(6pi,  JA(63)«  jiUoit  KFLAGOO) 
plMFNSION  CPHH^OOI,  CTMfT(500T 

OlMtNSlON  0ATfl(}60),  0ATr2{360).  06TY}()60it  0*TY4()60i 
OIMFNMQN  0(50).  |A(50)|  16(50),  1SC(50)*  MOl^OtA),  HoTSO) 
OlMtNSlON  L20(ttdl«  KGCnUO) 

CO>’M-N  IWt 

0IH(nS1UN  XNPI50).  YNP(50»t  2NP(60) 

CLMPIEX  C(  1630) 

COMPIEX  COATtI i00),CPAT2l 500)fCOATJ(50p) >CpAT6(500l 
COMPi  EX  C J(60)  lEPI  60T.EPP(60)f|T(60i  •€>  < (<>0) 

COMPLEX  CC04  50I  ,SGO(65)iCC(  IOO  | , VCI  100)  t Ko(  1001 


^oSplI!  EP^s!tpU[t^?ijETTS.EX,eY.E2 

" 0/6.eS6C-12,1.2544E*6/ 


DATA  EO.U' 
NGEN  • -I 
JOAO  - -J 
LOAD  • -I 
6M  > *1 
ICARO  • 0 
AM  " • 1 
(FLAG  • 0 


VOlTI  U 
GT 


HG 

NN  • 0 
NP  • 0 
. 0 


0. 


SIG2  • -1. 

T0|  • “1. 

SIG3  • -I 
ER3  - I 
T03  • 0. 

CMM  • 50. 

15^  ■ loo. 

INM  • 50 
iCJ  ■ 60 
WaItC  (6«74) 

DO  2 |•l•30 
2 RFLAGI I I - -1 

00  3 J-it  1 w 


000) 

0002 

0003 

0004 

0005 

0006 
0007 

oooa 

0009 

0010 
0011 
0012 

8§U 

0015 

0016 

mi 

0010 

0020 
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0022 

0023 

0024 

0025 
0020 
0027 
002  8 
0029 
0 0 } 0 

0031 

0032 
003  ) 
C0)4 

0035 

0036 
00)7 
0036 
00)9 
0040 

88tl 

004) 

0044 

0045 

0046 

0047 

0046 


ZLO( J»  • ( .Ot •0) 
J3  • J»tNM 

VG( JJl  • I .0. .0) 

3 2101 JJJ  • I .A# .61 


Nf  FP 
N61P 
NBAP 
AFFP 
AFf  T 
A8  IP 
AbIT 
AbAP 
A6AT 
STEP 


• 0 
• 0 
• 0 

« loop. 
• 1000. 
• 1000. 
• 1000. 
■ 1000 
- I0-- 
1. 


000. 


Y,YNP.7,2Ll6,2NP,STEP| 


S - 

ANM  • 0 

CALL  AEAO 
IftfLAOiAGe 
jFfT.AftlP. 

3SF,PhM 
«X.XNP,Y,T'^k, 

WRITE  (6,5^i 
)F  (MSG. I T.n  S3  TQ  5 

(MSG.f4.ll  rfAlTE  16.70)  KFlAGOO) 
IIFLXG.iG.4l  GU  TO  I 
(IFIAG.F4.5)  STOP 
IAM.lI.OI  hAITE  (6 
(AM.lT.O)  go  to  6 
( (rwi.Gr.OI.ANO.(NP.GT.O) ) CO  TO  7 
WAITE  1 6|116) 

IF  (1ELAC.E~ 

MSG  • 2 
GO  TO  4 
WAITE  (6*114) 


.1271 


0.1)  CO  TO  1 


0 04  9 
0050 

m 

OQii 

0p54 

0055 

0056 

0057 
0056 

0059 

0060 

006  ) 
00o4 

0065 

0066 
00  6 7 
0066 

0069 

0070 

007) 

0074 

0075 

0076 

0077 
0076 
0079 
0060 
0061 
0082 

008) 
0084 

0065 

0066 
0067 
0066 

ofli 

0092 

009) 


IF  IKFLAOUll.tJ.U  MRITF  (6.99)  T03 
IF  |KFLAGI26).NE.n  R I T C (6.12^) 

If  il  SSlfl  It’fHI 

IF  ( I (GRD.Gf . I ) . AMO.l KFLAGt 25) .EC.  ))  rfAlTE  (6.1261 
IF  <(  iGRl).OT.O).ANO.UFLAG(25).|Q.  ))  rfRlTE  (6.126) 
If  (AFLAG(2n.EO.()  MRITE  16, 121)  NT 
wpin  (6,111) 


ER6.S1C6 

HCT 


i I.AW9II  191  lOI 

wRlTF  (6.1(11 

IF  (KFLa6{2a).GT.O)  kRITE  (6,119)  ( I tOl  t ) • 2L(.0(  I ) • 1 • 1 » lO  AO) 

IF  (KFCAG(  16).GT  .0)  MRITE  (6,118)  (L20(  I ) , 2LL0(  1 ) . !•  1 .(.OAO ) 

WRITE  (6.111) 

IF  (AFl aGI 23  I .GT  .0)  WRITE  (6,120)  (KGEN(  1 ) , VOLT  ( 1 I , I -1 . NGCN) 

IF  (kFl A^( 1 3) .GT  .0)  WRITE  (6,97)  ( KGEN(  I I . VOLT  1 1 ) . !• 1. NGEN) 


WRITE  (6.111) 
WRITE  (6.116) 
WRITE  (6,98) 
WRITE  (6*112) 
F (RfLAG(22l 
(KFLAOl  15) 
(KFLAGl  16) 


n: 

E4. 


)AFlAGli;).EU. 
(V^HAGI  16)  .E^. 


WRITE  (6, 110) 

WRITE  (6,991 
WRITE  (6,100)  PHA|.PHAF,ThAI,THA 
WRITE  16,101)  PMIKPHIF.TM]]  ‘ ‘ 


i WRIT!  (6,102)  PHil.PHSF.THil 
) WRITE  (6,103)  (XnM I ),VNP(I) 
WRITE  16,i65|  AFFP 


Sfifl  I 

ssrri  IS: 

WR  " ‘ 

WRi.c  lO.J--- 

ANO.(lSCAt.LT.O))  WRITE  (6,73) 
' TO  129 


n IS: 


06)  AFFT 
09)  ABAP 
08)  ABAT 


iS^t 


Sll? 


(RFL AG( 19) .E 
(AFFP. LT. 500. 

(AFFT.L  r.500. 

(ABAP. lT. 500. 

(ABAT.LT.500. 

<AeiP.LT.500. 

( ARl T.L  T.500. 

( ( IB  ISC.GT.O)  . . . 

(kFLAGTaI.lT .1 ) GO  to  129 
OC  UB  1*1  • i NM 

128  1 SC(  I 

129  f H2  • FHC*1  .E6 
ONfGA  « 

IF  (S1G2.LT.0,)  EP2«ER2*E0wCMPLX( 1..-TD2) 

(T02.LT.0.J  EP2  • CMPLX ( 1r2 WEO.-i 1G2/6mEGA ) 
(SIG3.LT,0. ) EP3«ER36E06CMPLX(  U.-TOal 

li!  •;  BllF*-*'^**®!***! 

WRITE  (6,121)  INT 


I HEAR) 


111 : 

IF  (KFLAGTuI.NE.I)  go  to 


0097 

0098 

mi 

stsi 


m 


03 

06 

05 

06 

07 

08 
09 

0 

1 

2 
3 
6 

5 

6 

7 

8 
9 
0 

i\ 

n 

il 

il 

ii 

33 

36 

5 

6 

7 

8 

9 

60 

61 

S! 

66 


NPG  • NP 
NNC  • NM 


0165 


00  8 
XG( 
VG(  . 
8 ZG(  I 


l-l,NPC 

II : sill 


9 00  10  I "liNPa 
X(1)  • xGll) 

Y(l)  • VG(  1) 

10  2(1)  • 2G( 1) 

NP  • NPG 
NN  • NMG 

|WL  • 0 

If  (IGRO.LE.O)  GO  TO  15 
UT  UP  IWAGE  FOR  GROUND  PLANE 
2HIN  • 2( 1 ) 

R . 0 

IF  (21 r ).LT.2HlN)  2HIN«2(II 

00  (1  r-i|NP 

2(1)  • /Cl )*HCT 

IF  (2( 1 ).GT.1.E-60I  GO  TO  II 

] Wl  • 1 WL  * 1 

11  CONT INUE 

IF  (2N1N.GE.O.O)GOT012 
WRITE  ( 6.1  I 71 
IF  (IFLaG.EO.I)  go  TO  I 

IF  (Iflag.eo.2)  stop 

WSG  • 2 
GO  TO  6 

12  00  13  J-l,NH 
J6NI‘ 


13 


• J6NM 
lAIRl  • IA(J) 

IWLP  • IML*1 

00  16  I•IIILP•HP 
J • (♦NP-IWL 

?lil  : Jill 


0169 


26 


14  ZIJI 


II 


15 


le 

19 

20 
21 


22 
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KKM  ■ NN*  I 
HH  • 2*NH 

HP  • 2*np-;kl 

CALL  $0AI  UA.  15»i  1 1 I2«  13.  JA,  JB,MO.W«NMfNPtN«MAX,H|N«  ICJtlNMI 
If  (iAX.LE.AI  GO  TO  16 
MAITE  (6.71) 

IF  (IFLAC.rg.l)  GO  TO  I 

if  (iflag.fo.M  stop 

HSG  • 2 
GO  TO  4 

IF  (MlN.GC.l)  GO  TO  17 
WPlH  (6,>M 
IF  lIFLAG.eU.l)  &Q  I 
If  (IFLAG.fcU.H  WOP 
«SG  • 2 
GO  TO  4 
WRITE  (6.56) 

F (MA>[.6T.<».0R.MIN.LT.l.0R.W.GT.|CJ)  60  TO  54 

2 - ( 

(LOAO.GT.O)  go  TO  19 

00  18  l«l ,NW 
2L0( U • (0..3.) 


mn  I 

I’m 


IF  (NCFN.GT 


id  TO  21 


00  20  l-l.NX 

veil)  • ( 6..0. I 

RH  • NM 

IF  (IGRD.GT.OI  KN  - NM/2 
J • 1 

ANTENNA  CAtCUlATlONS 
IF  (LOAD. LE. 01  GO  TO  24 
IF  (RFLAGI 24)«GT .0)  GO  TO  22 

00  23  J«1.KN 

IF  I ( Ul  J>  .EO.KI  .*N0.  (KFLACI  UI.Q 
IF  liRFLA&(24I.Gr.0I.AN0JlG 

Continue 
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UOIKARN 


UUSlill 
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019? 
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0200 
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0204 
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0207 
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P2‘ 
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24  II  (NGEN.LT.OI 
RN  • NH 
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GO  TO  27 
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25 


00  26 
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25  00  26  I-l.NGCN 
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IF  ( (Ki LAGU  3).GT .OI.ANO.I i&RO.Gt.O)) 
If  l(  URO.GT  .01  .ANO.(KFLAG(  231  .GT.OU 

26  CCNT INUE 

27  CALL 
I fCGO 

If  IN  . 

IF  ((Flag. fg. 21 
NSG  • 2 

IF  (IFlAG.EG.l) 

GO  b 4 

28  If  (NGfN.lE.O)  GO 
WRITE  16,751 
WR  I TE  16,  76) 

WRITE  (6,77) 

NR17E  (6.82) 

CALL  .ANtj  ( lAijR 


0)J  VGI JI-VOiT( II 
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NP.AN.GN.C 
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IF 

MR 
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DO 

29 

IP 

■ 
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• 
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IK 

• 
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CALL 
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.LE .01 

751 
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771 

iNt  A< 


GO  TO  30 


:i)ll  ir.lll 

29  CONTINUE 
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pATVIl  1 I 
OATY^i  I I 


UlftO 
• 0 
• 0 
m ) 
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Tf  f*l 


Nf'l 
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CC  TO 
39  NPL  • 
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34  Pm  • 

00 

pH 


Tm  • 


APF  P 

• 0. 
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|F  ( I tCKO.CT.01  . AN0.MTh.6T  .901  • ANO.I  Tm.lT  .zrOM  ) 60  TO  35 
2CPFtO<  • •PH.$CD,SCS^fSCST,SPPMtSMM,UPN»STtM,TH,X,V,iilLOf  ZS»€fAfC 

\i 


I •rn.> 

3AN,BkA»  ICkOI 

jpSiE  * ^|$(!FFS. 

IF1CTHA6.6T.UE-I2I  PHSTH*  57. 295TT94AT  AN2  C A IMA6(  C TTSI  tKI  AL  IE  TTS II 


^-'11 


05)5 


FiEPNAC.CT 
F (NPl 
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CCNT INUE 
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If  1 IWA.GT.O)  GO  TO  40 
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|Pi£CSr»$CSP.5CST 
40  CPhHlI  • PM 
CtHflUl  • TH 
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COATXLI  • 6TP$ 
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OATV^IkI  • CAttSlfPTSl 


(4t41)  PH,TH»$PPMf$PTH«STPN«STTM«ACSP»ACST»€CS 


0ATY3IAI 
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CALL  POiPPr  IT.OATYn 
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CALL  POLPRT  lld.OATYAI 
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GO  TO  90 

ISURfSI 


29 


i* 


I 


BtSI  K 


I ’ ■ ■ ■•U  — *• 


ATY^lK)  • tA»51ETTSi 
f (NAl.NE.^I  CO  TO  5 


.^1  CO  TO 
CASSi  epp: 
CABSI  EPT< 
CABSiETP: 
CARSIITT 


90  IE  <NPl. 
OATyini 
0ATV2( 1 I 
OATYIU  ) 
OATYAi  I \ 

91  CONTiNOE 


WPITE  (6.92) 

IF  INPl.lE.Ol  CO  TO  92 
CAIL  POIPPT  IJfOATYll 
call  polppt  <A,0ATY2[ 
CALL  POLPPT  U>0ATY)| 
CALL  POLPKT  t6tOATY4| 
IF  iKfi  ACI lei .NE.I I CO 
92  WRITE  |<»,64l 


SI  CPHI  ( n .CTHCn  |),CO*TlUlfCOAT21ll,CO*TlllltCOA?4ni 


94  IE  IIELAC.EO.U  CO  TO  1 
ir  IIFLAC.eO.2)  STOP 

KKFLAG*0 
A JFI  A0*0 
KPFLAC-0 
KNFLAG-0 

IE  TaPlAGI 13) .cr .01  KKfLAC*! 
IF  (KPL ACI 23).CT .Oi  KJPLAC-I 
IE  UFLACl 14) .GT.Oi  KMFLAC*! 
IF  (KElACI 24). GT.OI  YNFLAG*) 
00  59  l«l.30 
99  KElACI ( ) • >1 


KFLAGiai  - I 
AFLAC120)  « I 
KFLAGI26I  • I 
IF  (KAELAG.oT.OI 
|F  IKjFLAG.UT.Oj 
If  (AMFLAC.Cr.O) 


, . . . J kflaciIaI- 

ir  hi; 

IE  (IFLAC.Io.4)  MPITC  Utll9l 

CO  TO  4 

It  f8S5i}  lloxl'TMt  RfOUTION  IRFICIENCr  IS  • . FIS.T// lOX,  • THE  Tldf- 

lERME  ROWER  IWRUE  IS  • .E 1 5.  E// lOX  t * tHE  RNTEMtA  IHREORNCf  tS 


\ 


ft'  NEAA-FIELO  ELECTRIC  FIELD  INTENSITY  AT  THE  OBSERV 

lATION  POINT  • ,El  I .9*S*  .E  U .9.  • , SEll.Sf  • CX*Y,Z  RESPECTIVELY)  ISt 
2*//i 

59  FORMAT  ( 20l • * EX • • » F 1 9 • 7 • • » J* ■ F 1 9 . 7/ 20X • • EY •• • F I 9 • 7 1 * AJ*fFl5.T/20 
lx. ’EZ •*  .E  19. 7i • ♦J».F|9.r////# 

60  FORMAT  (3X,E5.1,2X,E9.1.3X,El0.4«2X«EL0.4.20X.>(EI0.4»2X)fF6«ltU 

6l^^|>RMAr^O4^^*F0;  INClOENT*/T4lf»PLANE  WAVE 

62*EoSMAr *(• * IN^ i6ENT^iT2  7i*ECH4^ARlA  SIGMA* . T66 , • ABSORPT lON'iTBOt*  EX 
lTlNCTlUN*.rU4.*SCAfrER!NG'/*  PLANE*  , T 2 9 1 * UNC  I DENT  ^SC  AT  TER  EOl  * . I 
24X,3( 9X t 'CROSS  SEC TION* ,6X) /•  WAVE  < 1 9 2X , 3 (1  OX , * F OR • . 1 1 X ) /•  PhI 
3 THE  TA* . iX I • PHI-PHI  * 1 3Xf *PhI>ThETA*.4X«  * THETA'PhI * fZXf • ThETA- THETA 
4* ,31 9X, ‘PhI* , 7X, • THETA* ,4X) ) 

63  FORMAT  llX,2lF9.1. ]X|.l6(£l0.4,2Xi) 

64  FORMAT  (T94.«BACKiCATtER(NG*/*  I NC 1 DENT  * , T 3 7 • • EL ECTR ]C  FIELD  POLAR 
IIZATIJN  SCAtTfRlNG  MATRIX*/*  PL  ANE  * , T49  , • 1 j NCI  OE  NT*  K A?  tEREO)  V3X 
2.  * wave*  ,T  23*  * PHI  -PHI  * .T49,  'PHI  - The  TAST79, » TMETA-PmISTIoJiMMETA- 
3tHrTA*/ « pH  I theta* , 3X,4( 3X. •REAL*  rBX t * ImAG* .BX) ) 

69  FORMAT  llX,2IF9.1.iX),2X.4lEll.9i2X.ni.9i3X)J 

66  format  1T54. •0ISTAfj4*7Tj7.*ELFCtRlC  FIELD  POLARIZATION  SCATTERING 
1 matrix*/*  OBSERVATlQN'.TSOi* IINCIDInT-SCATTEREoI^/*  fOlNT*tlAXi 

1 • FmI-PhI  • .T49,*  phi-theta*  . t76,*  ThEU-PHI  * , TiO| » MhETA-THETA*  /•  P 

6 7^f4rmaT^  I*  68E4svAT]0N^f  .•IcHO^AR^A’^ilcMA*/*  POINT*  ,T29{  • IINCI  81 
IOENT-SCATTERED) * /*  PHl  THE TA* , T 14 , * PHI -PH I • »T24« •PHI -THETA*  ,T37,  09 

2 •IHFT A-PHl • ,T48i*TH€TA-THETA* I 

66  format  <IHt.9X, 'CONTINUE  EXECUTION  WITH  THE  FOLLOWING  ADDITIONS  AN 
lO/OR  CHANGES*//) 

69  format  l94X.*ELeCTRIC  FIELD  I NTE NSI TV  * / 9X • * DEGREE S* . I IX« • POWE R GA| 

IN* .2BX.*THEfA*.42Xt’PH|*/3X|*THETA*,9X.*PHt •,TX«^THETA*tBX«*PHl*»l 

, 2X.i<8X.*REAL'.Bx.*]MACStX.lMAGN*,9X«*RHASf  *1) 

70  format  <10x,*A«*44jRROR  IN  DATA  CARO  NUMBER  *»I2»*  EXECUTION  STOP 

I PE0**4*4** • ) 

Tt  FORMAT  I40X,*«  A WIRE  SEGMENT  MAVNOT  BE  SHARED  BY  MORE  THAN  FO 

lUR  •*/40X,*R  OIPOLE  MOOES CMECR  DESCRIPTION  DATA  CA 

2R0  ••/40X»*6  EXECUTION  STOMEO 


Tt  FORMAT  (40X.*«  i 

lUR  •*/40X,*R 

2R0  ••/40Xt*6 

3 ••  ) 

72  FORMAT  <40X, •• 

IS  ••/40X,*« 

2R0  ••/40X,*6 

3 ••  I 

73  FORMAT  I30X**A  BACR 
I CALL*//90Xt*RE0yES 

74  FORMAT  ^•l*•T90,JI^  ^ _ 

1 I90»*R  OHIO  iTATE  UN 

2 T90**P  ANTENNA  ANACVSI 


AN  ISOLATED  WIRE  MUST  HAVE  AT  LEAST  TWO  SEGMENT 

AND  Three  points checr  description  data  ca 

EXECUTION  STOPPED 

SCAUERINf  CALL  MUST  BE  INCLUDED  FOR  R BISTATIC 

... 


4 


BBT  AVAIUBIE  COPY 


i 150, •• 
h T$0.»* 
5 T50#’* 


NJOIFItO  FOR  use  AT  ••/ 

NAVil  FUSTGRAOUATE  SCHOOL  •*/ 

3 SEPTEHetR  19T5 

75  rckmil  t M*  t ISO*  791  <9*  l/TSOt  ***  iT7St«*M 

76  fORMAr  ( 750*  *•*  f 1 U|  *AHTENNA*  t T 7a,  •#• } 

77  FCRHAf  I ISO*' •** OX, 'calculations • ,776** •* /T SO, •••  , T 7 8 , • •• /T50 , 291 • 
I**  1 1 

78  FCRHA7  I r50**»* *9X,'N£AR  F If LO* # T 78# • • * > 

79  FORNAT  MSO,***  •9X**FAR  F | f LO*  • T 78  • I 

80  FORMAT  I T50, , 7x,*ttACK SCATTF AinG* tT7B , * I 

81  format  ITSO, • •• •6X,«aiSrAT IC  scattering* *T78* I 
8?  format  I////I 

83  FORMAT  lT30**FREQUENCr  ( I * . T 8 1 .£  1 1 , S I 
86  FORMAT  (TLO.'MME  RADIUS  i nHE  R S M .f  8 1 , 1 1 1 • SI 

85  FORMAT  <T30.*M|ae  CONOUCTIVITV  < Hf dAMHdS/ME TERI ' • 7 8 1 *6 1 1 . 5 1 

86  format  ir«0**w|RE  INSULATED  I NO/ V F 5 1 •, T 8 5 , * VE S * I 

87  FORMAT  ?t5d.*mUl  INSULATED  ( NO/ Vf  $)  • , T 85 . • NQ<  I 

88  format  < 7 )0* ' (NSULArrON  RADIUS  4 Mf TER S I * # 78  1* f It • 51 

89  FORMAT  < 7 30.*  insulation  conductivity  INhOS/ HUE  R I • , T 61 , f 1 1 • 5 1 

90  Format  (Tso** insulation  oielectric  constant  irel at  Wei • . T81, e ii .5» 

91  format  T>0,  ' INSULATI  on  loss  tangent*. rSiiEIKS) 

92  FORMAT  <TiO, 'EXTERIOR  MED  I UM* , T 8 I , • FR^E  SRACE'I 

93  format  <T)0**EXTERIOR  MEDIUM  CONDUCTIVITY  < MHOS/ METE R| • • T8 1 , E 1 1 . 5 I 

94  FORMAT  <T30. 'EXTERIOR  MEDIUM  DIELECTRIC  CONSTANT  1 RELAT I VE ) * , T81 * 

I Ell. 51 

95  FORMAT  C 7 30, 'EXTER lOR  M£0|UM  LOSS  TANGE  NT • ,78 1 ,E U • 5 1 

96  FORMAT  <T50,'M1RE  STRUCTURE *// T 20 ,* SE G* • 4X , 2< • NODE *• i9X ** LOG ATI  ON' 
1.  l8XI/T2i,'N0.'  •3X,i<  • NO.  «.9X.*X*,iSX,*Y*  .liX.'ll'.rxl/ITh,  I2,5X, 
22(  I2,5X,Ell.5.4X.En.5.4X|Ell.8,  IXM  I 

97  format  (T50|*ANT(nNA  F EEDS* / f 40 • <N00E • , I 6X , • VOt T S* /T4& , * NO.* , 12X , 
i 'REAL*,  7X,*  lMAGlNARY*/<r4l,12,6X,214X,Fll.5>U 

V8  format  (T50,'«',  ax, 'OUTPUT  RE OUE STE 0* , T 78 , • A • I 
99  FORMAT  < T 30, * S TR UC TORE  CURRENTS'I 
100  format  1730, 'FAR  F|El05  FO" 


format  1730, 'FAR  FIELDS  FOR  PHI  VARYING  FROM* , IX , F 5, 1 , • TO  *,FS.l, 
1 ‘AND  theta  varying  from  •,F5.i,*  TO  '.FS.!/ 

2150, 'IN  STEPS  OF  '.FS-l,*  OEGREES.*! 

FORMAT  <T30,' BACRSCaTTERING  FOR  PHI  VARYING  FROM  '.FS.I,*  TO  *,F5# 

m'  ANO  TMEfA  VARYING  FROM  'iFS.!,*  TO  •,fI.I/ 

0,'IN  STEPS  OF  •iF5.lj*  gi^EES.'l 
192  format  ir30.*8l$IAflC  SCATTERING  FOR  PNI  VARYING  FROM  *,F»«l,*  TO 
I'lfS.l,'  aW  theta  varying  from  •,F5.I,*  to  *,F5.l/ 

POINTS  IX,V,2r/904T40t3<El 

!85‘fSssll’!i|8::5t^8i  ja  lii  Mit8 

id*  FORMAT  ir$o,*Pt.or  for  tlSrATlC  scaTteAing-for 


^OR  rHETA«*,F9.1l 


107  format 
138  FORMAT 


<T30,*PLUI  FOR  BISTATIC  SCATTERING  FOR  PH1»',F5.II 
iriO.'PlUr  FOR  eACRSCAf TERING  theta**, F5. 11 
1130, 'PLOT  FUR  BACRSCATTERING  PHt*'tF5.ll 
it30,*no  output  ua  Plots  reouested'I 


• • , T T0,'*  » /T 


109  format 
I 10  format 

111  FORMAT 

112  format 

113  format  b 

114  format  I 750,291  ' •* )/T 50,'** , TTPi'R*  I 

I 15  format  110  ------  * - 

10  A (.HANUE 
I 16  format 

I * ' 


ifn 
1T50, 
1 T5J, 


.••»/T50i29<*R'l) 

, • I NPUT  DATA  • ,tT8, • •• I 


3*  I 

CPMAT  <//|OX,h01'R* |/I0X,*THE  DESCRIPTION  AND  THE  GEOMETRY  0 
. STRUCT jRE'/i3a, •must  81  SfATEO  IN  tHi  FIRST  DATA  6L0CK.*/I0X 
2*  EXECUTION  STOPPEO  •••»! 

format  l//|OXt*NO  PART  OF  THE  HIRE  STRUCTURE  CAN  LIE  SEtOM  THE  GRO 
I UNO  Pt ANf.'/lOX, •••••EXECUTION  STOPPEO^^^* * I 


8l6c; 


DOES  NOT  HAVE  A T ERM  IN  AT  ION  C AR 


GEOMETRY  OF  THE 

- 


117 

i UNO  PI  *Nt  .. 

118  FORMAT  I rsc, « STLUCTURE  L0A0S«/T40, 
imaginary*/! T4* 


NODE* •16X.'0HMS*/T41»*N0.* ,12X 


Lu«u 9’  / • ' , iOA  , ' un 

1 .'RE Ai  • . Tx,'  imaginary*/! T4 1 . I 2 . 6X . 214X . E 1 1 .5 1 1 1 

119  FfiPMAT  I tso.  * STAUC 7URE  L QAOS */ 7 39 , * SE GM^ NT  * , 1 4X , • OHMS' / T 4 1 , * NO* , 1 2 
lXi*R< AL*.  rx,*  IMAGINARY' /TT4l.l2,6Xi2l4X.f I 1.5111 

120  FfiPMAf  (t50,'AN7ENNA  F EEOS* / 1 39 , ' SEGME vf * • 1 4X • • VOLTS'/ T41 , 'NO. ' , 1 2 
lX.*RMi*,7X.'  IMAGINARY'/!  741 , 1 2, 6X  . 2 1 4X  . El  I .sll  I 

121  FORMAT  l//T30,*rHE  number  OF  iNfERvUs  FOA  CALCULATING  THE  ElEMENT 
1S'/T30,'IN  ThF  impedance  MATRIX  WITH  SIMPSONS-AULE  INTEGRATION  IS* 
2,/T30it3|*.  IF  CLOSED  FOAM  |NT|GRAT|5n  IS  REOUlREO  SET  |NT*0*///l 

122  format  CUOf'GADUNO  PLANE  I NO/vl  Si  • , 785  , *N0' I 

111  'f§kii  i{i8:;8so  s\i;icmrmifisf';xiJUw|...T.i.,ii.5/ 

i ria.'caguw  coNoucrfvIrr 
, FgSMit  HIO,»643UI«)  .L.Mf  •MK.Edr' I 


. IMAl 
l***ll 
END 


8?? 
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BLNK 


PURPOSE;  to  conpress  data  to  the  left  by  removal  of 

the  blank  spaces  on  the  input  data  cards. 

METHOD:  A (I)  character  is  compared  to  the  blank;  and 

if  it  is  true,  the  A(I*1)  character  is  shifted  to  tne  A(I) 
posit  ion. 


CALLED  BY;  READ 


CALLS  TO:  NONE 


tNf  BlNr  I Al 
OiM^NSION  A(80) 

OATA  SlANK/*  V 
A • 0 

00  1 I •1*80 
J • l-A 
A( Jl  • A<  II 

1 \f  IaII l•fO•8LANK)  K-R^l 

If  U.eO.Ot  AfTUAN 
AI6l*KI  • ftlANK 
Rf  TU«N 
€N0 


10 


i! 


14 
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CBES 


PURPOSE:  to  calculate  the  quantity  B01  where 

B01  = J (z)  / J (z)  . 
o 1 

METHOD;  If  the  absolute  value  of  the  argunent  for  the 

Bessel  function  is  less  than  12,  B01  is  calculated  via  the 
power  series  expansion  for  the  Bessel  function  in  the  DO 
LOOP  ending  at  statement  3.  If  greater  than  12,  the 
asymptotic  expression  is  utilized  at  statement  4.  if  the 
magnitude  of  the  complex  part  of  the  argument  for  the  Bessel 
function  is  greater  than  20,  B01  is  set  to  (0.,-1).  If  the 
complex  part  of  the  argument  is  negative,  the  sign  of  E01  is 
changed  prior  to  returning  to  the  calling  program. 

CALLED  BY;  SCANT 

CALLS  TO;  NONE 


c 


c 


sue^nur  c»f  s i / lAon 
tf  M . . St  I < 

CCHPlf I W I If  RM Jt  r Iff MN.N£2^  t 

data  pi  /3. 

KAffSWI.Ol.W.OI  GU  TO  h 
fact  Off  • 0.0 

m (Q.tO.  I 
• -0. ?!»•<•/ 

UffMJ  e < 1 .0*0.0) 


00  3 Nr  *1 .2 

N • NP»] 

JNINPI  • TfONJ 
« • 0 

1 Pt  • 

TfflMj  • TFt<*i|**<22A/ri.OATII|ff|N»T«)) 
IPtlhP)  • JNI  NPi  » T ) RNi 
IF  iNp.Nf. It  00  rn  2 
FACT'tff  > F AC  TO*  • ) .0/f  I OATIMl 
TfMMN  • TiPHNmT  f ff  MJ»F  AC  TOR 

2 IffHOH  • CAB^TTFHMJt 

If  iiTtffuTt.oT.i.oF-ioi  CO  ro  1 
J ffffNJ  • O.S*2 


BOI  • jN(ll/JN<2) 

A|  rUTTN 

h V • AINACin 

IF  U6SIT)  .0T.20.I  CO  TO  5 
AffC  • I .0* I . )•/ 

II  • CIXPIAAO) 

cc  • (xti  ./ex 

cs  • i.o.-UMiex-j./ixi 

SOI  • ICSffCCI/iCS'CCI 

•ttUAN 

S SOI  • 1,0. 't.l 

IF  lT,tf.6,l  BOI  • 1.0*1. I 

•ITUSM 

fNO 


J 

) 

A 

5 

6 
1 
B 
9 
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DSHELL 


PURPOSE:  to  calculate  the  mutual  impedance  term 

contributed  by  the  dielectric  insulation  on  the  surface  of  a 
thin  wire. 


METHOD;  The  contribution  to  the  impedance  matrix  is 

calculated  utilizing  the  equation  below 


z 

mn 


(6  - &)  ln(b/a) 

in  ’ jw&S  “ 

2 


/ 

m,  n 


P*  (1)  F*  (1)  dl  , 
m n 


where  z is  defined  in  subroutine  SCANT,  & is  the 
mn  2 

dielectric  constant  of  the  insulation,  b is  the  outer  radius 

of  the  insulation,  a is  the  inner  radius,  & is  dielectric 

constant  of  the  external  medium,  and  F is  the  sinusoidal 

expansion  function. 

CALLED  BY:  SCANT 

CALLS  TO;  NONE 


3 


EQUAL 


PURPOSE;  to  determine  position  (location)  of  the  equal 

symbol  on  input  data  card. 

HETHOD:  The  character  search  begins 

passed  to  the  subroutine.  On  returning  to 
program,  the  argument  passed  is  the  column 
equal  symbol. 

CALLED  BY:  READ 

CALLS  TO:  NONE 


in  the  column 
the  calling 
following  the 


EXPJ 


POBPOSE:  to  calculate  the  exponential  integral  with 

complex  limits. 

METHOD:  The  exponential  integral  is  defined  as: 

y -V 

W12  = / dv  = E (VI)  - E (V2)  j2nTi  , 

vi 

where  the  integration  path  is  the  straight  line  from  VI  to 
V2  on  the  complex  v plane  and 

E^(z)  = 

The  integration  path  is  a horizontal  line  in  the  w plane 
or  an  inclined  straight  line  from  VI  to  V2  the  v plane.  The 
integer  n is  zero  unless  this  path  intersects  the  negative 
real  v axis  at  a point  between  VI  and  V2.  When  there  is 
such  an  intersection, 

a)  n = 1 if  Im  (VI)  > Im (V2) 

b)  n = - 1 if  Im  (VI)  < Im  (V2) . 

The  term  j2nTi  is  caluclated  below  statement  12. 

CALLED  BY:  GGHH 

CALLS  TO:  NONE 
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BEST  AVAIUBIE  COPY 


c 


c 


c 


c 


^ChAOUtlNt  MTJ  IVitV?.WW) 

COHPI  E K r C .t  1 St  I , Vt  I V I ,V^  .M  1 2*  ; 

Ot*<FSSl(<N  V(/n.  Ml^n.  Utl6lt  FElbl 

DAI*  V/0.?^^H«i6irEaO.O.II«891>lK)l,U.  ?99?  M6JfOliO.SnS»‘.)^f01,0.9 
E6Wv6UEOlta.l‘>982dr<»E^O^lJ.9iJJr0liF«Ol  t0.<>9t!69l  UC  00.  0 . I / ( ^5  9^  0 

2 1 .0.;^^99•f9‘>E  01 . 0.  )66  76^2  /FOl  . 0 . S9? S ) 396E 0 1 ,0. 7S6S9  1 F 0 I « 0 . 1 QUO? 
>^eCO^tO.  1 3nO?3^t  o;  tO.  16bS9<t06F02.0.?0  7 76S7  9€02*0.2S6?3R99FO?«0.  >1 
9<.0/S  i9F  0 2*0.  )bS  3068iF  02  1 0.  <.80260e6F  02  / 

()A1  A W/0.<*^H96<»60F  ')0.0.<»1  700093FOOtO.  11337!  3flfOO  .0.  t0!99|9  7f-0t,0. 

1 2e  101  72  0E  -0  J,0.S98S<>  KlE*06t0.21823b87E00,0«  3<»22  1 01  7i00.0.26  302  Ks 
2F  00.0.  12b<«2SH?F  00.0.  «i020b66^E-0l  .0.05636  7787-02.0. 12  129  3b  IE -02. 0.  1 

3 I I6  79  90f  -0  »,0.6‘»5  9926  7E  -05 , 0. 2 2 263 I69F  - Ofc  ,0 .9 22  79  3097 -06 . 0. 39 2 1 09  7 
93F-I0,w.|9565l52€-I2.0.l9a)0270e-I5,0.l60059996- 19/ 

OAT*  0/0. 2299509  27  02.0. K9 I 1560  7 02,-0.91  93 15/67  03  ,-0. 7 0 7S93  39F 02.0 
1.  11259  /997  02,0.  l6O2l76UO3,-O.2  3 06M95eO3f-O.5oO9  960/FO3.-O.6  098>0 

2 5 97  02.0. 1 2 259  7707  02  . - 0.  10 1 6 1 9767 02 . -0. 9 7 2 1 9 59 17 01  . 0 . 79 7 2 9601 E 01 . -0 
3.210695  7«.t0  2,0.22C96  990fOi  ,0.09  72  0299701/ 

data  7/O.2nO31O>EO2.-O.3M5V/8  7FO3i-O.9  79092  2O7O2iO.l29OO6  7 27O)iO 
1 . 1 7999  2 267  , - 0. 129109317  03  .-0.  55  7055  797  03.0.  1352  90017  02  .0 . 1 969672 

21  70  3.0. I 79995 2 07  02, -0.32981 6197  00, 0.3 10 200^6702 ,0.0 165 765 770 1 .0. 2 2 
3|3696|702.0. 391 29 0427 02 , 0. 0 1 6 16 799e 0 1 / 

00  12  J|H«1,2 
* • «7At< II 
V • AlWAGl/l 
7 15  • 1 .0, .01 
A0  • C AHSi 2 1 
If  <48.70.0.1  GO  70  11 
IF  <X.^I.O..A9O.A0.GT.1O.)  go  TO  10 
t*  • ABSlTl 

IF  < X .1 7.0. . AND. VA.GT .10. 1 GO  TO  10 
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GANT1 


PU&POSE:  to  consider  the  wire  structure  as  a 

transmitting  antenna  and  calculate  the  input  impedance  and 
current  distribution. 

METHOD:  If  a wire  antenna  is  driven  by  a voltage 

generator  v located  at  one  of  the  current  sampling  points 
i 

1 and  if  displacement  currents  are  neglected.  Ampere's  law 
i 

yields 

V = v F (1  ) 
a i m i 

wnere  F is  the  sinusoidal  expansion  function.  Thus,  the 

excitation  voltages  V will  vanish  everywhere  except  where 

m 

V IS  not  zero, 
i 

The  DO  LOOP  ending  with  statement  2 uses  the  delta-gap 
model  defined  above  to  determine  the  excitation  voltage 
CJ(I)  for  all  the  dipole  modes.  These  are  stored 
temporarily  in  CG(I).  Then  subroutine  SQHOT  is  called  to 
oDtain  a solution  of  the  simultaneous  linear  equations. 

SJROT  stores  the  solution  (the  loop  currents)  in  CJ  (I)  . 

In  the  DO  LOOP  endinj  at  statement  6 , the  complex  power 
input  and  input  im pedaiice  (s)  are  calculated.  The 
time-average  power  input  (PIN)  is  the  real  part  of  the 
complex  power  input. 

Subroutine  RITE  is  called  to  make  the  transformation 
from  the  loop  currents  to  the  branch  currents.  If  IHR  is  a 
positive  integer,  RITE  will  write  out  the  list  of  branch 
currents. 
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Finally,  GANT1  calculates  the  radiation  efficiency  by 
calling  subrouinte  GDISS  to  obtain  the  time-average  power 
dissipated  in  the  lumped  loads  and  the  imperfectly 
conducting  wire. 

CALLED  BY:  BAIN 


CALLS  TO: 


GDISS 


SOROT 


tcukiex  rv  , ccfw 



COMMON  IML 
00  ) ■•liN 

CJII  » • I .0* .01 

K • JAUI 


00  2 KK-1,2 
HA  « UIk) 

KB  > IBIKI 

a : t. 

ij  rJ.is.-jf 

1 If  iKA.CQ.n 
JJ  • K«NM 
CJU)  • CJil 


CJf 

60.131111 


CJU)  • CJII)»FI«V6I  JJ) 
2 K • JBU  I 


3 C0HTINU6 


00  4 J-l.N 
4 CCdl  • CJtn 

call  SOI^OT  (C,CJt0.112.NI 
112  • 2 
Vll  • r.Ot.OI 
PllWN  ■ N 

IF  (IGAO.OT.O)  NHM  • (N4|tfl.)/2 

00  A {•UNNH 
HH  m lAUBlIll 

CGEN*CCEN22, 


vv  • CJt  I )4C0«IJC<CGEN1 


MAITE  Utl)  NN,211 
9 Yti  • rii frt 
A COMT iNUt 


IF  I IMF .GT.OI  MAITE  < A;?! 

CALI  RITE  nA.IB.lHM.TwR.lltUiM.fOvNOtNMvCJfCGilGROI 

Gc  • rEacivii) 

If  (igro.gt.oi  gg«2.4Realiviii 
Fin  • GG 

call  GOTsS  I AN.CGwCMNiO.OISStGANvNNrSGOvZLOfiSI 
RETURN 


r format  t50l« 'ANTENNA  BRANCH  CURRENTS*! 

• FORNAT  UOKv'THE  INRUT  INREOAnCE  AT  NODE  •«ll« 


1S*»FI9«?**  ♦ A*» 
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GDISS 


PURPOSE; 
dissipated  in 
lumped  loads. 


to  calculate  the  time-average  power 
the  imperfectly  conducting  wire  and  in  the 


METHOD:  The  time-average  power  dissipated  by  the  wire 
is  calculated  in  the  DO  LOOP  ending  at  statement  1 utilizing 
the  equation  below: 


P = 
d 

wnere  R 

s 

radius  of 


* 

I I dl 


is  the  surface  resistance  of 
the  wire. 


the  wire  and  a is  the 


The  power  dissipa 
by  the  DO  LOOP  ending 
perfectly  conducting, 
by-passed. 


t ■'d  by  the  lumped  loads  is  calculated 
at  statement  3.  If  the  wire  is 
CMM  <0,  the  first  calculation  is 


CALLED  BY:  GANT1 


CALLS  TO: 


NONE 
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SUf*OLUTlNe  (/OISS  (AH,COiCMM,OtOI$S|GA*4<NH|SGO*ZLOfZ$) 

CCMPLFX  CG(  M I SGO<  1 1 »Z(.0<n  «C JAtC li 
DiNfNSION  01  I I 
OAt  A Pi  /3.  1«.IS9/ 

OISS  • .0 

I  r umm.i E .0.1  GO  ro  2 

4LPH  • R( AL (GAHI 

beta  • AIMAOiGAMI 

PH  • PE  Ai.  ( Z S i/<  ^ .•PMAMI 

c 

00  I K*1.NN 
OA  • Oi A I 

OCN  • CAdSISGOlA  I )*»2 
E AO  • E A P( AL  Ph*OA  ) 

CAU  • ( EAOM  ./EAO)/^• 

CHO  - ( OSt  BE  Ta*OK I 
SAU  • Oa 

tE  (ALPH. he. 0.)  SAO«lEAO*t ./EAO)/(?.»ALPHl 
SBD  • OA 

IF  (BFTA.NE.O.i  ^ 80 > S I H ( BE T A* OR W BE T A 
FA  • BH*< SAD*CAO-SeOACBO)/OtN 
FB  • 2. APH*iCAO«SBD-SAO«CBD)/OEN 
C JA  • c G(Ki 
(.  • R»HH 

c je  • co(i I 

1 OtSS  • 0ISS*FA*(CABS(CJA)P«2»CAa$(CJ8)**?)AFBP(PEAL(CjA)PP€AL(CJBI 
1»AIMAGICJA)*A|HAGICJ3I  ) 

c 

c 

2 00  3 J-I.NM 

R • jthM 

3 OISS  • 0ISS*PEALUL01J))A<CABS(CGUnAA2lAPEALIZL0tRnp|CABSICG(Rl 
IM«2> 

C 

PE TuPN 
END 


1 

2 
3 
A 

5 

6 
I 

a 

9 
10 
1 1 

11 

lA 

15 

{? 
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19 

20 
21 
22 
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3A 
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GFf 


PURPOSE:  to  calculate  the  far-zone  field  of  a 

sinusoidal  electric  monopole. 


METHOD:  If  an  electric  line  source  has  length  d and 

endpoints  at  (x  ,y  ,z  ) and  (x  ,y  ,z  ) , then  the  coordinates 
111  222 

of  any  point  on  the  source  are 

X = X ♦1  cosx 
1 


y 


y ♦ 1 cosy 


z = y +1  cosz 
1 

where  cosx,  cosy,  cosz  are  the  direction  cosines  of  the  1 
axis,  and  1 is  the  distance  along  the  source  measured  from 
the  endpoint  (x^,y^,z^).  Let  the  current  distribution  on 

the  monopole  be 


I 

1(1)  = 

where  I and  I 


sinh  >'(d  - 1)  + sinh  yi 
sinh  /d 

are  the  endpoint  currents. 


1 2 

field  of  this  source  is 

E = (cosx  cos6  cos^  - cosy  cos9  sin4) 

4 

E = (-cosx  sin#  ♦ cosy  cos#)  E 


The  far-zone 


- cosz  sin9) 


E 

1 


where 


ge  ygi  /ft » > 

E = - - [ (e  - g sinh  yd  - cosh  yd)  I e 

1 4nr  (1-g*)  sinh  yd  ‘ 

yd  yf<*) 

♦ (e  ♦ g sinh  yd  - cosh  yd)  ® ] 
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f< » ) 


= X sin6  cos*  ♦ y sin6  sin*  ♦ z cos0 
1 1 1 

f<2)  = X sin0  cos*  y sin0  sin*  + z cos0 
2 2 2 

^ = cosx  sin0  cost  ♦ cosy  sin0  sint  ♦ cosz  cos0 

and  (c,0,t)  are  the  spherical  coordinates  of  the  observation 
point . 

In  this  subroutine  the  range  dependence  has  been 
suppressed.  The  far  field  vanishes  in  the  endfire  direction 
where  GK  = 0.  If  a ground  plane  is  present  (IGRD  > 0)  the 

equation  above  is  decoaposed  into  the  x,  y,  and  z components 

and  the  reflection  coefficients  are  applied  before  E and  E 

6 ^ 

field  components  are  returned  to  the  calling  program. 

CALLED  BY:  GFFLD 

CALLS  TO:  MONE 


4 3 


StSf^OUTlNE  GEF  (X  A.VA.ZA.XB  tVSt  1 
lA.ETl .Ft;,ePl,EF2»lGx6*EAxT 
cNpiEX  |MR.Ay.ftH.KR.EX»EY.EItEE 
CCNPLEX  ETl.EriiEAl.EP2f6AM«EtA 

iSSStli 

coMPttx  csr 

FP  • 12.56632 


A.2AiXB»VBf28*O.C60iS60fCTH»STHfCPH«SPHfGAMtCT 


XAB  • XB-XA 
VAB  ■ TB-YA 
ZAB  • is-ZA 
CA  • XAB/0 
CB  • YAB/0 
C6  • ZAB/O 

C • <CA«CPH*CB*SPH|«STH^CG6CTH 
GA  • 1.  -G*C 
ETL  - ( .0«.0) 

|t|  - j .Of. 01 
|pl  - i .Ot.O) 


ETL  - ( .0«.0) 

|t|  - j .Of. 01 

|pl  - i .Ot.oi 

£P2  • L .0*.0j 

IF  IGK.LT ..OOl 1 GO  TO  3 

FA  • UA«CPH»VA*$PH|#STHa2APCTH 

FB  • lXB«CPH«rB«SPF41*STHAZBPCTH 

EGFA  - CEXP(GAR*FA| 

EGFB  • CexP<  CA«I*Fb1 

EGGO  • CEXPL  GAM*G*6i 

CST  - ETA/<&A*SG0PFPI 

ESA  - CST*eGFA«TEGGO>GPSGO-CGOI 

ESB  • CST^IgfbaII./EGMaCPSGO-CGOI 

IF  < IGAO.LE.OI  go  to  2 


AV  • <-L.f01 


Is 

Py  • '(ERAPCTH-AAI/rERAPCTHAft 
AH  • (C Th-AA 1/(C THaAAI 
I EX  • CAPESA 

h : h:ili 

EE  ■ <EXPSPH*EYPCPH)P|AH-AV| 
EX  - EXPAV»EEPSPH 
Et  • EYPpy-EEPCPH 
E2  • -E2PAV 

ESA  •-EKPCA-EYPCBpEZPCG 
EX  • CAPESB 

\\ ' 

EE  - TexpS^ETPCPHIpIBH-AVI 


EX  • EXPRV»EEPSPH 
EV  • EYPAy-tE*CPH 
li  - -ezpAY 
ESB  •'EXPCA-EYPCB* 
T • LCAPCPHACBpSPH 


2 T - LCAPCPHA 
P • -CAPSPH# 

ETl  • TPCSA 

ihilli! 


f 


GFFLD 


POHPOSE:  to  calculate  the  fac-field  for  the  thin  wire 

structure. 

HETHOD;  The  far-field  for  the  structure  is  calculated 

from  the  loop  currents.  The  loop  currents  are  either  the 
currents  produced  by  the  transmiiting  antenna  calculations 
of  subroutine  GANT1  or  the  currents  produced  by  an  incident 
plane  wave. 

If  the  incident  field  is  generated  by  a distance  source 

with  spherical  coordinates  (r  ,0  ,♦  ) , the  excitation 

0 0 0 

voltages  induced  by  a incident  plane  wave  are 


m 


where 


E = E exp(K  r-f  ) 
i 0 0 

where  E is  a vector  constant,  ? is  a vector  from  the 
0 0 

coordinate  origin  to  the  distance  source,  and  r is  the 
radial  vector  from  the  origin  to  the  observation  point. 


The  field  E is  generated  by  test  dipole  m when 

n 

radiating  in  the  homogeneous  medium.  Osing  the  vector 
potential,  the  field  at  the  distance  point 


E 

m 


4nr 


-/r 

e 


0 


0 


exp(K  r-r^) 


dl 
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where  the  radial  coaponent  is  to  be  suppressed.  From  the 
above  equations, 

Unr  yr 

V = - fl-  e 0 e e 

a j«u  0 a 

If  an  antenna  gain  calculation  is  desired,  INC  is  set  to 
zero.  PH  and  TH  denote  the  spherical  coordinate  direction 
of  the  distance  observation  point.  The  phi-polarized  (EPPS) 
and  the  theta-polarized  (ETTS)  components  of  the  electric 
field  intensity  are  returned  to  the  calling  program. 

If  INC  = 1,  a backscatter ing  calculation  Is  desired.  In 
this  case  PH  and  TH  denotes  the  incident  angles  for  the 
incident  plane  wave.  These  are  also  the  spherical 
coordinates  of  the  distance  source.  The  outputs  returned  to 
the  calling  program  include  absorption,  extinction,  and 
scattering  cross  section  for  each  polarization;  scattered 
electric  field;  and  echo  areas. 

If  INC  = 2,  a bistatic  calculation  is  desired.  In  this 
case  PH  and  TH  denote  the  spherical  coordinate  of  a distance 
observer.  Since  this  calculation  uses  the  induced  loop 
currents  (EP  and  ET)  , a backscattering  call  must  proceed 
this  calculation.  The  outputs  returned  to  the  calling 
program  consist  of  the  scattered  electric  field  components 
and  echo  areas. 

EPP(I)  and  ETT  (I)  denote  the  phi-polarized  and 
theta-polarized  far-zone  fields  of  dipole  mode  I with  unit 
terminal  current.  In  a backscatter ing  situation,  the 
excitation  voltages  EP(I)  and  ET  (I)  are  obtained  by 
multiplying  EPP  and  3TT  by  the  constant  CJI.  Then  calls  are 
made  to  SQROT  which  stores  the  solution  (the  induced  loop 
currents)  in  EP(I)  and  £1(1).  RITE  is  called  for  the  branch 
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currents  CG  (J) , and  GDISS  is  called  for  the  tine-average 
power  dissipated  in  the  imperfectly  conducting  wire  and  the 
lumped  loads.  This  power  is  denoted  PDISS  and  TDISS  for 
phi-polarized  and  theta-polarized  incident  waves, 
respectively . 

In  scattering  problens,  the  incident  plane  wave  has  unit 
electric  field  intensity  at  the  origin.  GGG  denotes  the 
time-average  power  density  of  the  incident  wave  at  the 
origin.  ACSP  and  ACST  denote  the  absorption  cross  sections 
for  the  phi  and  theta  polarizations. 

PIN  and  TIN  denote  the  time-average  power  input  to  the 
wire  structure,  delivered  by  the  equivalent  voltage 
generators  VP  and  VT  at  the  terminals.  PIN  and  TIN  apply  for 
the  phi  and  theta  polarizations,  respectivity.  The 
time-average  power  input  is  regarded  as  the  sum  of  the 
time-average  power  dissipated  and  the  time-average  power 
radiated  or  scattered  by  the  wire.  ECSP  and  ECST  denote  the 
extinction  cross  sections  and  SCSP  and  SCST  denote  the 
scattering  cross  sections. 

The  distance  field  is  calculated  in  the  DO  LOOP  ending 
with  statement  7 for  scattering  situations,  and  in  the  DO 
LOOP  ending  with  statement  9 for  the  antenna  situation. 

The  radar  cross  sections  (echo  areas)  SPPfl,  SPTM,  STP^I, 
and  STTH,  are  defined  as 


a = limit  u 
r 

where  S and  S denote  the 
s i 

the  scattered  and  incident 


2 2ar 

nr  e S / S 
s 

time-average  po 
fields  evaluate 


i 

wer  densities  in 
d at  the  origin. 


For  an  antenna,  the  following  definition  is  employed  for 
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the  power  gains: 


2 2ar 

G {^»^)  = limit  4nr  e S(r,9,^)  / P 
p r --  i 

where  P , GG,  denote  the  time-average  power  input  and 
i 

S(r,6,^)  is  the  time-average  power  density  in  the  radiated 
field.  GPP  and  GTT  denote  the  power  gains  associated  with 
the  phi-polarized  and  the  theta-polarized  components  of  the 
field,  respectively. 

The  use  of  the  variables  JPLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 

CALLED  BY:  MAIN 

CALLS  TO:  GDISS 

GFF 


RITE 


SQROT 
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GGHH 


PURPOSE:  to  calucate  the  mutual  impedance  between  two 

filamentary  monopoles  with  sinusoidal  current  distribution, 

i 

METHOD:  As  stated  in  subroutine  SGANT,  the  mutual 

impedance  of  coupled  dipoles  may  be  expressed  as  sum  of  four 
monopole-monopole  impedances.  This  subroutine  calculates  the 
mutual  impedance  with  closed-form  expressions  in  terms  of 
exponential  integrals. 


For  skew  monopoles  it  can  be  shown  that  the 
monopo le- oionopole  mutual  impedance  is  given  by: 


t -z  2 

i+j  n mm 

z = (-1)  Bfe  (F  -e  G +0  G) 
ij  j1  12  22 


-t  -2  Z 

n ai  IS 

-e  (F  -e  G ^e  G)] 
j2  11  21 


where  m=  2/i,  n = 2/j  and 


16  TT  sinh  d sinh  d 
1 2 


The  functions  F are  defined  by; 
ik 


qz  cos  |[ 
i 

F =2  sinh  d e E (R  ♦ qz  cos  | 

ik  » i i 


- qt) 


where 

being 


k 

q = (-1)  » 

considered. 


d and  d are  the  lengths  of  the  monopoles 

12 

The  functions  G are  defined  as  follows: 
ik 
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I 

I 


G = E (B  ♦ 

<32  _ 

♦ q't  - jq”) 

♦ E(R 

♦ qz 

♦ q't 

♦ jq") 

ik  2 

2 

2 

2 

- E (r  ♦ 

qz 

♦ q't  - jq”) 

- E (R 

♦ qz 

♦ q't 

♦ jq") 

1 

1 

1 

1 

i l( 

where  q = (-1)  » q * = (-1)  t and  q”  = qb  ♦ q'c  with 
D = c cos  f and  c = d/sin  f.  The  angle  f is  the  angle  formed 
by  the  apparent  intersection  of  the  two  monopoles.  This 
will  be  discussed  later  in  detail. 


ik 

In  the  above  equation  for  G , t denotes  the  position  of 

an  observation  point  somewhere  on  raonopole  2.  B and  a are 

12 

the  distances  from  the  endpoints  of  monopole  1 to  this 
observation  point.  Finally,  the  £ functions  are  defined  as 
follows: 


E(a  ♦ jq") 


dw 


where  1 and  q"  are  real  quantities  with  dimensions  of 

length,  a is  a function  of  t,  a = a (t  ) , a = a(t  ) and 

1 12  2 

/ = jw  /u6.  The  integral  above  is  evaluated  by  subroutine 
EXPJ. 


PARALLEL  SKEW 


To  explain  the  input  data  for  GGHN,  refer  to  the  above 
figure.  If  the  monopoles  are  parallel,  then  the  new 
coordinate  system  is  defined  such  that  the  new  z axis  is 
parallel  to  the  monopoles.  The  coordinate  origin  may  be 
selected  arbitrarily.  SI  and  S2  denote  the  z coordinates  of 
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the  endpoints  of  the  test  aonopoie,  T1  and  T2  are  the 
coordinates  of  the  endpoints  of  the  expansion  aonopoie,  and 
D is  the  perpendicular  distance  (displ aceaent)  between  the 
monopoles.  The  mutual  impedance  of  parallel  monopoles  is 
calculates  in  the  last  part  of  GGHM  below  statement  5. 

For  skew  monopoles,  let  the  test  monopole  s lie  in  the 
xy  plane  and  the  expansion  aonopoie  t in  the  plane  z = D.  (D 
is  the  perpendicular  distance  between  the  parallel  planes.) 
If  the  monopoles  are  viewed  along  a line  of  sight  parallel 
with  tne  z axis,  the  extended  axes  of  the  two  monopoles  will 
appear  to  intersect  at  a point  on  the  xy  plane.  Let  s 
measure  the  distance  along  the  axis  of  the  test  monopole 
with  the  origin  at  the  apparent  intersection.  SI  and  S2 
denote  the  s coordinates  of  the  endpoints  of  the  test 
monopole.  Similarily,  let  t measure  the  distance  along  the 
axis  of  the  expansion  monopole  with  the  origin  at  the 
apparent  intersection.  T1  and  T2  denote  the  t coordinates 

of  the  endpoints  of  the  expansion  monopole.  Let  s snd  f be 

unit  vectors  parallel  with  the  positive  s and  t axes, 

respectively.  Then  CPSI  = = cos  . The  monopole  lengths 

are  d and  d . 
s t 

The  output  data  from  GGMM  are  the  impedances  P11,  P12, 
P21,  and  P22.  In  defining  these  impedances,  the  reference 
direction  is  from  SI  to  S2  for  the  current  on  monopole  s, 
and  from  T1  to  T2  for  the  current  on  monopole  t.  In  the 
impedance  P , the  first  subscript  is  1 or  2 if  the  test 

ij 

dipole  has  terminals  at  S1  or  S2  on  monopole  s.  The  second 

subscript  is  1 or  2 if  the  expansion  dipole  has  terminals  at 

T1  or  T2  on  monopole  t.  The  monopoie  lengths  d and  d are 

s t 

assumed  positive  in  defining  the  input  data  CGDS,  SGD1  and 
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SGD2 


For  parallel  raonopoles,  CPSI  = 1 or  -1.  SI,  S2,  T1,  and 
T2  are  cartesian  coordinates  for  parallel  monopoles  and 
spherical  coordinates  for  skew  monopoles.  For  skew 
monopoles,  the  radial  coordinates  Si,  S2,  Tl,  and  T2  tend  to 
infinity  as  the  angle  p tends  to  zero  or  n.  Therefore,  if 
the  monopoles  are  within  4.5®  of  being  parallel,  they  are 
approximated  by  parallel  dipoles. 

CALLED  BY:  GGS 

SGANT 

CALLS  TO;  EXPJ 
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GGS 


PURPOSE:  to  calculate  the  mutual  impedances  between 

two  filamentary  monopoles  with  sinusoidal  current 
distributions. 

(1ETH0D:  The  monopole- monopole  mutual  impedance  as 

defined  by  SCANT  is  calculated  using  the  equations  defined 
in  subroutine  GNF.  The  endpoints  of  the  axial  test  monopole 
s are  (XA,YA,ZA)  and  (XB,YB,Z3),  and  the  endpoints  of  the 
expansion  monopole  t are  (X1,Y1,Z1)  and  (X2,Y2,Z2) . DS  and 
DT  denote  the  lengths  of  monopoles  s and  t,  respectively, 
CAS,  CBS  and  CGS  are  the  direction  cosines  of  monopole  s, 
and  CA,  CB  and  CG  are  the  direction  cosines  of  monopole  t. 

The  effects  of  ground  for  vertical  co-linear  monopoles 
are  applied  in  a slightly  different  manner  than  mentioned 
previously.  As  with  self  impedance  calculations,  the  test 
monopole  and  the  expansion  monopole  are  laterally  displaced 
by  the  wire  radius.  This  lateral  displacement  is  used  to 
determine  the  angle  of  incident.  This  technique  is  applied 
at  statement  8. 

If  INT  = 0,  GGS  calls  GGMM  for  the  closed  form  impedance 
calculations.  Otherwise  GGS  calculates  the  mutual  impedance 
via  Simpson ' s-rule  integration  with  the  following  number  of 
sample  points:  IP  = INT  1.  If  the  monopoles  are  parallel 

with  small  displacement,  GGS  calls  GGMM  to  avoid  the 
difficulties  of  numerical  integration. 

Since  the  point  {X,Y,Z)  of  subroutine  GNF  lies  on  the 
expansion  monopole  t,  T is  the  integration  variable  and  is 
measured  from  (X1,Y1,Z1).  Cl  is  the  current  at  T for  the 
mode  with  terminals  at  (X1,Y1,Z1),  and  C2  is  the  current  at 
T for  the  mode  with  terminals  at  (X2,Y2,Z2).  C denotes  the 
Simpson  * s-rule  weighting  coefficient. 
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Below  stdtenent  7,  GGS  performs  some  analytic  geometry 
in  preparation  for  calling  GGMM.  The  remainder  of  this 
section  is  concerned  with  this  preparation. 

Let  s denote  a unit  vector  in  the  direction  from 

(XA,YA,ZA)  toward  (XB,YB,ZB).  Also  let  I denote  a unit 

vector  from  (X1,yi,Z1)  toward  (X2,Y2,Z2).  Then  s«E  = 

cos  6 = CC  where  0 is  the  angle  formed  by  the  axes  of  the 

two  monopoles.  Let  monopole  s lie  in  one  plane  P and 

s 

raonopole  t lie  in  another  parallel  plane  P . CAD,  CBD  and 

CGD  are  the  direction  cosines  of  the  unit  vector  5=^x5/ 
sin  y which  is  perpendicular  to  both  planes.  To  obtain  the 

distance  DK  between  the  two  planes,  a vector  I?  is 

11 

constructed  from  (XA,YA,ZA)  to  (X1,Y1,Z1)  and  take  DK  = 

17  .3. 

1 1 

A line  is  constructed  from  (X1,Y1,Z1)  to  the  test 
monopole,  such  that  the  line  is  perpendicular  to  the  test 
monopole.  SZ  denotes  the  s coordinate  of  the  intersection  of 
this  line  with  the  test  monopole,  and  the  cartesian 
coordinates  of  this  intersection  are  XZ,  YZ,  and  ZZ.  The 

direction  cosines  of  s x 3 are  CAP,  CBP,  and  CGP. 

From  the  point  (X1,Y1,Z1)  in  plane  P , a line  is 
constructed  perpendicular  to  the  point  (XP 1 , YP 1 , ZP 1 ) in  the 


5b 


plane  P . This  line  is  parallel  with  3 and 
s 

Let  1?  represent  a vector  from  (XZ,YZ,ZZ)  to 

PI  denotes  T?  (s  x 3)  . SI  and  T1  are  defined 
GGNN. 

CALLED  BY:  SGA NT 
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XA)*(XB*XA)*( VB-YA)AIYB- 


YB'YAj/AG 
TIRCAS 

tNcbs 

XX«$PH-CYY»CPH|«|RH1-AY1I 
XX*RV|»Ef •SFM 
eVY*Avl-E£PCPM 
-Et  IPKVIACGS 

'EX1PCAS-EyI*C0S»E21«CGS 


:?? : Vriliu 

EE  • (EXX*SPH- 


EYVPCPH|P|KH2*AV2) 


• EXX*AV2*EE*SPH 

• Evy*sv2-EE*CPM 

--Exi*CA^-E?2*CBS*E22*CCS 

'h1Sjli2S;|Sf^?|r' 

• C* CEGO'l./EGOI/2. 

• CEXPIGA**«T» 

- C*(EGCH1.^EG0I/2. 

• PIUETIPCI 

• PU*£n*C2 

: yAilWin 
Mhh 

• ^SGN 


ill : H 

Pf TUAN 

sh 


-ETA*0ELT/(3.*FPPSG0S*S&0T) 
CSMPi  J 


T*P12 

t*p5‘ 

T#P2 


OAl 

‘ih 


SZ l ‘Oixtc 

iO*r  U «]i-«»-JZZ«C*S»»»J*tT2-TA-$I2»C»SI»»2»«  t2-I*-Sll*C6il*» 

I 0«1>0K2 1/2. 
f IDQD.Gr.2>).*AN.ANO.INT.er.OI  «o  ro  t 

ALl  CCnh' .S?2>OOOtC«>S>SGOSi  SCOT , 1. ,CT* ,6*N,»1| ,f l2tP2l 


,L 

‘li'ilaS;li!:ll  !S’V5". 


!8 : liHSBilSilSi 

A • (X1*XA>*CA0AIYI*YA 
A • ABSIOAI 
F IOA.tftA9)  OA  • AN 

• kaasIkas 

• TA«s1*CBS 


PCBOPI21-2AIPCGO 
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GNF 


PURPOSE:  to  calculate  the  near-zone  electric  field  of 

a sinusoidal  electric  oionopole. 


E 

e 


METHOD:  An  electric  line  source  is  located  on  the  z 

axis  Kith  endpoints  at  z and  z as  shown  in  the  above 

12 

figure.  Let  the  electric  monopole  have  the  following  current 
distribution: 

I sinh  /(d  - 1)  I sinh  yi 

1 (1)  = _ . 2 

sinh  yd 


where  I and  I are  the  endpoint  currents,  y is  the  complex 
1 2 


propagation 

source  lengt 
E(e)  = 0 and 

E(e)  = 


constant  of  the  medium,  d = z 
h.  The  cylindrical  components 


4n6  sinh  yd 


-yR 

[ (I,  e t 


I 

2 


2 


e 


- z is  the 


of  the  field  are 


-yp 

2)  sinh  yd 


-yR 

♦ (I  cosh  yd  - I ) e ‘ cos  0 

' 1 2 I 

-yR 

4-  (I  cosh  yd  - I ) e 2 cos  6 ] 

'2  I 2 
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-/B 

[ (I  - I cosh  yd)  e 

‘ 2 ' ' H 

2 

-yR 

(I  - I cosh  /d)  e 1 ] 

2 1 R 

‘ 1 

where  is  the  intrinsic  impedance  of  the  medium  and  where 
denote  the  cylindrical  coordinates  in  a coordinate 

system  centered  at  the  endpoint  of  z^. 

These  expressions  exclude  the  field  contributions  from 
the  point  changes  at  the  endpoints  of  the  line  source,  since 
these  charges  disappear  when  two  monopoles  are  connected  to 
fora  a dipole. 

Let  the  coordinate  s measure  distance  along  the  test 
monopole  with  the  origin  at  (XA,YA,ZA) . From  any  point 
X,Y,Z,  a line  is  constructed  perpendicular  to  the  aonopole. 
S2,  denotes  the  s coordinate  of  the  intersection  of  this  line 
with  the  nonopole.  The  length  of  the  line  is  the  radial 

2 

coordinate  6,  and  RS  denote  6 . R1  and  R2  are  the  distances 
from  (XA,YA,ZA)  and  (XB,yB,ZB)  to  the  point  (X,Y,Z). 

In  the  statements  above  statement  1,  the  above  equations 
are  solved;  and  after  statement  1,  the  cartesian  components 

(E  ,E  ,E  ) of  the  field  are  determined.  If  a ground  plane  is 
X y z 

present  (IGFD>0)  the  reflection  coefficients  are  applied  to 
the  cartesian  components  before  returning  to  the  calling 
program. 

CALLED  BY:  GNFLD 

CALLS  TO;  NONE 


E(z)  = ^ 

4ti  sinh  yd 
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%UhMOUT  INe  GNF  I XA.VA,  ;*|XB.YBW8.XfY*Z  I AM,  OS  f CGOS  i SGO  S • E T A , GAM , E X 
tGAO|EAfli 

CCMPLtX  FKf-'.UVI  .AHl 


C(  **PUX  L JA,(  JB,  E J1  Ml  ,ER2  »£Sl 

CUxPtiX  tXl«EYi,E£l,fcX2,EV2fE£2 
DATA  PJ  ' JhJ  !>^/ 

CAS 
C8S 


>SG0S.GAM,CST,CG0S,ETA 


il} 

ttl 


Yri 
III 
AS 
Al 
I JA 

IJI  • £ JA/PI 
B7  • SGftT(SSfi£7**7J 


I JB 
E J2 
ESi 
E S7 
( 0 I 
EP2 
A»S 


C E X P( 'GAM*  R2  t 
t JP/«2 

t J7-tJl*CGOS 
E Jl-f J2*C&DS 
j,.o] 


I .0, 

( .0.  .0) 

Am*am 

1 »•  4Q  S.  I T .AHS)  GO  TO  I 

crni  • ^n/«i 

crH7  • 

ERl  • lEjA»SOOS*EJA*CGOS*CTHl-EJb*CTH2|/RS 
t-E Jb*SOOS*fc JB*CG0S*CTh2-£ JA*CTHI)/BS 
f TA/U.*P|  •SCl>S» 
cst*asi*cAS*ER  i*xxn 
CS«*irSl*CBS*f  Al*vvn 
CSr*(fcSl*CGS*EA 


ih2 
I CST 

E XI 

E rl 
ill 
EX2 
EY2 


CSf*(f S7*CA$*€R7*XXZ} 

. _ C$T*itS2*CeS»tR2*YYi 
fl2  • CSI*IES2*C0S*EA2*ZZ2I 
IT  (IGXO.LE.Oi  XETJAN 
R Vi  • I • 1 . ,0) 

RHl  > 4 '*1  • tO) 

RV2  • 

RM2  • 4 - I . tOl 

IT  < ICAO.eO.I I GO  TO  2 

Ml  • SORT! 4 XA-X » • ( XA-Xi ♦! YA-V) • I YA-Y I » 
R2  • SQRT4UB'X|*IXd‘Xi*(YB-y|*IYB-Yn 
THI  • ATANlRl/l2A-ii) 


I xP  XA ) /OS 

7 

4 YB-YA  I/OS 

0 

\ ZB-ZAI/OS 

9 

|X-XAI*CAS»I V'YAI ACBS*  U-2AI*C&S 

10 

si-  DS 

li 

XXA-S7*CAS 

13 

Y-Y4-$z*ces 

16 

Z-ZA  sz^c^s 

15 

XXZ**2*yyz«*Z*ZZZ**2 

SgRMaS*ZZi*«2) 

17 

Cexp<-CAN*R|  I 

10 

S9 
0 
I 

\\ 

It 

26 

27 

20 

29 

30 

l\ 

33 

36 

35 

36 

37 

30 

39 

60 

n 

63 

66 

45 

66 

67 

60 


*TAMP2/iZ8-in 


CSWRTUU- 

csqrturr 


NITMU* 


-liN4TH 


-4f RM'COSI Tm| 

(C0S(Tm1|-<R- 
-iE«R*lOS4  7H 
ICOS4TH2I-RR,. 

$gMT(ixA-x0i4ixA-xe)*4 vA-ve)RiYA-yeii 

8 

T (RG.IT • l.E*32 I GO  TO  3 
|X0>XA>/RG 
iVB-YAl/R^ 

lEXlASPH-EY  l*CPH)«4IIHl*ftVn 
-EXI*RVIaE|6SPH 
-fr|PA Vl-EE6CPH 

'hlMlir.iR!'*'*”'*'"'*’ 
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GNFLD 


PURPOSE:  to  calculate  the  near-zone  electric  field 

intensity  at  a given  point. 

METHOD:  This  subroutine  calls  GNF  for  the  near-zone 

field  of  each  wire  segaent,  and  suns  over  all  segnents  to 
obtain  the  near-zone  field  of  the  wire  antenna.  FI  is  used 
in  a Banner  similiar  to  FI  of  subroutine  SGANT.  CJ(I)  is  the 
loop  currents  calculated  by  subroutine  GANTI. 

The  use  of  the  variables  JFLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 

CALLED  BY:  MAIN 

CALLS  TO:  GFF 


vfrltCi 1,1x2 gam 

COMPLEX  era 

12111.  Mill.  Dili.  XIII.  VIII.  Ill 

DIMENSION  H0IINN.4I.  NOI  t I 
OXTX  PI.TP/).I*M9.6.2«}1I/ 

I?:  1:§’:8[ 

n • I .OfO) 

00  2 K-L,NM 

KA  • lAlRl 

Re  • leiKi 

NMO  ■ IGRO 

If  IR.LE.NM/21  ICAO— I 

MOR  • NOIRI 


I'l  !'laii!!!|iinr'-' 


iniifiiiiiir' 


2 ciMTliWC 

5JTUNN 


000 

oco 

000 

OOOA 

8882 

883.' 

0009 

881? 

ooli 

0013 

OOIA 

ools 

88lt 


0020 
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LEFT 


PURPOSE:  to  determine  position  (location)  of  the  left 

paren  symbol  on  the  input  data  card. 

METHOD;  The  character  search  begins  in  the  column 

passed  to  the  subroutine.  On  returning  to  the  calling 
program  the  argument  passed  is  the  column  following  the  left 
paren  symbol. 

CALLED  BY:  READ 

CALLS  TO:  NONE 


c 


c 


St^MOUl  INC  ItF  T 
CLHMcN  /4/  AltOI 
04r*  PLf^T/M»/ 

K • N 


INI 


00  I l•Kf•0 

7f*U1  I I.CQ.^IE^T  1 GO  TO  I 
1 CONTlPMJf 

N • 1 
(NO 


> 

4 

5 
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LIMECK 


PURPOSE:  to  insert  grid  characters  on  the  polar  plot. 

METHOD:  The  peroid  character  (ISM (2))  is  inserted  in 

the  proper  position  in  the  statements  above  statement  4.  In 
the  statements  after  statement  4,  the  grid  numbers  labels 
are  inserted  on  the  horizontal  axis. 

CALLED  BY:  POLPLOT 


CALLS  TO:  NONE 


c 

c 

c 


c 


c 


SUfiROUTINe  LINECK  IX, Y| 

THIS  subroutine  insures  all  grid  charactors  lie  on  the  polar  grid 

CCMN(,N  ISTN.LINE 
INTEGER  Y 

dimension  ISYMIIAI.  LINEUSOI 
IE  lY.EU.OI  GO  to  ) 

K • 0 

IE  (X.l r.lO.OI  GO  TO  5 


SET  UP  AREAS  OE 
I 


"PERIOD"  POLAR  GRID  POINT  CHARACTERS 


I NT! XI 
IABS(  I I 
I • ABSIX) 

IE  I W-  1 1 .GT  .0.5  I I •! 

1 IE  (I.IT.10.0.0R.2,GT.111.0I  GO  TO  2 
L INI  4 IT  • ISvMf  2 I 

LINE<6UI  • ISVM13) 

LINE46^I  • ISVMlil 
X • R*| 

IE  IKaO.2)  GO  TO  2 
I • 12?-l 
GO  TO  I 

2 L INE  4611  • I STM4  2 I 
IE  4V.NE.0I  CO  TO  5 

) on  •*  K-  1 1 * U 1 

LINE4XI  - ISYNI21 
A CONTINUE 


I 

? 

3 

A 

5 

6 
7 
6 
9 

10 
1 1 
12 
n 

lA 

I 7 
18 

19 

20 

li 

2i 

2A 

25 

26 
27 

l§ 

30 

31 

32 

33 


EILL  IN  GRID  NUMBER  LABELS  ON  HORIZONTAL  AXIS 

3A 

39 

i INE  4 11  1 

• ISVH47I 

36 

L|NE120j 

• 1 SYMl 10) 

37 

1 INE 1 21  1 

• ISYM49) 

36 

LINEI2?| 

• ISYMil 

» 

39 

1 |NE  4 30  ) 

• ISYNI9 

AO 

1 INE  31  1 

• 1 SVMI 5 

A1 

1 INE  I 32  1 

• 1 SYMl  1 1 ) 

A2 

L INE  AO) 

• ISYN40) 

A3 

1 INtlAl  I 

• ISYM43 

AA 

L1NE(A2I 

• 1 SYMl  1 

) 

AS 

IINEI90) 

• ISYMI 71 

A6 

L 1 NE  1 9 1 1 

• ISVMI9) 

A7 

L1NEI92I 

• 1 SYMl  1 1 ) 

AB 

L INE  161  I 

• ISYM41 

A9 

LINEI70I 

• 1 SVM4  7 

90 

L INC  4 71  1 

■ ISYM43) 

91 

L INC  4 72  ) 

- 1 SYM(  1 

) 

92 

L INE  4 801 

• ISYMIB 

S3 

L 1NEI81  ) 

• ISYM49 

9A 

LINEI82i 

• 1 SYM4  1 

1 

99 

LINE  4 901 

• jSYM49 

96 

LINE19I) 

• l$YM<9 

97 

L|NE492> 

■ ISYM4  1 

) 

98 

LINECIOO) 

• ISYMI 

01 

99 

LINE  4 ion 

• ISymH 

) 

60 

limiioH 

• ISYMI 

l» 

LINEIltil 

• ISYMI 

7) 

62 

CONTINUE 

63 

RETURN 

6A 

END 

69 
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NUMB 


( 

PURPOSE:  to  place  degree  numbers  on  the  polar  plot. 


METHOD;  The  current  line  which  is  being  printed  is 

passed  to  the  subroutine  in  the  calling  argument.  If  this 


line  contains  degree  numbers,  these  numbers  are  placed  in 
the  correct  position  by  the  IF  statements. 

CALLED  BY:  PTPLOT 

CALLS  TO:  NONE 


suBPuuriNe  NUMB  iyi 

This  suriKouriNE  purs  DE&«ee  numbers  on  polap  grid 

CCMNf'N  ISYH.LINE 
INffCER  Y 

OIHENSION  iSYMiUI,  efNE((30l 
if  «v.NK17l  GO  TO  I 
L INE i • I SVMl  7 ) 

L INE  I 3<»l  • I SYHl  91 
IINM  ibl  • ISVM<6I 
L INE I 0 7 ) • I SyHT  6 » 

L INE I 00 ) > i $V«1  \ i I 
L IN( T 091  • I SVM! 6 7 
if  GO  TO  Z 

I tNf) • istmi  r» 

L INt  11)1  • I SVM(  1 U 

t INM  19  I • I SYM»  t>  I 
I iNf ( JOB  » • I SYMl 6 I 
I INtl  1091  • |SYM( 9) 

I iNt I I I Oi  • i SYMl 6 i 
lY.NE.O)  GO  TO  3 
LINE! n • ISYHl 71 

I NE10}  • isrmi)) 

I INE  T 9i  • I SVMl 0 ) 

LINE  1 M 31  • t SYMI 6 1 


L I NE  U I 9) 
L I NE  It  I 9 7 


ISYMI 6) 
1 SYM167 


1^  IY.NE•*^1 ) GO  rO  9 
L INE  I • J SYMl  0 ) 

I INEt  I )l  • ISVMl  7) 

L INEl  19»  • |SYMl6i 
I INM  i08)  • I SYMl  9) 

I INf 1 1091  • ISYMl 9l 

LINE! 110)  • |SYM(6I 
17  lY,NE.-i;(  GO  fO  5 
I (NT  L 3 ) ) • I SYMl 8) 

I (NE i )9 ) • ISYMl l6) 

L |NH  )S  I • ISYMl  6 I 
llNEltt?)  • ]SVMi9| 
LlNf«88f  • iSYMiB) 

L INE I 89l  • I SYMI 6) 
CONT INUt 
RE  TURN 
END 


NUMBER 


PURPOSE:  to  convert  alpha-numeric  numbers  to  floating 

or  fixed  point  numbers. 


METHOD:  After  initially  determining  the  sign  of  the 

number,  the  DO  LOOP  ending  at  statement  6 scans  each 
character  beginning  at  N1.  The  DO  LOOP  ending  at  statement 
3 terminates  the  outer  DO  LOOP  if  the  character  being 
compared  is  not  an  alpha-numeric  number.  The  DO  LOOP  ending 
at  statement  5 converts  the  alpha-numeric  number  to  an 
actual  number.  Below  statement  7,  the  multiplier  correction 
is  applied  to  the  floating  point  number  before  returning  to 
the  calling  program. 


CALLED  BY:  READ 


CALLS  TO:  NONE 


c 


c 


c 

c 


c 


c 


SUnROUT(N€  NUMBER  <NliN2»X«!XI 
CO‘<N(.n  /a/  Aieoi 
0I^( NS lUN  H( 101 
OAU  8/  '0  • , * I * , 

DATA  ANNOStPLUS 
DATA  AK,AN,AU/»K* ,*«*,*U*/ 

N • Nl 
NSI&N  • 0 
1 I - *1 
IX  • 0 
ISM  - 0 

If  (A(N).EU«PLUSI  N-N*l 
If  (AINJ.Nfc. ANNUS!  CO  TO  t 
NSiON  • 1 
N • N»l 


1 00  6 I 'Nf 80 

If  (All I.Nf. POINT}  60  70  2 
I Sf  T • I 

GO  fC  6 

^ If  itsM.f 0.11  It  • n»i 

00  ) It  10 

If  (A< I ).60.e(Rl I CO  TO  4 

3 CONT |NU€ 

CO  TO  7 

4 00  S 8«ltl0 
KR  . K- I 

If  t A<  I I .EO.eiR)  I NONB*KK 
9 CONT INUE 

lx  • NJH8a10«1X 

N^  . IM 

4 CONTINUE 

7 If  (NSION. f 0.1  I I X • •IX 
r • IX 

If  I i I.LT  .01  II  • 0 
X • Y/(  10P*I I I 
If  (A(N2} .eo. POINT!  N^•N2♦1 
If  (AInIi.EO.AR!  X • X*!000. 

if  (aIn}) .EO.AH)  X • X*0.001 

if  (ainIi.eu.au  X • XAO.oooool 

If  l(A(Nn.EO.AX!.OR.U|N2!.IQ.AM|.OR.(A(N2!.EO.AU!  ! N2«N2»I 
N|  • N2 


1 

2 

3 

4 

5 

6 
J 


9 

to 

I I 
12 

13 

14 

!• 

18 

19 

20 
21 
22 

23 

24 

25 
?t 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 
43 

46 

47 

48 


f NO 
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POLPRT 


PURPOSE:  to  control  the  plotting  of  the  polar  plot, 

METHOD:  This  subroutine  is  the  main  subroutine  in  the 

polar  plot  package  and  is  responsible  for  calling  the 
various  subroutines  of  the  package. 

The  scale  factor,  S,  must  be  changed  according  to  the 
printer  character ist ics.  The  scale  factor  in  this  subroutine 
is  set  for  ten,  10,  characters  per  inch  for  the  abscissa  and 
eight,  8,  characters  per  inch  for  the  ordinate  axis. 
Therefore  S = 10. /8. 

After  initializing  DATAX,  DATAY,  and  X,  the  input  data, 
Y,  is  scanned  to  determine  the  normalizing  factor.  If  this 
normalizing  factor  is  less  than  1.E-32,  an  error  statement 
is  printed  and  the  plotting  is  abortted. 

In  the  DO  LOOP  ending  with  statement  8,  each  line  of 
the  polar  plot  is  printed  after  a call  is  made  to  PTPLOT  to 
establish  the  ploar  grid  information.  The  variable,  DIM,  is 
used  to  as  a scaling  factor  for  the  polar  plot.  The  value  of 
1,0  will  cause  all  of  input  data  to  be  plotted,  however,  if 
only  the  values  less  than  one-half  of  the  norualizing  factor 
are  of  interest,  then  DIM  can  be  set  to  .5.  This  will 
enlarge  of  the  center  of  the  polar  plot. 

CALLED  BY:  MAIN 

CALLS  TO;  PTPLOT 

SART 
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SUBHOUflNf  POlPftT  (N4NF,YI 
Cu**Ht'N  ISY»l,4.|Nt 

OlPfNSION  X{)60l.  r()60i,  0ATAIO60),  0ATATn60lf  llNE(l>Olt  ISVMI 
I I ^ I 

OIM(N$ION  TITlAUl,  T|ri2l2) 

DAf A T 1 Tl A/ •PHI  • , • THE  T»  A 
H • 360 

OIM  • 1.0 

NSf  • I 
AST  - t 

S IS  EACTOA  Of  PRiNTEAt 

ABSCISSA  CHAR.  PER  INCH  / ORDINATE  CHAR.  PER  INCH 

S • 10.0/6.0 

ZERO  OATAX  AND  OATAY 

00  1 IA*1,N 
0 • lA-  1 

DATA  X(  Ul  • 0.0 
DATA  YU  A I • 0.0 
1 XIIA)  * 0*3. 1A159?7/180.0 


FACTOR  IS  THE  NORNAC U 1 NG  DIVISOR 


00  I ia-?,n 

2 IF  (FACTOR.LT.YdAII  FACT0R*Y1IAI 

IF  INANE. EO. II  T I TLl-TI TLAI n 
IF  INANE. Eg, T|TL|«T|TLAI2I 

IF  I INANE  .Ey.  )|  .OR.  INANE.  £0.6  I .OR.  inane.  EO.  71. OR.  INANE.  E0.6  M 7 I 71. 
12  I II  I Tc A| I I 

IF  I INANE .EO.S I .OR. 1NANE.EQ.6I .0R.INANE.E0.9I .OR. I NAME .EQ. 10  I I TIT 
lL2lil«T ITLAI2I 

IF  (INANE .EU. 3 I .OR. inane. EQ.51 .OR. inane .EQ. 71 .OR. I NANE.EO. 9 1 1 71 7C 
WUj-TI  IL  All  I 

IF  1 1 NAME.  EQ  .6  I .OR.  I NANE.EO  .61  . OR  • INANE.  EO  • .OR  . I NAME  • EQ.  10  I I T|T 
U2I?)«T  ITLAl?) 

IF  IFAC TOR.Or. I .6-32)  GO  TO  3 
IF  INANE. LE. 21  MRITE  l6t9)  TlTCl 
IF  INANE. GE. 31  WRITE  UtlOl  TITL2 


IF  INANE 
RETURN 


NCRNACI ZE  DATA  TO  ONE 

3 on  6 lA'l.N 
6 Yl  lAI  • Vl  lAl/FACTOR 

IF  INANE. IE. 21  WRITE  16,111  TITLlffACTOR 

IF  1 1 SANE.  GF  .3  I . AND.  I NANE  .CE.b  I 1 WRITE  16,131  Tni2,FACT0R 
IF  INANE. GE.ri  WRITE  16,121  TfTL2, FACTOR 
Flic  OATAX  AND  OATAY  ARRAY  FROM  X AND  T ARRAY 

DO  3 |A»1  .N 

OATa  XIIaI  > V(  lAMCOSiXl  lAI  I 
5 OATA  YllAl*  Yl  1 A IPSINIXI 1 A1 1 

SORT  OATA  BY  QNOlNATE  NAGNITUOE 

call  sart  ioatax,oatay,ni 

OATAX  ANO  OATAY  ARE  SORTED  BY  OESENOING  MAGNITUDE  ON  ThE  OAFAY  VAC 

SET  UP  FUR  Plotting  polar  grid  with  data 

DO  B IVV-1  ,61 
CALL  PTPlOT  I IYY.SI 

LINE  IS  RETURNED  WITH  POLAR  GRID  INFORMATION 

SET  UP  ’Y*  BIN  SIZE  UPPER  ANO  LOWER  LIMITS 
ULL  IS  the  lower  bin  LIMIT 

UL  IS  the  upper  bin  limi t 

BIN  • OIM/BO.O 

ULL  • 01M'12*tYY*l IPSIN 

UL  • ULL»2*B|N 


cycle  Through  data  to  fino  which  ones  fall  :n  *v*  bins 


IF  INST.GT.NI  go  70  7 
00  6 JJ«NST,N 


69 


\4 


I 

c 


c 

c 

c 

c 

c 

c 


IF  lOAVAirl  JJI.ir.ULL)  GO  TO  7 
• JJ 

A*«AC  • soft  T(0ATA)(<  JJ)«OATAXI  JJI  ♦OATAYIJ  JF^OATAYI  JJIJ 

CneCK  that  magnitude  is  nut  over  dim 

IF  (AMAG.GT.01M)  CO  TO  6 

OK  IS  THE  FINAL  LINE  POSITION  FOR  TmE 

Ok  • DATAXI JJ)«S«A0«0/0IM»61.0 
IF  IJK . i T . 10.0)  GO  TO  A 
K • I NT  I Ok  I 
K • lABSlK) 

OK  • AHSlOK) 

IF  UOK>K  ) .GT  .0.  S I K«K*1 

IF  (OK. I T . lO.O.OR .OK.CT .1 11 .0)  GO  TO  6 

LINEIK)  • ISVM(A) 

6 CONTInJE 

T CONTINUE 
NST  • KST»l 

PRINT  OUT  ONE  LINE  OF  PLOT 

mRITE  <6, IN)  LINE 
B CONT INUE 

RE  TURN 


•9  format  I)0X,1AN.»  CQMPPNFNT  of  THE  ELECTRIC  FIELD  IS  lESSVJOX, 

1 'THAN  l.E'BNf  therefore  THIS  FIELD  MAS  NOT  • /lOI • *PL0T TEO.  tXEC 
Motion  4|ll  continue  as  normal.*//) 

10  format  (10X,*rHE  maximum  value  of  the  bistatic  pattern  FOR  •/ 

1 lOX. IAN. •'* , 1AN» • ( INCIOENT-SCATTEREO)  IS  LESS  THAN  */ 

2 lOX,  • l.E*30.)  POLAR  PLOT  NOT  CALLED.*///) 

11  format  1*1*. IAN.*  electric  field  antenna  PATTERN  FOR  SPECIFIED  PL  A 
lNE.*,9X.*N0AMAllnNG  FACTOR*  *•610. S) 

12  FCRmAT  (•laiSTATlC  SCATTERING  PATTERN  F OR  * . t AN . • - * , I AN , * ( I NC I DENT- 
ISCATTEREQ)  polarization.* .9X.*N0RMALI2INC  FACI0R«*,E10«5l 

13  FORMAT  F UACKSCATTER  IMG  PATfERM  FOR  * > 1 A N « ' • 1 AN  • M (MC  rOEMT- SC  A T T 
lEREOI  P0LARUAtl0M.*t9X,*M0RMALUING  F AC  TOR  ■' f E 1 0 . SI 

IN  FORMAT  (IX.130A1  ) 

EMO 


>>? 

1)1 
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! 


I 


PTPLOT 


PURPOSE:  to  establish  the  grid  information  for  the 

polar  plot. 

METHOD:  In  the  DO  LOOP  ending  at  statement  1 the 

alpha-numeric  characters  are  transferred  to  ISYN  in  order 
to  pass  via  COMMON  to  other  subroutines.  In  the  statements 
following  statement  2,  the  equations  for  the  plotted 
concentric  circles  are  established.  Below  statement  7 the 
grid  marks  on  the  090-270  axis  are  inserted. 

CALLED  BY:  POLPRT 

CALLS  TO:  LINECK 

NUMB 
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I 

C 

c 


c 


c 


c 

c 

c 


c 


SUAAUUtiNf  »>TOLUV  ilVY,S) 

THIS  StlffROuriNf  SfTS  UP  POLAA  GPIO  INFORMATION 
CCMMiN  ISVH-iiNfc 

01  Mt  NSinN  L iNE  4 t 301  • ISVHU^K  ISYNtt^l 

0*1  A I SVN/  |H» . 1H«  , IH  . IH* ,1H/, 1 HOf  IHI •1H2, lM3tlHA  » 1 H5 • I H6 t I H0 , 1H9/ 
iNffcGIA  VtVVfM 

SET  UP  iSVM  FROM  ISYN  FUR  COMMON 

00  1 K*  t . 

ISYMIKI  • ISYN(M 
1 CCNTlNUe 


CICAA  LINE  AND  SET  TO  BIANK 


OC  2 I • 1 . I 30 
2 L INI ( n • ISYMI 1 ) 


1 

2 

3 

A 

5 

6 
t 
6 
9 

10 

IJ 

I 3 

lA 

U 

19 

20 

i\ 

a 


IF  (r«E0.01  GO  TO  7 

25 

26 

SET  UP  EUUATIONS  FUR  CONCENTRIC 

CIRCLES 

27 

28 

Yt  • Y*v 

29 

2 • IY/*2.5/2MS 

30 

* • 6|.0»S0«n  2500.0-2l 

31 

( ALl  L I NICA  1 » , Y 1 

32 

IF  lY.  1.1.32.  OR. Y,lT. -321 

GO 

TO 

3 

n 

M • 6l.0»SUPn  1600.0'2I 

3A 

CALI  LINECA  4X,Y) 

35 

3 

IF  ( V .GT  .2A.0R. t .L T 2A| 

GO 

TO 

A 

36 

X • 6l . 0»S0Rf  4900 .0*2 1 

37 

CALL  LINICa  (X.Y) 

36 

A 

IF  4V.GT.L6.0P.Y.LT.-16I 

CO 

TO 

5 

39 

X • 61 .0*  SUMT I AOO.O-I F 

AO 

call  L 1 NICa  I K . V I 

Al 

5 

IF  lY.GT.e.PR.Y.LT.-B)  GG 

TO 

6 

a2 

A • 61 .0»S0RT1 100-2 1 

A3 

CALL  LINICa  U*Y} 

AA 

SET  UP  EUUAflONS  FOR  MULTIPLES 

OF  30  DEGREES 

A5 

6 

X - 6I.0*1.732051*Y*S 

A6 

CALL  LtNECK  4X,Y1 

A7 

X • bl .O^YAS/l. 732051 

A8 

I 

c 

c 

c 


1 C*IL  LJNEC*  (X.v) 

PUT  IN  PUIAR  PLOT  NUMBER  LABELS 

CALL  M)MP  ivi 
W • lAOST  Yl 


FILL  IN  POiAM  PLOT  AT  OOOt  090*  180*  ANO  270 


IF  (M.NE. 
I INEISSI 
I I N(  ( 5 7 ) 

L INC (S9) 

L INF  4 63  ) 

I INt  4 6S  I 
L INF  i 67  I 
IT  IM.NE, 
I INI  I 56  I 
L INI  I • B I 
L INF  16 J i 
L INI  I 6?  ) 

I INI  t6«»l 
I INL 1 66  I 
IF  (W.NE. 
I iNi  (57  1 
L |7«  159  1 
L INI  160  I 
L INI (62  I 
LlNf  (631 
I INE  465  I 


AUl  00  rc 

• ISYMI2I 

• ISVH(2i 

• SVM(21 

• I SYM4  2 I 

• ISY«(2| 

• ISTNI 21 

321  GO  rc 

• ISVA(2j 
- ISVM(2I 

• ISrH42| 

• ISVNljj 

• ISYN(2I 

• fSYN(21 
2a1  go  T( 
■ I SVMI 2 I 

• I SVM( 2 I 

: insiil 

• ISyN(?' 

• isYNi; 


10 

IF  IM.NE.I6 

1 GU  T 

i |NE( 58  1 • 

i SVMI  2 

L1NE4601  • 

is»"(  e 

L1NE<62J  • 

ISVN(2 

L1NEI6A1  • 

1 STM( 2 

11 

IF  IW.NE.Ofl 

1 GO  1 

LINE! 591  • 

jSVMIJ 

1 1 NE  1 6 3 1 • 

1 SVMl} 

12 

CONT INUE 

RE  TURN 

END 

8 


9 


10 


11 


12 


A9 

50 

51 

52 

53 
5A 

55 

56 

57 

58 

59 

60 

63 

6A 

65 

66 

67 

68 

69 

70 

73 

7a 

75 

76 

77 

78 

79 

80 
81 
82 
63 
BA 

85 

86 

87 

88 
89 
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READ 


PURPOSE:  to  interpret  and  translate  the  input  data 

cards. 

METHOD:  The  program  utilizes  free  format  for  the  data 

cards,  that  is,  the  program  uses  character  recognition  to 
determine  which  parameters  are  being  read.  In  the  IF 
statements  containing  A(1),  A (2),  A (3)  , and  A(4),  the  first 
four  characters  on  the  data  card  are  compared  to  the  first 
four  letters  of  the  key  words.  This  will  determine  the  type 
of  parameters  that  card  contains.  The  other  IF  statements 
determine  which  parameters  are  being  read. 

Subroutine  BLNK  is  called  to  remove  the  blank  spaces  on 
the  parameter  cards.  Subroutines  EQUAL  and  LEFT  are  called 
to  determine  the  position  of  the  equal  character  and  the 
left  paren,  respectively.  Subroutine  NUMBER  is  called  to 
convert  the  alpha-numeric  characters  to  numbers,  either 
fixed  or  floating  point.  This  numerical  value  is  assigned 
to  the  parameter  just  determined. 

A detailed  explanation  of  the  data  cards  is  found  in 
appendix  II  titled  "USERS  MANUAL". 

CALLED  BY:  BLNK 

EQUAL 
LEFT 
NUMBER 
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StBROUT 1 NE  fftn  I I A. 1 B, IBISCi 1C  ARQ, 1 GA I N , I GROi INF  A Rt INT , I SC A T . I MR • 
I 1 FI  AGf  Kf  lAb.KC.F  ^,1  OAO  fl  ZO.MSGvNBAP.NBI  P iNFf  P.NGENtNHtNP  , ABAP  , ABAT  , 
?AFFP.AFFt.ABlP(AR|l , AH , 6M, C MM • E P 2 1 1 R Jtf P A • F mC ■ HO f i PHAF  • PHA I iPHlF .P 
)Hi I , PhSF .pHSI  . tHAF . fHAl.TMIF, fHl I . ThSF, T hS 1 , $ I G2 • S I C )• i I OA • T 02 » t6 ) 


trH>|  , inar  |iHai*inir|ini  i* 

A^VUt  I ,ifXNP,ir,YNP.Z.2LL0,ZNPtSTtPl 


[ HHOH  /a/  At  4o» 

C; “Pi fcx  voi n n , 21 1 01 
OlHt  NSIUN  ' " 

Dl Hf  NS  1 ON 

DATA 

|Ab,Aa^'a't'P'.' 

DATA  PL AN< .COMMA ,M|NUStPLEFT«PUINT.R IGHT • SLANl/*  * t * t * » * ‘ * • * 1 * 


MtNSIUN  Udlf  IBtM.  Xilli  riUt  Kilt  KGENlllt  XFLAGU) 

MfNSION  XNPlUt  YNPtlti  2nP(II.  LZOtll 

TA  AA.AB,AC.AOtAe,AFtA/^tAH,AItARtAltAMA.AN,AOtAP,AQtAP|AStAr,ALlt 
,Ai/*A',*B*.*f*.*0*.*f*.*F*,*&*t*M*»M»t»R*,*L*t*M',*N*t*0*#*P*t 


I . * I ' , •/  •/ 

Pad  • S J.29^ / 19 


IN! 
IBISL 
IGAIN 
INFAR 
I SCAT 
( «R 


- I 
•I 


- I 


If  t|FlAG.f0.6l  GO  TO  2 
If  IMSG.NE.OI  GO  TO  A 

1 READ  <S|7A.EN0>12I  A 

2 IF  K AM  I .NE.AC  I tUR.  t A(2KNE.BI  ANKI.OR.I  A(  ) ) • NE  • BL  ANK  I . OR.  I A I Al  . NE 
t . dl  ANR 1)  GO  TO  3 

MRITE  l6t7Al  A 
GO  TO  I 

3 MRl TE  16. 7S) 

GO  TO  5 » 

A RE  AO  1 St  76fEN0>  72  I A 
5 KARO  • ICAOO*l 

^RME  t6»77l  ICARU.A 

IF  l|MSG.NF.OI.AN0.11A(n.CQtAE).ANO.lA(2}tEO.AN|.ANO.(A(3I.EO.AOI 
III  GO  TO  70 

IF  M MSC.NF .01 . AND. 1 1 A( 1 I .EQ.AS I .AND.! A( 2 l.EQ.AT 1 .ANO. 1 At  3) .EQ.AOI 
I . AND.  ( A t<.)  . Eg.APli  I CO  TO  A9 

IF  (1  AM  1 . EQ.  AC  I .ANO.  I A<  21  • EQ.  BL  ANK  I . AND  • 1 A 1 i I • EO.BL  ANK  | .ANO.  (A(  Al 
I.IG.BLANKI I GO  TO  73 

IF  imsg.gt.oi  go  to  a 

CALL  sink  1A) 

N • A 

INSULAT ION 

IF  1 1 At  M .NE.AI I .OR. 1 A121 .NE.AN) .OR.l Al 3) .NE.ASI .OR. lAI A) .NE  tAUI  I 
1G(J  TO  10 
KFlAGUOl  • 1 


call  left  ini 

t If  .OR.  IA(N. 2I.NE.ro  I .OR.  I RIN»  )l  .NE 

I.AIII  GO  TO  7 
RFIAGM)  • I 

call  equal  ini 

Cali  number  lN,N2tXi,|x| 

pM  • I i 

IF  U|N2l,Eg.RlGHTl  CO  TO  A 
IF  lA(N2I.Nf  .SlANTI  go  to  71 
N • N2*  1 
GO  TO  6 

' MAIINI  .NE.AOl  tOR.lAINFlI.NEtAII.OR.UINFjl.NE.AEl.OR.UINASI.NE 
I.AIII  GO  TO  6 
RFiAOlA)  • I 
call  EOUAl  INI 
call  number  IN.N? .XI , IX ) 

ER2  • XI 

IF  IAIN2)  .EQ.R  IGHTI  GO  TQ  A 
IF  IAIN2I.NE.SLANTI  CO  TO  71 
N • N2* I 
GO  TO  6 

8 IF  MAINI.Nt.ACI.0R.IAlN»||.NE.A0l.0RtIAINF2I.NE.AN1.0R.(AIN*3).NE 
1 .AOI I GO  TO  9 
XFLAGISI  • I 
call  E..JAL  INI 
call  NUMBER  |N.N2tXltlX) 

SIG2  • XI 

IF  IAIN2I  .EQ.R  IGHT  I GO  TO  A 
IF  fA|N2I.NE.SLANT|  GO  TO  71 
N • N2»l 
GO  TO  6 

9^|F^MA|N|.nE.ALI  .OR.l  AINFII.NC.  AOI. OR.l  A (N*2I  .NE.ASI  .OR.l  AINF  31. NE 

RFiAGlM  • I 

CALI  equal  ini 

call  number  IN,N2tXltlXI 

T02  • XI  ' 

F IAIN2l.fQ.R|GHl|  GO  TO  A 
IF  IAInK.NE. SLANTI  GO  TO  71 
N • N2» I 
GO  TO  6 


WIRE 
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3 
A 
5 
t 
7 
B 
9 
10 
1 1 
12 
13 
lA 

ii 

16 

19 

20 
21 
22 
23 
2A 

25 

26 

n 

29 

30 

n 

33 

3A 

35 

36 

37 

38 

39 
AO 
Al 
A2 
A3 

AA 

A5 

A6 

A7 

AB 


A9 

50 

51 

52 

53 
5A 

55 

56 

57 
5B 

59 

60 
61 
62 
63 
6A 

65 

66 
67 
6B 

69 

70 

?i 

71 
7A 

75 

76 

77 
7B 
79 

50 

51 

B3 

BA 

85 

66 

B7 

SI 

90 

93 

9A 

ll 


•V 


c 


c 


c 

c 

c 


c 


c 


< A(2).NE.A|) .OR*!*! 3>.N€.AA). 

IAIN«11 .NE*AA).0ft.lA<N*2I.NE. 

• IX) 

E8  TO 


I IX) 


OR.IAUl.NE.AE)) 

AO>.OA«UIN*)|,HE 


10  IF  n At  1 ) .HF .AW) .OR. 
l&O  TO  i i 

CAIL  LFfT  IN) 

11  If  U A| N) .NF . ARI .OR* 
l.AD)  00  10  12 

• KfLA0(2)  • 1 
CALI  rguAL  (N) 
call  NOMPFR  (NtN2,Xl 

AM  • N I 

If  U(N2)  .(^.RISHT) 

IF  (AIN?). NF. SLANT) 

N ■ N2* I 
CO  TO  1 1 

12  If  ( I A(N) .NF.AC) .OR.! A(N»1I .NE. AO) •OR.  I A ( N»2 ) .NE • AN) .OR. I A IN » 3) .NE 
l.AO))  00  TO  Tl 

KflAOn)  • 1 
call  FoUAL  IN) 
call  NJMHFR  (N*N2(X1« 

CMM  - X I 

If  IA|N2).E0. RIGHT) 

If  IA(N2I. NF. SLANT ) 

N • N2* 1 
CO  TO  I 1 

FXTFRNAt  MEDIUM 

1 ) If  ((All  I.NF.AC ) .OR 
IGO  TO  1 7 
XFlAG(e)  ■ I 
call  LFfT  IN) 

14  If  1 1 A INI .NF .AC  I .OR.! AINrI) .NE .AO) .OR. I A (N^2I .NE. AN) .OR • I A I N»3) .NE 
1 .AO) ) GO  TO  1^ 

KfLAGI9)  • 1 
CALL  fCUAL  INI 
CALL  NUMBER  IN,N2»Xli 
SIG}  ■ XI 

IF  I AIN2)  .FO.R IGHT I 
If  IA(N2).NC. SLANT) 

N • N2«  1 
00  TO  14 


CO  TO  A 
GO  TO  71 


I Al 2) .NE.AX ) .OR.I Al )).NE.AT) .OR.IAIAI.NE.AE) I 


> IX) 


GO  TO  A 
00  TO  71 


15  If  (lAINl.NE.ADI .OR 
1 .AL) ) GO  TO  16 
RFLAGIIO)  • I 
CALL  EQUAL  IN) 


>( AIN»1) .NE. Al ) .OR. I A|N#2) .NE. AE) .OR. I A|N«)).NE 


97 

98 

99 
1 00 
101 
102 
103 
iOA 

105 

106 
107 
iOB 
109 
I 10 
1 1 1 
112 

1 14 

115 

1 16 
1 1 7 
118 
1 19 
1 20 
121 
122 
123 

125 
I 26 

P' 
128 
129 
I 30 

131 

132 

133 
I 34 
135 

36 

37 
36 

39 

40 

41 

42 

43 
144 


C 


call  number  (N,N2tXl,lX) 

ER  3 • XI 

If  <A(N2).FQ. RIGHT)  GO  TO  4 
If  IA|N2).NF. SLANT)  GO  TO  71 
N • N2» I 
GO  TO  14 

^‘i1as)!*oo'to^7i‘‘**°**'**‘'^*‘**^^'*°**°'*’***^*^**^^*‘'*’®*’**‘***’^***^ 

Kf  I AGU  I ) • 1 

call  equal  IN) 

call  NUMBM  IN,N2»X1«1X) 

T03  • XI 

If  ( A| N2) .FO.R IGHT ) GO  TO  4 
If  IA|N2).Ne.SLANT  GO  TO  71 
N • N2* 1 
GO  TO  14 


LOAD 

”l&0  }i*|  •*0'*0»*  I.O*.  (*(♦!. NE.  *01 » 

XflAGI 14)  • 1 
GO  TO  19 

‘*1  GO*  TO  22“^^'*** 

NfLAGI24)  • 1 


iAU  Vf  I 


LFfT  IN) 

number  IN.N2|Xl.|X) 

LEOUJ-i'iS'  '' 

N • N2»  1 

CALL  NUMBER  (NiN2*XlilX) 
RMAG  • XI 
N • N2*  1 

call  number  «N»N2«XUIX) 

‘SISr  • *1 


lEAL  • RMAG*COS(ROE 


G/RAO) 

G/radI 


R 

R _ 

RIMAC  • RMAGRS lNiR5E£ 

/ILOII)  • CNRLXIRREAL«R|MAC) 

RIGHT  ) CO  TO  4 
UANT)  GO  TO  71 


?A|n||«E9. 
jr  IA|N2).NE. 
I • 1*1 

N • N2*  1 

GO  TO  20 


145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

158 

159 

160 
161 
162 
163 

64 

|65 

66 

67 

68 

69 

70 

i! 

74 

75 

78 

79 
I BO 

ll 

B5 

RB 

B7 

SB 

B9 

90 

91 
192 


75 


Zl  KflAQtZ^I  • -1 
lOAO  • -I 
OU  TO  4 


ZZ  If  I I Al  I J .NE.Af ) .OR.l A«2) .NE.AA) .OR.I A( )) ,Nf .AE I .OR. U(4) •NE .AQi I 

i(*o  fu  zi 

Af l A0( 11*1  , 

t Al  I t M f IN) 

CALL  NlJMbER  |N,H2,Xl,  |X) 
f Nf.  • t\ 

GO  TO  A 
PLOT 

23  tf  UAllJ.Nr.AP)  .OK.  ( A(  21  .NE.AL  I .OR.  ( A<  31  .NE.AOl  .OR . 1 A ( 4 I .NE  . AT  J 1 
IGO  TO  31 

RfLAGU2)  • 1 
CALL  Lfft  IN) 

24  If  ( I Al N) .NE.AF 1 .OR. I A(N«l) .NE • AA ) .OR. < A | N» 2 ) .NE • AR 1 .OR . I A I N4 3) . N£ 

1 .Af I I CO  TO  25 
IGAIN  • 1 
NffP  • I 
CO  TO  2 T 

25  If  ( I A( N) .NE .Aft) .OR. (AIN*  n .NE. Al ) .OR.I A|N»2) .NE. AS) .OR.( AIN*  31 .NE 
(.AT) ) CO  TO  26 

(RISC  • I 
NftlP  • 1 

GO  TO  27 

26  If  U A( N) .NE . AR ) .OR. I A IN*  n .NE. AA ) .OR.( A (N*2) .NE. AC) .OR • ( A( N* 3 ) .HE 

i.ak))  go  to  n 

(SCAT  * I 
NEAP  - I 


27  UO  29  I •N,9a 
K • ) • 1 

|f  (All I.EQ. SLANT ) GO  TO  29 
2ft  CONI INUE 


GO  TO  71 
29  N • A 

If  1 1 A(  N) .NE  .AT ) .OR. I A( N* 1 > .NE. AHI.OR.l A 1N*2)  .NE  «AE  ) .OR.  I A(N»3).NE 


(.All)  GO  TO  TO 
CALI  FCilAl  IN) 

C Ai L NUKRFR  (N,N2(X1,IX1 
I T INf f P. EO.  I ) Af f T«X 1 
If  ( Nftl  P.  t U.  1 I Ai5l  T*X  I 
If  INRAP. EU.  (I  ARAr«Xl 
If  I A(N2) .EO.R IGH  T I GO  TO  4 

If  (A(N2).Ne. Slant)  go  to  n 

N • N2* I 

GO  TO  2**  • 

30  If  I I At N) .Nf . API .OR.I AIN* 1 ) .NE. AHI .OR.I A IN* 2) .NE.AI  ) I GO  TO  71 
CALI  ( UUAL  IN) 

CALL  N*)N»EH  IN,N2,XI  , |X) 

If  INf » P .Eg, i ) Af f P«X  ) 

If  ( NRI  P.f  U. 1 ) AR  |P«X 1 
If  INftAP.fg.l)  AHAP*Xl 
If  IA|N21.fU.M:>NTI  GO  TO  4 
If  (AIN2) .NE .Slant ) gu  to  71 
N • N2* 1 
GO  TO  24 


31  if  I I A| J| .NF .AO) .OR. I A(2) .NE.AU).OR.( A( 3) .NE. AT ) . OR. ( A ( 4 1 .NE • AP) ) 

IGU  TO  *4 
KflAG(72)  - 1 
CALL  LEFT  |N) 

32  If  I 1 Al N) .Nf .Aft) .OR. I A|N*1 ) .NE.Al ).0R. I A<  N*2) .NE • AS  1 .OR . I A| N* 3) . NE 
l.AT))  CiJ  TO  33 

AfLAGI ) ft)  • I 

leisc  « 1 

CALL  ICUAL  IN) 

CALL  NOKRER  IN|N2«X1»|R) 

PMSI  • X) 

N • N2*  I 

call  number  IN,N2»X1,|X) 

PHSf  • XI 
N • N2*  1 

fN,N2.Xl,|X) 

N • N2* I 

CALL  NUMftER  IN.N2fX|,|X) 
tMSf  • XI 

If  (A|N2f.E0. RIGHT)  GO  TO  4 
If  IA|N2).NE. SLANT)  GO  TO  71 
N • N2*  k 
GO  TO  32 


76 


33  M A( Nl .NF .AF ) .UR. U<N»1) .Nf .AAl •OR*  t A < Na2I .NE • AR I .OR • ( A (NA 3 1 .NE 
I .AF  M 00  f 0 3A 
RFLAO(iC>)  • 1 
I&AIN  • I 
CALL  ECUAL  (Nl 

call  number  <N,N2*XUII) 

PHAI  • XI 
N • H2* I 

CALL  NUMBER  IN,N2tXUIX) 

PHAF  • XI 
N - N2»  I 

CALL  NUMBER  IN»N2«X1,1X| 

IMAI  - XI 
N ■ N2*  i 

CALL  NUMBER  <N«N2«X1,|X) 

Thaf  • XI 

IF  U(N2I.EU. RIGHT)  00  TO  A 
IF  < Ai N2I .NE .Slant ) cc  to  ti 

N • N2* I 


3A  IF  UAIN)  .NE.AN)  .OR.U(N*II.NE.AC  ).0R.(A1Na2)«NE.AA).0R«(A(N»3I.NE 
I .ARt ) GO  TO  AO 
R F L AC  I I 9 ) • I 

inear  • ? 

CALL  EUUAL  ini 
IF  (AINI.EQ.PLEFT)  GO  TO  33 
INEAR  • I 

I • I 
00  TO  36 


35  00  37  L«1*50 
I • I 

N > N»1 

36  CALL  NUMBER  IN*N2»XU1X) 

XNPI  11  •XI 

N • N?*  I 

call  number  IN*N2«XUIXI 

TNPI I I • XI 
N • I 

CALL  NUMBER  IN,N2fX|,IXI 
2NPI 11  - XL 

IF  llNEAR.EO.il  00  TO  39 
INEAR  • L»1 

IF  1 A1N2I .EQ.RIOHT  I 00  TO  38 


289 

290 

291 

292 

295 

296 

295 

296 

297 
296 
2 99 

300 

301 

302 
305 
30A 

305 

306 

307 

308 

309 
510 

311 

312 
315 
31A 

315 

316 
31  7 

318 

319 

320 

323 

32A 

325 

326 

Hi 

330 

351 

332 

333 
33A 

335 

336 


N • N2 
37  continue 


GO  TO  7 I 

38  N2  • N2*| 

INEAH  a |^4fAR-t 

39  IF  1 A1 921 .EO.R lOHT I 00  TO  A 
IF  <A1N2I.NE.SLANTI  00  TO  71 
N > N2*  I 

00  TC  32  I 

AO  IF  ||A(N|.NE.AB).OR.(AIN»l).NE.AA) .ORa I A I Ma 21 .NE • AC ) .OR • I AINA  31 .HE 
l.ARII  GO  TO  A1 
Kf  LAG! I 71  • I 
tSCAf  • I 
call  EOUAL  ini 
call  number  IN,N2*XUIX) 

PHI  I • Xl 
N a n2a  I 

CALL  NUM«IR  (N.N2tXWlXI 
PHIF  • » 

N • N2a  I 

CALL  NUMBER  iN»N2«Xl,IXI 
fHll  • XI 

N • N2»  I 

call  number  IN.N2»XI.|XI 
iHiF  - XI 

If  ia<n2i.eq.right|  ro  a 
IF  IAInH.NE.S  ' 

N • N2A  I 
GO  TO  32 


.SLANT!  00  TO  71 


A1  IF  1 1 AINI .NE.AC) aOR.IAINAll .NE.AUI.OR.I AINA2I .HE.AR  l.aR.IAlNAllaNE 

l.ARI I 00  TO  A3 
RFLAOf  I 51  • 1 
INR  • I 


??  )iiici?{'?KisMri  CO  TO  . 

N • RAl 

IF  I AIR t.EO.SLANT I 00  TO  32 
A2  CONtlNUE 

00  TO  71 


357 

358 

359 

!:! 

3A3 

3AA 

365 

366 
67 

368 

369 
350 
>51 

ill 

35A 

355 

356 

357 

358 

359 
560 
36  1 

362 

363 
566 
565 
366 
36? 
568 

369 

370 

371 

373 
3 76 
J T5 
3 76 
JTf 

378 

379 

Hi 

382 
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) IF  ( 4 At  Ni .Nl .AS)  .UR. I A<N»II .NF. AT).UR.( AIN«2).NE. AE) .OR. 1A(N»3I .NE 
4 .AP) ) on  TO  71 
CALI  EQUal  iNt 
CALI  NUMH4R  IN«N2iAl»lX) 

SUP  • *) 

IF  4 ALS? ) . L W.H lOHT  ) GO  TO  A 
IF  4A4 S2) .NL .SlANf)  GO  TO  71 
N > N7* I 
GO  TO  i? 

FFEO  POINT 

k IF  4 4 A|  1 ) .Nf .AF ) .OH. 4 A1 2) .NE.AE ) .OR*  4 A4 i) .NE.AE) .CP.lAl^) .NE .AO) I 
IGO  TO  Mb 
RFL A041 3)  • i 
GO  TO  46 

i IF  4 4 A4  I ) .NE . AG ) .OR. 4 A4  2 ) .NE. AE ) .OR. 4 Al 3 ) .NE.ANI .OR. 1AI% I .NE .AE) ) 

ion  T(j  !»<; 


AFiAG42i)  « 1 

N04N  • 0 

CALL  LIFT  4N) 

tail  NOM0{R  4N.N2,XI,IX) 

NOEN  • NOEN»| 

AOEN4NutN)  • IX 

IF  ( A4  N2»  .EQ.R  UhT  I GO  TO  4 

N • N2* I 

call  number  4N,N2*XUIXJ 
VMAG  • XI 
N • N2*  I 

CALL  NUMBER  4N«N?.X1. IX) 

VOEG  • Xl 

VREAl  * vHAg*CO$4 VOEG/RAOl 
V|MA^  • VMAG«SiNI VnEG/RAOl 
VOLTINOEN)  • CNPIX4VREAL.VIMAG) 

IF  4 AINU  . 4 O.R  IGHT  ) GO  TO  A 
|F  4 Al n}) .NT . SL ANT ) GO  TO  71 

IF  4 4A4N2I  .EO.Sl  ANT  } .and.  4 A(N2»1 1 .EG.ei  ANKM  GO  TO  AB 

N • N2* I 

GO  TO  A 7 

RE  AO  I A . 76  i A 

tCARO  • ICARO»l 

WRITE  46.771  ICARO.A 

N • I 

CALL  8LNK  lA) 

CO  TO  A 7 


OESCR IP  T 1 ON 

a9  if  ( 4 AI  1 ) .NI  .AO) .or. 4 A(2I .NE.AE I .or. 4 A4  3) .NC.AS) ,0R.1A< AI.NE.ACn 
IGO  70  b^ 

■FLA^4W»  • 1 
J • 0 

call  LEFT  4N) 

30  call  number  4N.N2.Xl.IX) 

J • JM 

NW  • j 

I A4  J)  . IX 
N • N?*  I 

CALI  NUHBFR  IN.N2(X1.1X) 

I HI  J)  • IX 

If  4A4N2).EU.R|CHT)  GO  TO  A 
IF  AIN2) .NE .SLANT)  GO  TO  71 

IF  UAIN2).fQ.UANt).ANO.U4NAM.EO. BLANK))  GO  TO  51 
N * N2* I 

GC  TO  50 


jCARO 


t'lt 


AR0»1 


call  6LNK  4Ai 

. rh  ICARO.A 


WRITE  4b 
N • I 
GO  TO  50 


52  (F  4 4 A(  LI  .NE  .AGl  .OR.  1 44  21  .HE.4E1  .OR.  4 401  .HE.  40)  .OR.  1 4 1 A)  .HE  • 4M4U 
I GO  T Cj  55 
RFLAG4I2I  • 1 


JJ  • 0 
Call 


- If F T 4N) 

53  call  number  (N.N2,X1,IX) 

JJ  • JJ*1 
NP  • JJ 
XIJJ)  • XI 
N • N2a 1 

CALL  NUN0ER  4N,N2tXl.|Xl 
VIJJ)  • XI 
N • N2* 1 

^ALL^NUMBER  IN.N2fXl,|X| 

88 18  T, 

If  44A|N2l.EO.SLANT).ANO.lA(N«ll.EQ.BCANKn  GO  TO  54 
N • N2* I 
GO  TO  51 


385 

386 

387 
386 

389 

390 

393 

39A 

395 

396 

397 
?98 
399 
AOO 
AOI 
A02 
a03 
AOA 
AOS 
A06 
A07 
a08 
409 
AlO 
Al  I 
A12 
Al  3 
AlA 

a15 

Alb 

AI7 

AlS 

a19 

A20 

A22 

A23 

A2A 

a25 

A26 

a27 

A28 

A?9 

A30 

A3l 

A32 


A33 
A 3A 
A)5 
A36 
A37 
A38 
AJ9 
AAO 
AAl 
AA2 
AA3 
AAA 
AA5 
AA6 
Aa7 
AA8 
AA9 
A50 

A5) 

A5A 

A55 

A56 

A57 

A5S 

A59 

A60 

A6l 

A62 

A63 

A6A 

A65 

A66 

A67 

A6i 

469 

A70 

472 

473 

474 

475 

476 

47B 

479 

4*0 
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Q 


RIAO  (5. 76)  A 
ICARO  • ICARO«l 
MftMfc  I6,r7)  ICARO.A 
CALL  BLNK  <A) 

N • \ 

CO  TO  5 ) 

INTFAVAL  for  CALCULATION 

IF  n AM  I .NF.AI  ) .on.  I AU)  .NE.AHI  .OR.I  Al  ))  .NC.AT  ) 
ICC  to  56 
KFLAGI^U  - I 
CALL  LfFT  IN) 
call  NUMBFR  1N«N?«X1,|X) 

I NT  • IX 

IF  IA(N2).E0.AICHTI  GO  TO  6 
GO  TO  71 


OR. IA<6) .NE .AE  I I 


t»  IF  M AM  ) .NE  .AG)  .On.l  A(2)  .NE.AR)  .on.  I A|  3)  .NE.AOl 
IGU  10  66 
AFLAGM5)  • 1 
AFLA0M6)  • 1 


iGfcO  • Z 
CALL  LEFT  (N) 

J IF  I 1 A(  N)  .NE  .API .on. ( AIN»l I .NE.AE ) .on. t A(N»2) .NE. 
l.AFM  GO  10  bB 
IGRO  • 1 
GO  TO  6R 

B IF  UAlNT.Nf  .AC).0n.lA(N«l).NE.A0).0n.(A(N»2).NE. 
1 .AO) I GO  TO  59 
E“R  • iO. 

5iGR  • .02 
GO  TO  6R 

9 IF  M A( N) .NE  .AP I .OR.I AIN» 1 1 .HE. AO) .OR. I A (N«2) .NE. 
I .AB)|  GO  TO  60 
ERR  • R. 

SIGR  • .001 
GO  TO  6R 

0 IF  ilAI  N)  .NE.ASI  .0R.UIN«1)  .HE  • AE  I .OR.  ( A ( Nr  2 ) .NE  « 
ERR  • 80. 

SIGR  - R. 

GO  TO  6R 

1 IF  M Al N) .NE .AH) .OR. I AlNRl ) .NE. AE ) .OR.I A<N«2) .NE . 
1 .AG) ) GO  TO  62 

CALL  equal  IN) 

CALL  NuneER  IN,N2,X1«IX) 


OR.(AU).NE.AU)  ) 


AR).0R.|A<N»3).NE 


AO) .0R.(AINr3I.NE 


A0I.0R.IAINr3).NE 


AA) ) GO  TO  61 


A1).0R.(ACNr3).NE 


hgt  • XI 

IF  I A(N2l .ro.p IGHT ) GO  TO  R 
IF  I A< N2) .NE . SLANT ) GO  TO  71 
N « N2* I 
GO  TO  57 

2 IF  I I A I N) .NL  .AC) .OR.! AINrL ) «NE • AO) .OR. I A I Nr  2) .NE . AN)  .OR  • I A (Nr  3)  .NE 
1 .AO) ) GO  TO  63 

CALI  equal  IN) 

CALL  NONttER  IN,N2«X1,)X) 

SIGR  • XI 

IF  I AIN?) .E Q.R IGHT I GO  TO  R 

IF  IAIN?). Nt. SLANT)  GO  TO  71  • 

N • N2*  I 
GO  TO  57 

3 IF  (lAINl.Nf  .A3).0R.(AINR|).NE.A1).0R.I  A|Nr2).NE.AEI.0R.IAINr3I.NE 
l.AM)  GO  TO  71 

CALL  E CUAL  IN) 

CALL  NJN8ER  4' 


CALL  NJNB 
ERR  • XI 


IN.N2fXl.tXl 


IF  I A (N>) .EU. RIGHT  I CO  TO 
IF  (AIN?) .NE .SLANT)  GO  TO 


IF  IAIN? 
N > N2r 1 
GO  TO  5t 


00  65  K»N.80 

IF  IAU).N. RIGHT)  GO  TO  R 

N • ARl 

IF  (AIK ).E0. SLANT ) 60  TO  57 
CONTINUE 


b IF  II  AMI  .NE  .AS)  .OR.I  A(2I  .NE.ATI  .0R.(  A(  3).NE.A0)  .OR.(AIR).NE.AP) ) 
IGO  TO  67 
IFLAG  - 2 
RETURN 

7 IF  ilAI  |).NE.AC).0R.IA(2).NE.AH).0R.(A|  3 1 .NE . AA ) . OR. U I R ) .NE • AN) ) 
IGO  TO  68 

TflAG  • 3 
RE  TURN 

8 IF  II Al 1 I.NE.AE I .OR.I ACai.NE.ANI .OR. lAI 3).NE.A01 I GO  TO  71 


IFLAG  • 1 
RE  TURN 
IFIAG  • 5 
RE  TURN 
IFLAG  * R 
RETURN 
HSG  • I 

RFLAGMO)  • ICARO 

UFLAG.NE.5)  write  l6f7B) 
IFLAG  • 5 
RETURN 

IFLAG  • 6 
ICARO  • ICARO-I 
RETURN 


l//^/5xI loATA  CARDS*//! 


ELiRNAT  IIOAll 

FORMAT  |6X  f U f 2X  t iOAl  ) 

format  I • ilBlB  END 


CARO/STOR  CARO  NISSING* 

79 


RITE 

PURPOSE:  to  generate  a list  of  branch  currents  from 

the  input  loop  currents. 

METHOD:  The  generation  of  branch  currents  is 

accomplished  in  the  DO  LOOP  ending  at  statement  2.  The 
branch  currents  are  stored  in  CJ  (I)  by  the  latter  part  of 
the  DO  LOOP  ending  at  statement  3.  If  the  branch  currents 
are  requested  for  output  (IHR  positive) , the  DO  LOOP  ending 
at  statement  5 accomplishes  this. 

CALLED  BY:  GANT1 

GFFLD 

CALLS  TO:  NONE 


80 


•k 


c 

c 


I 


c 

c 


c 

c 


c 


SUBAOUT iNe  RITE  < I A 1 1 B . I NN,  1W«. 1 1 • I 2 f I ) t HO • NOtNH , C J i CG t IGROI 
CCHPifi  c JiDrCGli  ),CJA,CJB 

OIMFNSION  1411»»  IBlll.  inn*  IZIll,  Dili*  M01INM*A1*  NOIlt 
AMAX  • .0 


00  ) K* 1 «NH 
RA  • lA (R I 
RB  • IBIR) 

C JA  • I .0, .01 
CJB  • I .0*.0I 
HOR  • NOIRI 


00 


2 M •!  ,NOR 
1 • HOI  R,  In 
M • 1. 

Jf  IRB. 10. m I ))  CO  TO  1 

n IRB.  10.  in  I M 1 . 

CJA  • CJA*F|«CJII I 
CO  TO  2 

IF  iRA.fO.m  ID  F(»*U 
CJB  • CjB»F|*CJi I I 

continue 


CGIRI  • CJA 
KR  a R*RN 

CCIRRI  • CJB 
ACJ  - CABSICJA) 

ecj  * cabsicjbi 

IF  (ACJ.OT.ANAXI 
IF  (BCJ.GT.AMAX) 
CONTINUE 


AMAX-A 

AHAX-B 


J 

J 


IF  I IWA.CT.OI  GO  TO  A 
return 

A IF  (ANAX.LE.O.I  amax-u 
MRITE  16*81 
NNC  • NM 

IF  IIGRO.GT.O)  NMG  • NH/2 


00  5 R* 1 1 NHG 
CJA  a CgU) 

RR  a R*NH 

CJB  • CGIKRI 
CCJA  • CABSICJA) 


1 

2 
3 
A 

5 

6 
1 
8 
9 

!i 

13 

lA 

15 

IS 

19 

20 

i\ 

\l 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 

t\ 

A3 

AA 

A5 

A6 

AT 

A8 


i 


CCJ6  • CABSICJBI 
ACJ  • CCJA/AMAX 
BCJ  • CCJB/ANAX 
PA  « . 0 
PB  • .0 

IF  lACJ.GT.O.)  PA 


IF 


I8CJ._  . . . 

WRITE  <6*7)  R*lAIRlfCJA,CCJA.ACJfPA 


• 57.29578*Ar AN2(A1HA0ICJA) .REAL  I CJA 1) 


.CT.O.I  P8  • 57.29578RArAN2lA 


AlMAGiCJe. . . . . 

,U|RI*CJ6*CCJB*8CJiPB 


/eitPEAllCJB)  I 


WRITE 

return 


I 6*61 


6 FOR  NAT  MHOl 

7 FORMAT  l2X*t2*2l2K.12*2X.Ell.5»lX*€n*5*U,Ell.5*lX*Elt.5.U*F6.l) 
1) 

• FORMAT  U2IA6X  ••NCRNAtUEO*  .BX  I/*  SEG*«2I'  NODE  * ■ AX  * 'RE  AL  • *6X  * • 1 MA 
IGINARV*  •3X,*HAGNt  TuOE*  * 3X *• HAGN  ITUOE* , 3X **  PHASE  • ) I 
END 


A9 

50 

51 

52 

53 
5A 

55 

56 

57 
56 

59 

60 
61 
62 
63 
6A 

65 

66 

67 

68 
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SART 


PURPOSE:  to  soct  data  for  polar  plot. 

METHOD:  This  subroutine  sorts  the  values  of  the 

points  to  be  plotted  by  the  polar  plot  package  starting  with 
the  greatest  positive  value  of  y to  the  greatest  negative 
value.  In  the  DO  LOOP  ending  at  statement  1,  the  value  of 
(x  ,y  ) is  interchanged  with  the  value  of  (x  ,y  ) if  y is 

i i j j j 

greater  than  y . 

i 

CALLED  BY:  POLPRT 

CALLS  TO:  NONE 


c 

c 

c 


c 


c 


$Ut<fiOUMN€  SASr  lOATAX.OJkTAY.Nl 
OlHFNSinN  6aTmY(»00) 

This  RuuriNf  sorts  data  ih  oatay  sy  magnituoc 


HN  • I 

00  2 i • > »NN 

NM  - I ♦ I 


00 
tf  . 

ST[>A  - OAfAYl  I ) 

OAl A V(  n • OATAYI J) 
DATA  Y1 J)  • STOA 
SrOR  • OATAXCll 
5aTA  XI  I)  • OAMRU) 
DATA  X<  m srOR 
1 CONTINUC 

I CONI INUE 


1 J>NH  N 

IOATAV  M .oe^OATAYl J)|  GO  TO  1 
Uyi* 


Rf TORN 
(NO 


I 

i 

A 

s 

6 

T 

S 

9 
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SCANT 


PURPOSE:  to  calculate  the  Mutual  iapedance  between 

filamentary  monopoles. 

2 
C 


1 

2 


SOURCE 

DIPOLE 


D 

1 


RECEIVING 

DIPOLE 


Z = P11  ♦ P12  ♦ P21  + P22 


METHOD:  In  the  induced  emf  formulation,  the  mutual 

impedance  of  coupled  dipoles  is 


Z = - 


/ 


I (t)  E (t)  dt 

2 I 


where  denotes  the  current  distribution  (normalized  to 

unit  terminal  current)  on  dipole  2,  and  E^(t)  is  the  field 

of  dipole  1 when  it  transmits  with  unit  terminal  current. 
Distance  along  the  axis  of  dipole  2 is  denoted  by  the 
coordinate  t.  E^  may  be  expressed  as  the  sum  of  the  fields 

from  each  of  the  monopoles  comprising  dipole  1.  Furthermore, 
the  integral  is  the  sum  of  the  integrations  over  each  of  the 
monopoles  comprising  dipole  2.  Thus,  the  dipole-dipole 
mutual  impedance  may  be  expressed  as  the  sum  of  four 
monopole-monopole  impedances. 


It  may  be  convenient  to  draw  the  above  figure  in  terms 
of  monopoles  with  the  current  distribution  shown  as  dotted 
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lines 


(The  Bonopole  letters 


remain  the  same.) 


A 

' P11  / 

c 

A 

\ PI  2 > 

; \ 

» \ 

J 

B 

; P21  / 

• 1 

’ 1 

c 

B 

! P22  ' 

t 1 

/ \ 

2. 


The  surface  impedance  is  calculated  just  above  statement 

B0 1 denotes  J / J where  J and  J are  the  Bessel 
0 1 0 1 


functions  of  order  zero  and  one  with  complex  argument,  ZAPG. 
It  is  assumed  that  all  the  wire  segments  have  the  same 
radius,  conductivity  and  surface  impedance. 


In  the  DO  LOOP 
the  segment  length 
of  the  shortest  an 
the  segment  length 
thin-wire  theory, 
by  RETURN  to  the  m 


ending  with  statement  3,  SCANT  calculates 
s 0 ( J)  . DHIN  and  DMAX  denote  the  lengths 
d longest  segments.  If  the  wire  radius  or 
s are  clearly  beyond  the  range  of 
N is  set  to  zero  at  statement  4 followed 
ain  program  to  abort  the  calculation. 


At  statement  5,  the  program  selects  a segment  K,  and  a 
few  statements  below  this  it  selects  another  segment  L.  K is 
a segment  of  test  diple  I,  and  L is  a segment  of  expansion 
mode  J.  The  mutual  impedance  between  segments  K and  L is 
obtained  by  calling  subroutine  CCS  or  GGMH.  In  statement  18, 
this  impedance  is  lumped  into  C(NNH).  The  mutual  impedance 


Z between  dipoles  I and  J is  the  sum  of  four 
segment-segment  impedances. 


The  variables  IFLAG  and  JFLAG  are  used  if  a ground 
plane  is  present  for  the  calculation  of  the  mutual  impedance 
elements.  If  IFLAG  is  equal  to  JFLAG,  the  mutual  impedance 
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teras  will  not  have  the  effects  of  a ground  plane  since  both 
monopoles  lie  on  the  same  side  of  the  ground  interface.  If 
the  monopoles  are  on  the  opposite  sides  of  the  interface 
(IFLAG  not  equal  to  JFLAG) , the  reflection  coefficient 
correction  must  be  applied  to  the  mutual  impedance  elements. 
This  same  technique  is  applied  in  subroutines  GNFLD  and 
GFFLD. 

In  SGANT,  segment  K has  endpoints  KA  and  KB,  and  segment 
L has  endpoints  LA  and  LB.  It  is  convenient  to  think  of  KA 
and  KB  as  points  1 and  2 on  segment  K,  and  LA  and  LB  as 
points  1 and  2 on  L.  The  four  segment-segment  impedances 
can  be  defined  as  P(IS,JS).  The  first  subscript  IS  refers 
to  the  terminal  point  on  segment  K,  and  the  second  subscript 
JS  refers  to  the  terminal  point  on  L.  Thus  IS=1  or  2 if 
dipole  I has  its  terminal  point  12(1)  at  KA  (point  1)  or  KB 
(point  2) , respectively.  Similarly,  JS=1  or  2 if  mode  J has 
its  terminal  point  12  (J)  at  LA  or  LB.  The  impedances 
P(IS,JS)  are  defined  with  the  following  reference  directions 
for  current  flow:  from  point  1 toward  point  2 on  each 
segment.  If  dipole  I has  this  same  reference  direction  on 
segment  K,  PI=1;  otherwise  FI=-1.  Similarly  FJ=1  or  -1  in 
accordance  with  the  reference  direction  for  mode  J on 
segment  L.  In  statement  18,  P(IS,JS)  is  multiplied  by  FI  and 
FJ  before  its  contribution  is  added  to  2 

ij 

Subroutine  GGHB  calculates  the  impedances  Q(KK,LL)  which 
are  like  the  P(IS,JS)  but  have  different  conventions  for 
reference  directions  and  subscript  meaning.  The 
transformation  from  the  Q impedances  to  the  P impedances  is 
accomplished  in  the  DO  LOOP  ending  with  statement  13. 

If  the  wire  has  finite  conductivity,  the  appropriate 
modification  is  applied  to  the  impedance  matrix  just  above 
statement  15.  The  terms  arising  from  the  dielectric  shell 
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oil  an  insulated  segment  are  obtained  from  subroutine  DSHELL 
just  above  statement  16.  Finally,  the  lumped  loads,  ZLD,  are 
added  to  the  diagonal  elements  of  the  impedance  matrix  in 
the  DO  LOOP  ending  at  statement  23. 

K is  a segment  of  test  dipole  I,  and  L is  a segment  of 
expansion  mode  J.  When  the  segment  numbers  K and  L are 
equal,  SCANT  calls  GGMM  to  obtain  the  mutual  impedance 
between  two  filamentary  electric  monopoles.  These  monopoles 
are  parallel  and  have  the  same  length.  Honopole  K is 
positioned  on  the  axis  of  the  wire  segment,  and  monopole  L 
is  on  the  surface  of  the  same  wire  segment.  Thus,  the 
displacement  is  equal  to  the  wire  radius.  The  two  monopoles 
are  side-by-side  with  no  stagger. 

when  segments  K and  L intersect,  SCANT  again  calls  GGHH 
for  the  mutual  impedance  between  the  two  filamentary 
monopoles.  Monopole  K is  situated  on  the  axis  of  wire 
segment  K,  and  monopole  L is  on  the  surface  of  wire  segment 
L.  The  axes  of  segments  K and  L define  a plane  P,and 
monopole  K lies  in  this  plane.  Monopole  L is  parallel  with 
plane  P and  is  displaced  from  it  by  a distance  equal  to  the 
wire  radius. 

CALLED  BY:  MAIN 

CALLS  TO:  CBES 

DSHELL 

GGMM 

GGS 
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gtsi  COPY 


SUBHOUT INE  SGA^T  U A . t B • I NH • i4t . 1 SC • 1 1 • I 2 • I 3« JA t JS t MOt N . NO .NM.NP • A 
IM.BS.C  tCG0.CM*4.0tEP2t(^3,ET  A|FHKGAH,s60.XtV«2,k0f  2S*(M«  1ga6| 
CCHPLt*  tfiB 

COMPc(Ji  iGfZH.^S.EGDiGOtCGOStSGOStSGpr  .BQI 

CUMPiE«  PI  1 .pI?.P2  1.^22,01 1 .M2«0h«02itCP2«EP.ETA*GAM,EP) 

CCHPiEx  f PStLAfdMEA.BETA.iAAG 

COMPi  E X P(2, 2i .Ql 2fi) fC06l 1 I.SGOt l),C( 1) .ZlOl II 
01**ENSI0N  xUli  ^Ulf  Ml).  6tl)(  lAU).  Idlllf  MQIINM.AI 
OIMCNSION  iiul.  iMl).  1)111,  JAM).  J&uT,  NOtllf  isdlll 
data  EO,TP.UO/B.0SA€>li.6«2B)leil.2566E-A/ 


tP  • £P3 

ICC  • IN*N*N)/2 


00  1 l«  1 , ICC 
C(l)  • 1 .0..0) 


2S  • I . 0. .01 
If  ICHN.ie.O.)  GO  10  2 
ONEGA  • fP*fH2 

CPSIIA  - CMPiUE0,-CNN*l.e6/0N£GA) 

CWEA  • « .0. I • I •QN^GA*EPS1 LA 

BEtA  • OMEGA*SORT(UO)ACSOAr(CPSlLA-EP) 

2APG  • Of  rA*AN 
CALL  CBES  (ZARGiBOl ) 

2s  • erTA*eoi/CriEA 
2H  • ;s/i rP*AM*GANI 
OnIn  • I .t  )0 
ONAX  • .0 

00  ) J* I .NN 
K - lAljI 
L • iB«  J) 

0(J)  • SQ«T((X<x|-X(t))**2A|V<lC|-yUI)PP2#l2MC)-2ILMPP2) 
If  I0< Jl.LT.OMlN)  ONIN-O(J) 

If  lOlJI.GT.ONAX)  0HAX*0(J) 

EGO  • CEXPIGAN*0( J)) 

CGOIJ)  - (E&0»l./EG0l/2. 

SGOlJl  • UG0-t./E00l/2. 

If  (0M|N.Lr.2.*AN|  GO  TO  A 

If  ICABSl GAMAAN) .GT.0.061  GO  TO  A 

If  ICABSIGANaONAXI.GT.).)  go  to  a 

If  IAM.GT.O.I  go  to  5 

CONTINUE 

N*0 

NRirC  <6.2A)  AN,ONAXfOMlN 
WAITE  C6.2SI 


S OU  19  K >1 «NW 
IfLAG  ■ 0 

If  n IGRO.GT .0) . ANO. IK.GT.NN/2 ) ) IfLAG*! 
NOR  • NOtKl 
KA  • UU) 

RB  • IBIK ) 

OR  « 01 K ) 

COOS  • CCOIR) 

SGOS  • SGOU) 

00  19  I *1  .NN 
JfLAG  • 0 

|f  11 IGRO.GT.O) .ANO.IL.GT.NM/2) I JfLAG*! 
NUl  • NOIL) 

LA  • lAU  I 
LB  > IBIL  1 
OL  • OU ) 

SGOr  • SGOILI 
NIL  • 0 

00  19  11*1  I NOR 

1 • MOl K, II > 

NN  « 11-1 lAN-l 191-11/2 
fl  • 1. 

If  IRB.  EQ.  121  in  00  TO  6 
If  IRB.EQ.IIIIM  fl*-l. 

1$  • 1 
GO  TO  1 

h If  IKA.EO.Dl  II)  f 1*-1. 

IS  • 2 

7 00  |9  jj«l ,NOL 
J • NOU  f JJI 

MNM  • NNAJ 

If  I I .GI.JI  GO  TO  19 
f J • 1. 

If  IlB. FQ. 121 J) I GO  TO  B 
If  ILB.EQ.Itl J) I f J*-l« 

JS  • I 
GO  TO  9 

B If  (LA.CQ.IIIJ)  I f J— U 


B If  (LA.CQ.IIIJ)  I f J— U 

9 ff  INic.NC.OI  GO  TO  IB 
nil  • 1 

If  IK.IO.L)  GO  TO  lA 

INO  * UA-XA|9Ui-XA|9(LA-XB|9lLi-RB) 
NGAp  • JGAO 

If  UfLAGaO.JfLAGI  IGA0--1 
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BtSUVAriS'i  COPY 


( INp.eo.O)  GO  to  10 
>rGMfNr$  K AND  L SHAIli  NO  ^OtMTS 


im-AWiili 

igap  • m:;I6 
GO  ro  16 

SEGMENTS  K AND  L SHARE  ONE  ROIHT  (THEY  INTERSECTI 
10  KC  • 0 
JM  • KB 
JC  • KA 
KE  • I 

INO  « (KB>LAI*<KB-IB| 

Ie  1 IND.NE.OI  GO  TO  ll 
JC  • K8 
KE  • -I 
• KA 
KG  « 3 
tl  lO  • 3 
JR  • LA 
Lf  • -I 

IF  UB.EQ.JCI  GO  TO  U 
JP  • LB 
LF  - I 

12  k^N*-^Kf*i.F 
CESl  • 


UM  JCP*1Y(  JRl-VCXI  l•UUN|-y^JCl'•♦U 

SGOTtCPSl.ETRtGAMtOUf  U iQllfa 


CESl  • liXUP)>XUCM*(Xl  JMI-Xl  jcn» 
l<  JP»-2  UClMlh  JN^ZI  JCMl/JOKApL} 
CALL  GGMm  (•OtOK,«OfOL«AN«CGOSi S46S* 
1 ) .0(2*1  ) .OU,2)  I 


00  13  KK-1,2 
AP  • UBSUK-KC) 


DO  13  Lfl.2 
LP  • UaSlLL-lGl 
PUP.LP)  • SGNAQ(KK«LLI 
13  CONTINUE 


: GRO*NGRO 


GO  TO  10 

K»L  (SELF  KiACtlON  OF  SEGMENT  Kl 

lA  on  • (.O..Oi 

012  • } .0..0{ 

IE  (CMM.L|.0.)  GO  TO  19 
GO  • ^NApK 

1^1*-  ?G*4 sS^S*^^0S-CDI/2. 
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90 

99 

00 

01 

_03 

lOA 

OS 

s? 

00 

09 

0 


n 

I! 

|0 

31 

II 

it 

39 

37 

1$ 

AO 

ii 


012  • £GP(uO*CGOS'SGOSl/2. 

Pl2  • 1 .0..C1 

IF  II SCK.EO.OI  GO  TO  16 

|F  IBM. IE. AN)  GO  TO  16 

call  OShELL  lAN*BN.OK*CGOStSOO$«EP2*EPtETA,GAM.Pll*P12i 

“ §{i : ?U:g!i 

CALL  GGNN  I.OiDK* .0«0K,AM,CG0StSC0$*SG0S*l.*ETA,GAM»Pll«P12iP2i«P2 
121 


on  • p) noi  1 

GW  • P1S012 

?}  :M  : 811 


II  P 

P 


„ iii : ip 

Pt2.2»  • oil 
IE  (kA.NE.LA)  go  TO  17 
GO  TO  ■ 


. - 

10  ClHMM)  • c 
19  CONTINUE 


MMN|aE|aEJ*PI rs*2S) 


23  1-  ,N 
• 41-  J*N- 

■ NK* 


DU  23 

MN 

1 J 

JJA  • JAII  I 

• J4^ 


I M'j-n/2 


II 


JJ6 


lll2.|^jlBUIII  Jl-JUNM 


JJJ  • JJA 


PQ  22  K*n2 
NOJ  • NOIJJJI 


00  21  iJ*l*NOJ 

u ,0  n 


V’ 

(, 
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BESl  AVArJSu  COPY 


C 

c 

c 

c 


SORT 


PURPOSE:  to  define  the  set  of  dipole  modes. 

METHOD:  In  the  DO  LOOP  ending  at  statement  3,  the  set 

of  dipoles  is  defined  by  filling  the  vectors  II  (I)  and  13(1) 
(the  endpoints  of  dipole  I)  ; 12  (I)  ( the  terminal  point  of 

dipole  I);  and  the  vectors  JA(I)  and  JB(I)  (the  monopoles 
comprising  dipole  I)  with  the  node  numbers  and  segment 
numbers,  respectively.  The  DO  LOOP  ending  at  statement  8 
determines  MD(J,K)  (the  list  of  dipoles  sharing  segment  J) 
and  ND  (K)  (the  number  of  dipoles  sharing  segment  J)  . 

CALLED  BY:  MAIN 

CALLS  TO:  NONE 


c 


c 


c 


c 


c 

c 

c 

c 

c 

c 


SUBROUMN(  sour  ( |A,lB»ll»l2tl3i  JAt  JBfNOtNO*NM,HP«N»MAJI,HINtlCJ«lN 
INI 

OINFNSION  JSPI20) 

12111.  nui.  JAIU.  JBllI 

NOUI.  NOCINM.AI 


IBllI 


OiNFNStON  mil 

OJMFNSlON  lAlll 
I • 0 

DO  > l.NP 
NJK  - 0 


00  I J-liNH 

INO  - 1 |A1JI-R)*1 IBIJI-KI 
If  I INO.NE.OI  GO  TO  1 
NJR  • NJK»1 
JSPINJKI  • J 
CONT/NUe 

NOO  • NJR-l 
IF  1N00.LE«0I  GO  TO  3 

00  2 I NO* 1 .NOO 

1 • I»1 

IF  tl.GT.ICJI  CO  TO  2 
I PO  • 1N0»1 
JAI  • JSPI IMOI 
JAIII  • JAl 
JBI  • JSPIIPOI 
JBlll  • JBI 
Mill  • Ta(JAP 
If  llAlJAr 


• EO.RI 
K| 


121  I » • R 
13111  • lAlJBll 
IF  ilA(JBIl«ld. 
CCNTINOE 

CONT INUE 

H • I 

00  A J«liNM 
NOUl  ■ 6 

00  A R*l»A 
NOIJ.RI  • 0 


111  D-IBl  JAI) 
I31I)«IBI JBI) 


GT.ICJ)  111 


ICJ 


1 

2 

3 

A 

5 

6 
7 
B 
9 

10 

U 

12 

13 

lA 

15 

16 
IT 
IB 

19 

20 
21 

I! 

26 

\l 

29 

30 

l\ 

33 

3A 

35 

\^r 

30 

39 

AO 

A1 

A2 

43 

AA 

A5 

A6 

A7 

40 


C 


c 

c 

c 


c 


00  B I- It 
J • JAI  11 


III 


00  7 L-lt2 
NOIJI  • N0(JI»l 
R • 1 

N ■ 0 

NJR  • NOIJ.RI 
IF  INJR.NE.OI  GO  TO 
N • 1 

NOIJ.RI  • I 

R • R*1 
IF  (R.GT.A 
IF  IN.EO.O 
’ J - JBI  I I 

I CONTINUE 

NlN  • 100 

NAX  • 0 


00 


CO  TO  2 
GO  TO  5 


0 9 J»l, 
QJ  * NOI 


7riN0j“&f Inaxi  NAX-NOJ 

If  inoj.lt.nin)  nin»ndj 


AETUAN 

ENO 


A9 

50 

li 

5A 

55 

56 

57 
50 

59 

60 
61 
62 

63 

64 

65 

66 
AT 


1 
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SQfiOT 


PURPOSE;  to  solve  the  set  of  simultaneous  equations  to 

determine  the  currents  on  the  thin  wire  structure. 

METHOD;  This  subroutine  considers  the  matrix  equation 

ZI  = V which  represents  a system  of  simultaneous  linear 
equations.  NEQ  denotes  the  number  of  simultaneous  equations 
and  the  size  of  the  matrix  Z. 

On  entry  to  SQROT,  S is  the  excitation  column  V.  On 
exit,  the  solution  I is  stored  in  S.  Z(I,i7)  denotes  the 
symmetric  square  matrix.  Also  on  entry,  the  upper-right 
triangular  position  of  Z(I,J)  is  stored  by  rows  in  C(K)  with 

K = (I  - 1)  ♦ NEQ  - (I  * I)  / 2 ♦ J - 

If  112  = 1,  SQROT  will  transform  the  symmetric  matrix 
into  the  auxiliary  matrix  (implicit  inverse) , store  the 
result  in  C (K)  and  use  the  auxiliary  matrix  to  solve  the 
simultaneous  equations.  If  112  = 2,  this  indicates  that  C(K) 
already  contains  the  auxiliary  matrix. 

The  transf romation  from  the  symmetric  matrix  to  the 
auxiliary  matrix  is  accomplished  in  the  DO  LOOP  ending  at 
statement  5.  The  solution  of  the  simultaneous  equations  is 
accomplished  in  the  remainder  of  the  program. 

CALLED  BY:  GFFLD 

CALLS  TO:  NOKE 
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1 


SUBROUTINE  SQAOT  1 C ■ S 1 1 MA » I 12 1 NE Ql 

1 

COMPiEi  Cl  1 1 tSI 1 1 tSi 

2 

N • NC  0 

• 3 

IF  1 IW.E0.2i  CO  TO  6 

4 

Cl  1 1 • CSOATICI 1 1 1 

5 

c 

6 

1 

DO  1 A-2iN 

C in)  • Cin»/Cl 1 I 

7 

8 

1 

10 

00  5 l•2tN 

11 

|H0  • 1-1 

U 

IPO*  1 » J 

13 

10  • 1 1 -1  MN-I  lAl-l  1/2 

14 

11  • lO^l 

15 

c 

00  2 l-l||NO 

\l 

tl  • U-l  1/2*1 

18 

2 

Cl  1 1 ) • Cl  1 1 l-CIL  1 )•CILI  1 

19 

c 

Cl  11)  • CSOATICI 1 1 )) 

^0 

21 

IF  1 IPO.GT .N)  GO  TO  5 

22 

c 

23 

00  4 J- |PO*N 

24 

IJ  • ID*J 

25 

c 

26 

00  3 |. IMO 

27 

MO  • IM-i t•N•|M•M-MI/2 

28 

Ml  • M0*1 

29 

MJ  ■ M0*J 

30 

3 

Cl  1 J)  • Cl  IJI-CIMJ  lACIMl ) 

31 

li 

c 

4 

Cl  IJ)  - Cl  IJ)/C< 1 11 

c 

34 

5 

CONTINUE 

35 

c 

36 

6 

SUI  - S(U/CU) 

37 

c 

38 

00  8 l•2tN 

39 

IMO  • 1-1 

40 

c 

41 

00  7 L»1 . IMO 

42 

LI  - Tl-I|aN-IL*L-L)/2*I 

43 

7 

Sil)  • SI1)-CILI)AS1L) 

44 

c 

45 

I J • ll-n*N-IMI-l)/2*l 

46 

8 

SID  - SI  tl/CI  1 1 1 

47 

c 

48 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NN  • ( N»  U*N)/2 

SINI  • StNI/ClNNI 
HHU  • H-i 


00  10  I •) *Nno 
H • N-J 
RPO  • K *1 

no  • <K-l 

00  9 l-KPOtN 

ni  • no*L 

9 S(nl  • S(ni>C(KLl*S(Ll 


nn  • no*n 

10  SIKI  • S(K)/CIRKI 

IF  ( IMR.IE.OI  CO  TO  13 
CNOP  • .0 


11 


SA  • CABsUt  ill 
IF  (SA.GT.CNOA)  CNOA-SA 

IF  (CNOA.L€.0«>  CNOA-U 

^0  I ri fN 


- SI  II 
• CABS 


A6SI SSI 


SNOA  • SA/CNOA 
PH  • .1 
IF  l‘ 

WAIT 


SA.GT.O.I  PH  • 97,299TBAATAN2IA|MAGISS)fAeALfSSM 


12  WAltE  I*. 141  1 .SNOAtSA.PHtSS 


WAITE  lArlBI 
13  AETGAN 


14  FOAHAT  UXilt3»lF10.3tlF19.7tlF10«Ot2Fl9.6l 
13  FOPHAf  IlH^I 
fHO 


49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

73 

74 

75 

76 
17 

78 

79 

80 
81 
82 

83 

84 

85 

86 
• 7 


V 


• BLNK 


M 

A 


, . EQUAL 

^EAD  • 1 1 : . . . 

• . , • • ‘LEFT 

* * NUMBER 


SORT 


.SGANT*  I } : * 


, ‘CBES 
. DSHELL 


GANT1 


GNFLD« 


GFFLD.  • : : ! 


. . GFF 


....  LINECK 

. , *POLPLOT«  • : . , 

• * • NUMB 


POLPRT* ; 


* •SART 


• • • •EXPO 
^ * • ‘GGM 


0. 1 ; , • 


, • GDISS 
• .RITE 
* • SQROT 


CALLING  SEQUENCE  OP  THE  SUBROUTINES 
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SYMBOL  DICTIONARY 


A 

ABAP 

ABAT 

ABIP 

ABIT 

ACSP 

ACST 

AFFP 

AKFT 

AN 

BN 

C 

CG 

CGD 

CJ 

CNH 

D 

ECSP 

ECST 

BFF 

EP 

EPP 

EPPS 

EPTS 

BP2 

EP3 

EP4 

ERR 

ER2 


characters  of  the  input  data  cards 

backscatter ing  phi  plane  angle  for  plotting 

backscatter ing  theta  plane  angle  for  plotting 

bistatic  scattering  phi  plane  angle  for 
plotting 

bistatic  scattering  theta  plane  angle  for 
plotting 

absorption  cross  section  for  phi  polarization 

absorption  cross  section  for  theta 
polarization 

far-zone  phi  plane  angle  for  plotting 

far-zone  theta  plane  angle  for  plotting 

radius  of  the  thin  wire  of  the  structure 

outer  radius  of  the  dielectric  shell  of  the 
insulation  of  the  wire 

elements  of  the  open-circuit  impedance  matrix 

branch  currents  for  the  structure 

cosh  /d  for  a given  segment 

loop  currents  for  the  structure 

conductivity  of  the  wire 

length  o a given  segment 

extinction  cross  section  for  phi  polarization 

extinction  cross  section  for  theta 
polarization 

radiation  efficiency 

loop  currents  induced  by  a phi  polarized  wave 

phi-polarized  far-zone  feild  of  the  dipole 
mode 

scattered  electric  field  in  the  phi  direct i 
due  to  a phi  polarized  wave 

scattered  electric  field  in  the  thet# 
direction  due  to  a phi  polarized  ws v 

complex  permittivity  of  insulatiom 

complex  permittivity  of  aabi^’t  r 

complex  permittivity  ot  jt 

EP4/EP3 

relative  di elect  tic 


ER3  relative  dielectric  constant  of  the  aibient 

■ediui 

£R4  relative  dielectric  constant  of  the  ground 

ET  loop  current  induced  by  a theta  polarized 

wave 

ETA  intrinic  iapedance  of  aabient  aediua 

ETPS  scattered  electric  field  in  the  phi  direction 

due  to  a theta  polarized  wave 

ETT  theta  polarized  far-zone  field  of  the  dipole 

■ode 

ETTS  scattered  electric  field  on  the  theta 

direction  due  to  theta  polarized  wave 

EX  near-zone  electric  field  in  x direction 

EY  near-zone  electric  field  in  the  y direction 

EZ  near-zone  electric  field  in  z direction 

EO  8.854E-12 

FHZ  frequency  in  hertz 

FHC  frequency  in  aegahertz 

GAR  intrinic  progration  constant  of  the  ambient 

■ediua 

GG  tine-average  power  input 

GPP  power  gain  associated  with  the  phi  polarized 

component 

GTT  power  gain  associated  with  the  theta 

polarized  component 

HGT  height  of  the  structure  above  ground  plane 

lA  first  node  of  a given  segment 

IB  second  node  of  a given  segment 

IBISC  indicator  for  bistatic  scatter  calculations 

ICARO  indicator  for  the  data  cards 

ICJ  diaension  corresponding  to  the  number  of 

simultaneous  linear  equations 

IFLAG  indicator  for  program  termination 

IGAIN  indicator  for  antenna  gain  calculations 

IGRD  indicator  for  presence  of  the  ground  plane 

IRC  indicator  for  the  type  of  far-zone 

calculations 

INEAR  indicator  for  near-zone  calculations 

INN  dimension  corresponding  to  the  number  of 

■ on  poles 

INT  number  of  integration  steps 
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ISC 

ISCAT 

IHR 

11 
112 

12 
13 
JA 
JB 

KFLAG 

KGBN 

LOAD 

LZD 

HAX 

HD 

BIN 

HSG 

N 

NO 

NGEN 

MN 

NPL 

OHEGA 

PH 

SCSP 

SCSI 

SGD 

SIG2 

SIG3 

SIG4 

SPPH 

SPTH 

STPH 

STTH 


indicator  for  the  insualtion 
indicator  for  bachscatter  calculations 
indicator  for  current  distribution  output 
endpoint  node  of  a given  dipole 
indicator  for  auxiliary  natriz 

terainal  node  nuaber  of  a given  dipole 

endpoint  node  number  of  a given  dipole 

first  segment  nuaber  of  a given  dipoile 

second  segment  number  of  a given  dipole 
print  indicator 

list  of  generator/feed  locations 

indicator  for  structure  load 

list  of  iapedance/load  locations 

maximum  of  the  nuaber  of  segments  connected 
to  any  one  given  node 

list  of  dipoles  sharing  a given  segment 

mininum  of  the  number  of  segments  that 
connected  to  any  one  given  node 

indicator  for  error  printout 

number  of  simultaneous  linear  equations 

total  number  of  dipoles  sharing  a given 
segment 

indicator  for  antenna  calculations 
number  of  segments 
indicator  for  polar  plot 
angular  frequency 

phi  angle  for  far-zone  calculations 

scattering  cross  section  for  phi  polarization 

scattering  cross  section  for  theta 
polarization 

sinh  yi  of  a given  segment 
conductivity  of  insulation 
conductivity  of  the  ambient  medium 
conductivity  of  ground 

echo  area  phi  incident- phi  scattered  wave 
echo  area  phi  incident-theta  scattered  wave 
echo  area  theta  incident-phi  scattered  wave 
echo  area  theta  incident-theta  scattered  wave 
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TD2 

loss  tangent 

of 

the  insulation 

TD3 

loss  tangent 

of 

ambient  medium 

TH 

theta  angle  for 

far-zone  calculations 

TP 

2n  (6.28318) 

UO 

1.2566B-6 

VG 

antenna  coaplex 

driving  voltages 

VOLT 

list  of  VG's 

X 

X'coordinate 

of 

each  node 

XNP 

list  of  XP's 

XP 

x-coordinate 

for 

near-zone  calculations 

Y 

y-coordinate 

of 

each  node 

YNP 

list  of  YP's 

YP 

y-coordinate 

for 

near-zone  calculations 

Y1 1 

complex  power 

input 

Z 

z-coordinate 

of 

each  node 

ZLD 

complex  load 

at 

a given  node 

ZLLD 

list  of  ZLD's 

ZNP 

list  of  ZP's 

ZP 

z-coordinate 

for 

near-zone  calculations 

ZS 

surface  impedance  of  the  wire 

Z1 1 

antenna  input 

impedance 

\ 
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USER'S  HANUAL 


The  Anteanas-Scatterecs  Analysis  Prograa  (ASAP)  for  thin 
wire  structures  in  a hosogenous  conducting  sediua  perforns  a 
frequency  dosain  analysis  of  antennas  and  scatters.  The 
prograa  is  applicable  in  the  presence  of  either  a perfect  or 
a finite  ground.  This  appendix  vill  describe  and  explain 
the  data  cards  necessary  to  execute  the  coapute  prograa. 
Although  the  prograa  was  written  for  the  IBM  360  coaputer 
systea,  it  can  be  executed  on  another  systea  with  ainor 
aodif ications. 

The  prograa  utilizes  piecewise  sinusoidal  expansion  for 
the  current  distribution  with  Kirchhoff  Current  Law  enforced 
everywhere  on  the  structure.  If  the  structure  contains  end 
points,  the  currents  at  these  points  are  assuaed  to  vanish. 

I.  Prograa  Liaits 

The  thin  wire  assuaptions  are  questionable  and  the 
accuracy  and  convergence  deteriorate  if  the  radius  of  wire 
utilized  for  the  structure  exceeds  0.01  of  a wavelength,  if 
the  longest  sequent  is  greater  than  one-fourth  of  a 
wavelength,  if  the  length  ratio  of  the  longest  and  shortest 
segaents  exceeds  100,  or  if  the  total  wire  length  is  less 
than  30  tiaes  the  wire  diaaeter.  If  a wire  is  bent  sharply 
to  fora  a snail  acute  angle  (less  than  30  degrees)  , the  thin 
wire  Bodel  is  questionable.  It  is  assuaed  that  the  wire 
conductivity  greatly  exceeds  the  conductivity  of  the  aabient 
aediuB.  for  insulated  wires,  the  dielectric  layer  is  assuaed 
to  be  electrically  thin. 

II.  (liniana  Data 

The  ainiaua  data  necessary  to  execute  the  prograa  is: 
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a.  description  of  structure 


b.  radius  of  wire  used  for  the  structure 

The  prograe  will  default  to  the  other  parameters  necessary. 
The  default  parameters  are: 

a.  wire  for  the  structure  is  copper 

b.  frequency  of  operation  is  300  nhz 

c.  homogneous  medium  is  free  space 

A more  detailed  explanation  of  the  defaults  will  be 
discussed  when  the  data  card  for  the  parameter  is  described. 

III.  Outputs 

In  antenna  problems,  the  output  includes  structure 
currents,  impedance(s)  of  f eed  (s) , gain,  polar  radiation 
plots,  and  near  field  calculations.  In  bistatic  scattering 
problems,  the  output  includes  structure  currents,  complex 
elements  of  the  polarization  scattering  matrix,  polar 
reradiation  pattern  plots,  and  echo  areas  produced  by  a 
plane  wave.  For  backscatter ing  problems  the  output  includes 
absorption,  scattering  and  extinction  cross  sections  in 
addition  to  the  outputs  of  bistatic  scattering.  Host  of  the 
outputs  are  suppressed  and  must  be  requested.  Since  the 
program  can  produce  a large  volume  of  output,  care  should  be 
exercised  until  the  user  is  familiar  with  the  outputs. 

IV.  Data  Cards 

The  Analysis  Program  utilizes  free  format  for  the  data 
cards,  that  is,  the  program  utilizes  character  recognition 
to  determine  which  parameters  are  being  read.  Data  placement 
(location)  on  the  input  card  is  not  crltial.  Blank 
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characters,  oa  all  input  cards  but  the  COHHENT  data  card, 
are  ignored  and  naj  be  used  at  the  discretion  of  the  user. 
Since  character  recognition  is  used,  only  the  first  four 
characters  of  the  hey  words  nust  be  present  and  correct. 

The  fornat  for  the  COHHEMT  CABD  utilizes  standard 
POBTBAV  fornat  (i.e.  'C  in  colunn  1 followed  by  at  least 
four  blanks) . The  COHHENT  CiBD  is  the  only  type  of  input 
card  that  position  in  the  data  block  is  critial.  This 
(these)  card(s)  nust  be  placed  at  the  beginning  of  a data 
block.  A data  block  is  a series  of  related  data  cards. 
Several  data  blocks  nay  be  used  to  define  an  analysis 
problen.  This  will  becone  clear  when  the  ternination  cards 
(END,  STOP,  or  CHANGE)  are  discussed.  There  is  no  linit  t(.i 
the  nunber  of  consent  cards  that  nay  be  used.  As  a check  for 
the  user,  all  input  data  cards  will  appear  on  the  output  as 
they  appear  in  the  input  deck. 

The  fornat  of  other  data  can  be  of  one  of  two  forns: 

a.  type  of  card  (option  1/option  2/ ) 

b.  paraneter  (value)  . 

The  type  of  fornat  to  use  will  be  apparent  as  the  individual 
data  cards  are  discussed. 

The  nunerical  values  for  the  paraneters  nay  be  stated  in 
any  one  of  the  fol lowing  forns.  The  progran  will  translate 
the  nunber  to  the  proper  forn  for  the  specified  paraneter, 
either  fixed  or  floating  point.  All  of  the  following 
exanples  have  the  sane  value. 

0.0001  or  .0001  or  100.0  or  1000  or  .1H  or  0.1H  or  .0000001K 

-6  -3  3 

0 » 10  H » 10  K » 10 
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1*  JIEl  This  card  is  used  to  define  the  paraaeters 
associated  with  the  wire  utilized  by  the  thin  wire 
structure.  Two  options  are  available  and  are  defined  as: 

RADIOS=value  of  the  radius  of  the  wire  in  aeters 

COMDUCTIVITY=value  in  aeganhos  per  aeters  . 

The  wire  data  card  aust  appear  in  the  first  data  block  to 
define  wire  radius.  The  default  value  of  the  conductivity 
is  50  aegaahos/aeter  (copper) . 


IRE(  RADI0S=.001/  CONDU=28.5) 


msULATION  This  card  is  utilized  to  define  the  paraaeters 
associated  with  the  insulation  of  the  wire  used  for  the 
structure  to  be  analyzed.  If  this  card  is  onitted,  the 
program  assumes  that  the  structure  is  uninsulated.  Four 
options  are  availaole  and  are  defined  as: 

RADI(JS=value  of  outer  radius  in  meters 

CONDOCTI VITY^value  in  aicromhos  per  meter 

DI ELECTRIt:=value  of  relative  dielectric  constant 

LOSS  TANGENT=value  . 

The  conductivity  and  either  the  relative  dielectric  constant 
or  the  loss  tangent  (but  not  all  three)  options  may  be 
stated. 

( RADIDS=.015/  COND=7./DIEL=5) 
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3.  EXTERIOR  MEDIUM  This  card  is  utilized  to  describe  the 
homogeneous  medium  surrounding  the  structure.  If  the  medium 
is  free  space,  this  card  may  be  omitted.  Three  options  are 
available  and  are  defined  as: 

DI ELECTRIC=value  of  relative  dielectric  constant 

CONDUCTIVITY=value  in  micromhos  per  meter 

LOSS  TANGENl=value  . 

As  with  INSULATION  card  state  either  conductivity  or  loss 
tangent. 


XTE (LOSS=.  45) 


4.  DfiSfRiPTION  This  card  is  utilized  to  describe  the  shape 
of  the  wire  structure  to  the  program.  The  user  must  divide 
the  wire  structure  into  segments  of  the  appropriate  length 
and  number  each  node  starting  at  one.  A node  is  a point 
where  a segment  begins  or  ends.  A maximium  of  four  segments 
can  meet  at  any  given  node.  An  isolated  wire  must  contain 
at  least  two  segments  and  three  nodes.  Thus  the  DESCRIPTION 
CARD  must  show  at  least  3 consecutive  nodes  for  all  portions 
of  the  wire  structure.  The  structure  is  described  by 
stating  the  node  numbers  that  each  segment  connects.  The 
description  of  a sguare  loop  might  appear  as; 

DESCRIPTION ( 1-2/2-3/3-4/4-1)  . 

The  description  of  a dipole  and  reflector  might  appear  as; 

DESCRIPTION (1-2/2-3/3-4/4-5/6-7/7-8/8-9/9-10)  . 

If  the  description  will  not  fit  on  one  data  card  continue  on 
the  next  card  as  if  the  previous  card  were  longer.  The 
dipole  example  might  appear  as: 

DESCRIPTION ( 1-2/2-3/3-4/4-5/ 

6-7/7-8/8-9/9-10)  . 

Note  that  the  last  character  on  the  card  to  be  continued  is 
a slant  (/) . As  many  cards  as  necessary  may  be  used.  The 
maximum  numuer  of  nodes  permitted  is  fifty.  If  ground  plane 
is  present,  the  maximum  number  is  twenty-five.  If  a ground 
plane  is  present  and  the  structure  touches  the  ground  plane, 
the  lowest  node  numbers  MOST  be  used  for  the  touching  nodes. 
That  is,  if  the  structure  touches  the  ground  plane  at  two 
points,  node  numbers  1 and  2 MOST  be  assigned  to  these 
nodes. 


/ DESCR (1-2/2-3/3-4/4-1) 

/ 
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5.  GEOHETRY  This  card  is  used  to  state  the  physical  location 
in  rectangular  coordinates  of  each  node  of  the  DESCRIPTION 
CARO  . The  rectangular  grid  is  in  units  of  neters.  If  node  1 
is  located  at  x1,y1,z1  and  node  2 at  x2,y2,z2  and  node  3 at 
x3, y3 ,z3 ,etc. , the  GEOHETRY  CARD  night  appear  as: 

GE0METRY(x1,y1,z1/x2,  y2,z2/x3,  y 3 ,z3/ ) 

As  with  the  DESCRIPTION  CARD,  continuation  cards  are 
permitted. 


1/-. 1,0.1/-. 1,0-. 1/. 1,0,-. 1) 
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FEED  For  antenna  analysis  the  feed  point (s)  and 
voltage(s)  nust  be  stated.  In  the  foreaention  dipole  and 
reflector  example  if  the  feeds  were  at  node  2 with  a voltage 
source  of  .5  at  an  angle  of  -90  degrees  and  at  node  4 with  a 
voltage  source  of  .5  at  an  angle  of  ■••90  degrees  the  FEED 
CARD  night  appear  as: 

FEED (2,. 5, -9 0/4,. 5, *90) 

The  order  of  the  information  for  each  voltage  source  is  node 
nuaber,  magnitude,  and  phase  angle.  This  order  is  repeated 
until  all  sources  are  stated.  If  the  source  infornation  will 
not  fit  on  one  card,  use  another  card  similiar  to  the 
initial  one;  that  is,  repeat  the  word  "FEED”.  If  only  one 
voltage  source  is  applied  to  the  structure,  only  the  node 
number  must  be  stated.  In  the  dipole  example,  if  the  drive 
is  at  node  3,  the  FEED  CARD  might  appear  as: 

FE£D(3) 

A default  source  of  one  volt  at  zero  degree  phase  is 
assumed.  Voltage  sources  should  only  be  stated  for  nodes 
with  only  two  segments. 


FEED (2,. 5, -9 0/4,. 5, ♦90) 
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7.  LOAD  This  card  is  used  to  describe  the  loads  to  be  placed 
at  various  locations  on  the  structure.  The  format  for  this 
card  is  siniliar  to  that  of  the  FEED  CARD,  that  is,  the 
word  "LOAD"  is  used  in  the  place  of  "FEED".  The  order  of 
the  information  on  the  card  is  the  same.  Since  this  card  is 
frequency  dependent,  it  must  be  changed  if  the  frequency  of 
operation  is  changed.  No  default  parameters  are  available. 
The  structure  is  assumed  unloaded  unless  this  card  is  used. 
Once  the  structure  is  loaded,  it  will  remain  loaded  for  the 
remainder  of  the  data  blocic  series.  To  unload  the  structure 
the  following  card  may  be  used: 

L0AD(-1) 


LOAD (1,1 20, -45/3, 120, ♦45) 


B.  ODTPOT  This  card  is  used  to  request  output  data.  Host  of 
the  output  is  in  tabular  fora.  Hore  than  one  OOTPUT  CARD 
is  peraitted  per  data  blocic,  but  not  for  the  saae  type  of 
output.  If  only  the  antenna  input  iapedance,  antenna 
efficiency,  or  tiae-average  power  input  is  of  interest,  no 
OOTPtJT  CARD  is  necessary.  These  paraaeters  are  autoaat ically 
printed  if  a FRED  CARD  or  GENERATOR  CARD  is  utilized.  One 
or  Bore  of  the  following  options  aay  be  used  to  request  the 
various  outputs  available. 

PAR  PIELD=phi  initial,  phi  final,  theta  initial,  theta 

final 

This  option  gives  the  coaponents  of  the  electric  field 
intensity  in  the  far  field  as  phi  and  theta  varies  between 
liaits  specified  in  one  degree  divisions. 

BACKSCATERlNG=phi  initial,  phi  final,  theta  initial, 
theta  final 

This  option  gives  the  absorption,  scattering,  and  extinction 
cross  sections,  and  the  coaplex  eleaents  of  the  polarization 
scattering  aatrix  for  an  incident  plane  wave  illuainating 
the  structure  froa  the  spherical  direction  of  phi,  theta  as 
both  vary  between  liaits  specified  in  one  degree  divisions. 

BISTATIC=phi  inital,  phi  final,  theta  initial,  theta 

final 

This  option  gives  echo  area  and  the  complex  eleaents  of  the 
polarization  scattering  aatrix  for  an  incident  plane  wave 
illuainating  the  structure  froa  the  spherical  direction  phi, 
theta  final  of  the  backseat tering  output  option,  reradiated 
in  the  phi,  theta  direction  as  both  vary  between  liaits 
specified  in  one  degree  divisions.  A bistatic  output  request 
aust  be  accoapanied  with  a backscattering  request  in  the 
saae  data  block. 

STEP'walue  in  degrees 
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This  option  will  cause  any  of  the  above  output  options  to  be 
stepped  at  a different  interval  size.  That  Is,  if  one  of  the 
above  options  is  to  be  stepped  at  ten  degrees  intervals,  use 
this  option.  This  option  overrides  the  one  degree  stepping. 

CURRENT 

This  option  gives  the  currents  on  the  structure  which  are 
produced  by  the  feed/generator  voltages  and/or  the  incident 
plane  wave  of  the  backscattering  request. 

NEXR=x1,yl,z1 

or 

NEAR=(x1,y1,z1/x2,y2,z2/x3,y3,z3/etc ) 

This  option  gives  the  value  of  electric  field  cosponents  in 
the  near  field  for  the  antenna  at  the  point  or  points 
specified. 


'output  (FARF=4  5,50,25,50) 


£ii22  This  card  will  produce  aornalized  polar  plots  in  the 
specifed  plane  for  the  stated  option.  The  plane  is 

specified  by  stating  either  ''PHI= ” or  •*THETA= . The 

PLOT  CABO  overrides  the  limits  of  the  ODTPUT  CARD  for  the 
same  option.  If  only  a normalized  pattern  is  of  interest, 
only  a PLOT  CARD  is  necessary.  If  a table  of  values  and  a 
normalized  pattern  is  desired,  both  a PLOT  CARD  and  OOTPUT 
CABO  must  be  used.  Only  one  PLOT  CARD  is  permitted  per 
data  block.  The  following  pattern  plots  are  available: 

PAR  PIELD/plane 

This  option  will  plot  the  far  field  intensity  for  each 
component  of  the  electric  field. 

BACKSCATTERlNG/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix. 

BISTATIC/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix  produced 
by  the  incident  plane  wave  stated  by  final  limits  of  the 
backscattering  option  of  the  output  request. 
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10.  GBOQMD  This  card  is  used  to  describe  the  ground 
paraseters  if  a ground  plane  is  present.  If  no  ground  plane 
is  present,  the  structure  is  assuned  to  be  in  free  space  or 
the  hoBogeneous  aediui  of  the  EXTERIOR  HEOIOH  data  card. 
Seven  options  are  available  and  are  defined  as: 

PERFECT 

This  option  will  analyze  the  structure  over  a perfect  ground 
plane. 

GOOD 

This  option  will  analyze  the  structure  over  a good  ground 
plane  where  the  conductivity  of  the  ground  is  .02  ahos/aeter 
and  the  relative  dielectric  constant  is  30. 

POOR 

This  option  will  analyze  the  structure  over  a poor  ground 
plane  where  the  conductivity  of  the  ground  is  .001 
ahos/aeter  and  the  relative  dielectric  constant  is  U. 

SEA 

This  option  will  analyze  the  structure  over  salt  water  where 
the  conductivty  of  the  water  is  4.  ahos/aeter  and  the 
relative  dielectric  constant  is  80. 

HEIGHT^value  in  meters 

This  option  will  analyze  the  structure  with  origin  of  the 
GEOMETRY  Card  this  height  above  the  ground  plane.  The 
lowest  point  of  the  structure  aust  not  lie  below  the  ground 
plane.  It  aay  lie  on  the  ground  plane. 

C0I1DUCTIVITT=  value  in  ahos/aeter 

This  option  is  used  to  state  the  value  of  conductivity  of 
the  ground  plane  if  the  default  values  aentioned  above  are 
not  utilized. 
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DIELBCTBIC>  value 


This  option  is  used  to  state  the  relative  dielectric 
constant  of  the  ground  plane  if  the  default  values  nentioned 
above  are  not  utilized. 


'G  aOOHD(HEIG=1  0/CO  1ID=.  002/DI  EL  = 10) 


ll*  IlXfilil  fog  CALCaiATION  This  card  is  used  to  state  the 
nuebec  of  intervals  to  be  used  for  calculating  the  eleaents 
of  the  iapedance  eatrix  with  Siapson ' s-rule  integration.  A 
large  value  for  the  nuaber  iaproves  the  accuracy  at  the 
expense  of  greater  execution  tiae.  For  aost  probleas  a 
suitable  coabination  of  speed  and  accuracy  is  obtained  with 
a value  of  four,  the  default  value.  If  the  rigiorous 
closed-fora  iapedance  expressions  in  teras  of  the 
exponential  integrals  is  desired,  set  this  value  to  zero. 


IHTEH»AL=value 


12.  GENEBAlOfi  This  card  is  siailiar  to  the  FEED  CARD  in  use, 
except  that  the  segeent  nuabecs  are  stated  instead  of  the 
node  nuabers.  This  is  useful  if  three  or  four  segaents  neet 
at  a node.  The  positive  terainai  of  the  generator  is 
connected  to  the  specified  segaent  such  that  current  is 
forced  in  the  the  positive  direction.  The  positive  direction 
of  current  flow  is  froa  the  first  stated  node  nuaber  of  that 
segaent  toward  the  second  stated  as  ordered  on  the 
OESCBIPTIOM  CARD. 


i 
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13.  iHPgPAHCE  This  card  is  siaiiiar  to  the  LOAD  CAfiD  in  use, 
except  that  the  segaent  nuabers  are  stated  instead  of  the 
node  nuabers.  As  with  the  GENEPATOB  CABD,  this  is  used  if 
three  or  four  segaents  are  connected  to  a node.  The 
impedance  will  be  connected  to  the  positive  terainal  of  the 
specified  segaent.  The  foraat  of  this  card  is  the  same  as 
the  LOAD  CABO. 


/ IHPE(1,120,-45/3, 120,*45) 


14.  FREQUENCY  This  card  is  used  tO  specify  the  frequency  in 
aegahertz  if  it  is  to  be  other  than  the  default  value  of  300 
HHz. 


/ FREQ(12.5) 


118 


I 


15,  CHAMGg  This  card  at  the  end  of  the  data  block  signals 
the  progran  that  the  follotring  data  cards  are  changes  to  tne 
previously  read  data,  for  the  next  run.  If  a "CHANGE  CARD" 
is  used,  the  outputs  aust  be  requested  again  in  the  next 
data  block.  The  "CHANGE  CARD"  cannot  be  used  to  change 
"DESCRIPTION  CARD"  or  "GEOHETRY  CARD"  data  vhen  operating 
with  a "GROUND  CARD".  Use  an  "END  CARD"  to  make  changes  when 
a "GROUND  CARD"  is  used. 


16.  EJfD  This  card  signals  the  prograa  that  this  is  the  end 
of  a data  block  series  and  to  reinitialize  data  for  the  next 
problew.  An  "END  CARD"  cannot  be  used  with  a "CHANGE  CARD". 


17.  STOP  This  card  signals  the  prograa  that  all  of  the  data 
cards  have  been  read  and  to  terminate  itself  when  execution 
is  completed.  This  card  must  be  used  as  the  last  card  in 
place  of  the  "END  CARD"  of  the  last  data  block  series.  A 
"STOP  CARD"  cannot  be  used  with  an  " END  CARD"  in  the  same 
data  block. 
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C AN  EXAMPLE  PROBLEM 

C 

C V ANTENNA 

C 

WI  RE  (RADIUS’^!  M) 

GE0M(0,-. 18. 1 8/ Ot 09, f. 09/0, 0,0/0, 0.09,. 09/0,. 18,. I 8) 
OESC(  1-2/2-V3-4/4-5) 

FEED( 3) 

0UTPUT(FARF  = <»5,50,65,80/STEP  = 5) 

CHANGE 

0UTPUT(BIST  = 45,A5,45,<»5/BACK  = 0,0, 10,12) 

OUTPUT(CURRENT) 

CHANGE 

C 

C CHANGE  STRUCTURE  SHAPE  TO  OIPQLE 

C 

GE0M(0,-. 25 ,0/0,-. 125,0/0,0,0/0, .125, 0/0, .25,0) 
PLOT(FARF/PHI=90) 

GROUND! HEIGHT *.2 5/GOOO ) 

STOP 


THE  ABOVE  DATA  DECK  WILL  PRODUCE  THE  OUTPUT  ON  THE 
FOLLOWING  PAGES. 
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FAR  FIELDS  FOR  PHI  VARYING  FRQH  AS.O  TO  SO.OANO  THETA  VARYING  FROM  6S 
IN  SItPS  OF  S.O  DEGREES. 
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CALCULAf IONS 


I 


UKO'^eOO 

UJ#<N»00 

uJ9‘iA9*0 

ujr*j#r-o 

«/>•••• 

i/I  • • • t 

i/1  • • • « 

i/I  • • • • 

tfisoDGon 

« o-«^o 

< iQIDIDO 

4099*0 

Till 

TIM  1 1/1 

XIII 

TiM^irv 

a 

1 

a. 

a—  1 1 

O 0-^0 

o o«#-« 

O 0—0 

o o#— 

UJUiOOO 

UJUJOOO 

UJUJOOO 

UJUJOOO 

»*«iO 

•MO  1 

•MO 

KlO  1 

o 

M ^lUUJUJ 

M— lUlllllU 

o 

MOlilUJUi 

•— mojinuj 

• 

9*000 

• 

^♦-aiOdo 

ovcorsiO 

o 

<»«<MOr<4 

<^0*000 

«m/\J04M 

4•«C^9*0 

X Zi^om 

XZ9*#0 

M 

ZZOOO 

XZO'fMO 

QC0)»*'O<^ 

flCOCM^O 

QCO'^'OO 

«K09*i/l0 

H 

0<l^— ►-O 

0<9*»— O 

M 

o<^-— ^o 

5^9i^-BO 

< 

zx  • • • • 

zz  • • • • 

< 

zz  • • • • 

zx  • • • • 

oooo 

oooo 

V 

oooo 

OOOO 

UJ 

UJ 

I 

^fsjru 

#(^■# 

I 

rgfsjrg 

#9*# 

ujOOO 

UJOOO 

V 

UiOOO 

UJOOO 

O 1 f 1 

O I 1 1 

0 1*1 

O 1 1 1 

a 

DUiUJuJ 

OUJUiUi 

O 

OUJUJUJ 

OUlUiUJ 

z 

♦-O'-#®* 

^.iNOOrg 

Z 

♦-  Om-O 

w— 9rNi 

< 

< 

Mi/M/\in 

Z-^a)«u 

XmC'm 

Z-«co— 

2909 

o 

O 

0<MO*M 

0999 

4 

<r<|#mo 

• 

4miAmo 

o 

X • • * * 

X • • • ■ 

o 

z • * • • 

X • • • • 

oooo 

OOOO 

OOOO 

oooo 

H 

#— # 

N 

^9“# 

M 

>-ooo 

VO— o 

M 

vooo 

vooo 

T 

a 1 t 1 

oc  1 1 1 

X 

ec  1 1 1 

ec  1 1 1 

A. 

<UJUJUJ 

<UJUJUJ 

a. 

<UJUiUJ 

4UJUJUJ 

Z#fiO»»> 

ZJMOrg 

Z9‘if10* 

Zm— #■ 

• 

Mh.N.h» 

UJ 

— lO— >0 

» 

•— CT'Clj* 

UJ 

— #9# 

(/> 

oo-«o 

> 

Oh'i/'f- 

i/I 

> 

0 909 

UJ 

1 UJ 

<0'>#0' 

< 

<QOirtoo 

UJ 

UJ 

<9#9‘ 

4 

4—9  — 

1 > 

X'#^#0 

X 

XlMO'IMO 

> 

Z#'9#0 

X 

X m— mo 

C3 

• < 

•-••••• 

O 

< 

*•«••• 

Z 

X 

1 1 1 o 

o 

Olio 

z 

X 

1 1 I O 

O 

Olio 

UJ 

< 

Ul 

o 

•NfMfM 

rsj 

#9*# 

o 

•MfSl/M 

•M 

#9*9 

UJ 

ooo 

OOO 

u- 

OOO 

M 

OOO 

z 

NJ 

1 1 1 

oc 

1 1 1 

z 

•M 

1 1 1 

OC 

1 1 1 

UJ 

1 M 

UJUiUJ 

< 

UJUIUJ 

UJ 

UJUJUJ 

4 

UJUJUJ 

a 

) ec 

o 

oc 

_jf<irgit| 

•1 

^##9 

M 

< 

o 

<###• 

M 

< 

o 

49(M9 

u 

) mJ 

a. 

010090 

u 

Ulh-f^l*^ 

& 

UJ#^*’# 

z 

o 

OC^rrt^ 

er^#^ 

z 

o 

oc  (Ximoo 

a. 

< 

— #— o 

a. 

4 

—9—0 

• • • • 

K 

• t • « 

• • B • 

H- 

• • • • 

UJ 

1 M 

OOOO 

UJ 

1 OOO 

Ui 

OOOO 

UJ 

1 OOO 

X 

r 

I 

T 

X 

I 

K 

a. 

UJ 

UJ 

V 

a 

UJ 

V 

UJ 

Orum-^IA 

O'Mm-#!/! 

OlMfO#!/! 

0<Mm#9 

Of 

UJ 

o 

UJ 

O 

oc 

UJ 

o 

Ui 

o 

o 

X 

z 

X 

z 

O 

r 

z 

X 

z 

u. 

u. 

V 

»- 

UJO<0#<D 

ujO#r>JO 

UJOO*U^(^ 

ujO^#i^ 

o 

> 

i/1  • • t • 

>• 

i/l  • • • • 

O 

V 

i/I  • • • • 

V 

i/l  • • • • 

z 

<D 

4 oao<D  <x> 

o 

<00—9* 

z 

(9 

409000 

o 

40099* 

X III 

I ^ 

X lit 

I <M— 9 

K 

O 

a 

o 

a.  — 1 

CC 

Q 

o. 

o 

a.  — 1 1 

UJ  UJ 

UJ 

UJ 

UJ 

UJ 

H* 

O 

o 

►— 

U 

u 

D 

0^0 

o 

O#-^ 

N- 

3 

o-«o 

3 

09  — 

< 

O 

OUJ  ooo 

o 

OUJ  OOO 

4 

O 

Ouj  ooo 

Q 

OUi  OOO 

o 

1 z 

UjO 

z 

UJO  I 

o 

Z 

UiO 

Z 

ujO  I 

4/J 

•^O  UJUJUJ 

•M 

KJO  UJUJUJ 

i9 

♦*J3  UJUJUJ 

•M3  UJUJUJ 

<00'0 

9*000 

M—  COOCD 

mV  COJMO 

UJ  iA 

j—  rsjO'^ 

i/I 

9*(DO 

UJ 

t9 

OJO/M 

i9 

<^0*O 

> 

<Z  OO'^I 

K 

«Z  9*<#0 

> 

V 

4Z  Kiom 

V 

4Z  9*rMO 

< 

z 

TO  mo'^ 

z 

xo  9*mo 

< 

z 

XO  mom 

Z 

Zo  9*90 

X 

UJ 

oc 

UJ 

«<09*— •« 

X 

UJ 

gC40l^— 

UJ 

a40<^'^— 

1 

ox  • • • • 

« 

OZ  f • * • 

1 

OC 

OX  • • • • 

OC 

OZ  • • • f 

UJ 

ct 

z oooo 

oc 

z oooo 

Ul 

oc 

r oooo 

eC 

Z OOOO 

z 

D 

3 

z 

3 

3 

u 

4>jr*jAJ 

U 

< 

o 

iMrg/M 

O 

9«^9 

UJ  ooo 

UJ  ooo 

UJ  OOO 

UJ  ooo 

a. 

O III 

O III 

a. 

O III 

O III 

O ujujUJ 

3 UJUJUJ 

m UJUIUJ 

3 UJUJUJ 

T 

► ®‘#(^ 

^ — OOru 

z 

— 91*19 

u-  — 9fM 

— ### 

Oi/10 

V 

lAUM/l 

M m— m 

M 

Z o^dO-^ 

z 09*0 

M 

Z — <n— 

Z 909 

X 

O 

0#<<l 

X 

lO  fNjmrg 

O 999 

<Om-#«*> 

<0O#ff| 

4 0m#m 

40m9m 

o 

X • * • • 

Z • t • • 

o 

X • * • • 

z « • • * 

UJ 

OOOO 

oooo 

UJ 

oooo 

OOOO 

»• 

V 

'T— # 

< 

9^9 

M 

► ooo 

V 0—0 

•m 

V ooo 

> OOO 

o 

at  III 

« III 

V 

OC  III 

oc  III 

a 

4 UiujUJ 

< UJUJUJ 

o 

4 UJUJUJ 

4 UJUJUJ 

*/> 

Z #91^ 

Z IMOJM 

s/t 

Z (^t/>9* 

Z m— 9 

— <9— « 

*A 

— 9*90* 

— 999 

o o««o 

O 

< 

O 9^9 

< ‘>#0* 

4 Oino 

4 9*#9* 

< m9  — 

ico#-o^ 

ZOIMCMM 

i/l 

ZO#9«# 

zo(i>— m 

• 

— • • • • 

V 

— • • • • 

z 

O 1 1 1 

OOl  1 

Z 

O 1 M 

OOl  1 

UJ 

UJ 

« 

• 

<yrsi«M 

oc 

furuAJ 

99*9 

? 

000 

1 1 1 

??? 

5 

OOO 

1 1 1 

OOO 

u 

UJUJUJ 

UJUJUJ 

o 

UJUJUJ 

UJUJUJ 

^ rui««ril 

•J  Ooi# 

^ mium 

^ 999 

z 

« 

Z 

4 9fM9 

o 

144 

UJ  09*0 

u 

UJ  Kmr» 

Ui  9r»9 

z 

« <0^9 

Z 

g£  — m— 

QC  090 

O^#  — 

4 

om^m 

0«U9^ 

• • • • 

• • • • 

« 

t t • • 

• • • • 

m 

OOOO 

o 1 oe 

O 

oooo 

O 1 oo 

UJ 

O— 

Ul 

o— fU99 

9 

o 

o 

o 

z 

z 

z 

z 

^•4um9 

UJ 

Ul 

UJ 

UJ 

4.A 

*A 

i/I 

i/I 

126 


! 


UJ(^'00'0 

ui*^-^9‘0 

• • t • 

!/>•••• 

4D4D<DO 

Till 

TINJOJ® 

O 

1 1 

o o«^o 

o oo— 

UiUjOOO 

UJ®000 

"UO 

►UO  1 

o 

»«OUJUJUJ 

— ^UittiUi 

• 

•J^ooo 

o 

o 

zzooo 

ZZ0'(^O 

OCO'^'OO 

flCOOmo 

H 

5<»^— 

o<a— O 

zz  • • • • 

zz  • • • • 

oooo 

oooo 

UJ 

Z 

rvj<MiN 

ooo 

ujOOO 

®ooo 

O 1 1 1 

O 1 1 1 

o 

OUiUJUi 

o®  UJ® 

z 

Kirkoi/\ 

►-fMiMO 

< 

mm-CA'O 

— iNjrgr\j 

Z*4C0>^ 

Z(^■^^C^ 

O 

OO'OO' 

• 

<o,tTno 

<0000 

o 

X • • • • 

X * * * * 

oooo 

oooo 

N 

moo 

oao 

M 

>ooo 

VOQO 

z 

oc  1 1 1 

« 1 I 1 

o 

<UiUJUJ 

<®®® 

ZlAf*.^^ 

ZOOO 

• 

mm^O^ 

UI 

— rsiO<*u 

«/) 

> 

oo\^y> 

UJ 

UJ 

<C>(AO» 

< 

<ooo 

o 

> 

3 

XOIMOO 

o 

*«•••• 

— • • • • 

z 

X 

1 1 1 o 

o 

Olio 

UJ 

o 

IMIMCM 

*<u 

OO'O 

UJ 

ooo 

mm 

OOO 

z 

Kl 

1 1 1 

oc 

1 1 1 

UJ 

mm 

UJUJUJ 

< 

®®® 

o 

oc 

oiJ>otrk 

-lOOiA 

M 

< 

o 

<0—0 

o 

o 

uuf^o^* 

a. 

®Ou^O 

z 

o 

ao^o 

o 

0^00 

< 

fMoruO 

• • • • 

• • • • 

Ui 

M 

oooo 

UJ 

1 ooo 

X 

r 

Z 

►- 

o 

UJ 

w 

® 

OlMO<#U> 

QlMOOrf> 

ec 

UJ 

o 

UJ 

O 

o 

z 

z 

Z 

z 

u. 

UJOO>0(^ 

®0  — O' 

O 

>- 

tA  • • • • 

> 

m • • • • 

z 

o 

<00000) 

(fi 

<00i/>0' 

•• 

Z ill 

Z 

ac 

o 

o 

a — 1 I 

UJ 

UJ 

UJ 

►- 

o 

u 

o 

o««o 

3 

oo  — 

< 

a 

OUJ  ooo 

O 

OUJ  ooo 

o 

z 

UJO 

Z 

®o  1 

tA 

M 

*uo  Uiujuj 

mm 

®UJ® 

mm^  OOO 

— ^ 0—0 

UJ 

lA 

OOO 

tA 

O'— o 

> 

<Z  ooo 

<Z  ^O'O 

< 

z 

ZO  ooo 

z 

ZO  O'U^O 

'S 

UJ 

UJ 

ac<o^— 

1 

oc 

OZ  • • • • 

OC 

OZ  • • • • 

Ui 

oc 

z oooo 

€L 

z oooo 

z 

o 

D 

4 

u 

iMrsj'M 

U 

ooo 

o 

UJ  ooo 

® ooo 

o 

O III 

O III 

UJUJUJ 

O ®uj® 

z 

^ tnmtn 

^ ojruo 

'O^J^ 

mm 

mm 

z -i®— 

Z O'lr^ 

3 

o <MOfU 

O O'00‘ 

<oo<o 

<Oo®o 

O 

z • • • * 

z • • • • 

Ui 

oooo 

oooo 

o^**< 

0^0 

M 

> ooo 

> ooo 

V 

oc  III 

Ot  III 

o 

< UJUJUJ 

< ®®® 

•A 

z 

z ®trt® 

tA 

M .00® 

— ru®<M 

< 

O 

o mO'iA 

< 9«00» 

< ooo 

tA 

zo®o® 

zooMm 

»••••• 

— • • • • 

z 

O 1 1 1 

oo  1 1 

UI 

« 

J^JSOU 

OO'O 

% 

000 

1 1 1 

ooo 
* T 1 

000 

1 < I 
UJtUUJ 

OOD’d'O 

f ^•«fMiM 
> t • • 

LJ  ooo 


I 

< OoO 
»-  III 
UJ  UiiUU) 
X JOO»^ 

ujO'#  O 

CL 

• • • 

or>o 


M I flOoOQO 

oc  < ooo 

^ ^ III 

< Ui  U^UiUJ 

Z Z-IO'^^O 

0 uia>  09^ 

2 ac<7‘l/^f«^ 

•-  orsi  — 

ac  • • « 

01  oo  I 


uJZ^ 

^o^- 

►-— < I 

<»-o 

o<</> 

vA^s*  I 

ooc?^ 

<^UilAi 

O^OZ 
O— ►- 
0.0  I 

r« 

Q— Z 

-J-'O. 


000 

1 I I 
liiuju; 

OOOIM 

Z4’4'0 


^ iAm«A 
lu 


3 


UMAIUi 

«J 

ttl  0tf>0 

« o^>o 

• • • • 

O I oo 

<11 


iUOlUl 

f^ooO 

OOfVii^ 

Z«o»« 

I <n«i 

• • • 

**  III 
Z 

0 

1 

•-  ooo 

Z III 
O 

<tf>00 

Ui<Om 

aeN>^>« 


4000 
• • • 
mO*^«M 


oz>  ooo 
u33xooo 


m 

o 


• 4 

1 • « « 

► • 

« ^ O 
Ol  I 

• 

< Z.^UJ 

• 

Z 

« 

zo 

UiU> 

# 

at  • 

O ocrsj 

« 

Orsi 

Z m 

• 

M <M 

■■ 

» 

U 

CL  • 

ac*/* 

« 

Z N 

Ui  o 

UJZ 

• 

— < 

h>0 

« 

h- 

• 

OtUJ 

<UJ 

• 

TX 

UflC 

CJ« 

« 

i/)UJ  O 

# 

>-  1 

« 

oo 

Z^  Ui 

uu 

# 

z • 

o<  or* 

M_i 

• 

— o 

U*-U  <rM 

# 

oc 

<o 

« 

Ui 

1 mcD 

• 

^ » 

l/t 

• 

^^z<  • 

# 

<X 

^oCui^  1 

A 

• 

UA. 

— <OUJ 

• 

0-i**Z 

• 

l/> 

Oo^  ^ 

• 

az  1 o 

• 

MM  1 

• i 

> • • 1 

» • 

»-Ui 

OMX^Ui 

«> 

-i  a<co 

lAi  UJC^ 

M 

IL  •# 

tta> 

M 

Z 

u • 

flC^ 

••  1 

ac 

a.  CL 

2 

O • 

M mic 

>X 

«•«  o 

MOm  • 


COMTINUC  EXECUTION  MITH  THE  FOCLOMINC  ADDITIONS  AND/OR  CHANCES 


oooo 

• • • • 

oooo 


o oo 
o oo 


2 3 8| 

lAo 

Mih 

o • • • • 
W I ooo 


VIOOO 

UiOpO 

>000 

oo«> 


oooo 

• • • • 

oooo 


8* 


oooo 

• • • • 

oooo 


# ^ 

« « 

UI 

* a 

* ^ 

: 5, 

s 

* X 

00  0 

• UlUi 

00  0 

# 

• oe^ 

5 

0 

• 

• «»Ui 

X 

X 

UlUI  UI 

» 

« 

• X 

0 

00  0 

• 

* ^s* 

lA 

M 

00  0 

fU 

• 

* X«A 

K>OiA  tf> 

• 

• 

• <0 

UI 

< 

IA<M  fU 

• 

# 

« 

u 

»u^O-< 

0 

• 

M 

5 

• • • • 

«••• 

#^X^4C 

s 

1 1 00 

0 

UI 

* 

UiU>^ 

> ^Ul 

fU 

NnOmS 

• 

X«M 

• 

O 

00 

liJ>  Ui 
Ui  ^O 
Ik  <0 
ujO 
« Oto 

1 ■: 

Ui  o 


• # 

• • 

« • 

* • 

• o • 

SUI  « 

^ * 

t n s 

• o # 

• UI  • 

• « « 

9 S I 

* *.  i 

• a • 

• 0*0 
• • • 

» •p 

• *0 

• * 

• # ■ 

• 

•••••X 


£ou 

888^ 

ssP 


oooo 

• • • • 

oooo 


UIX 


li 


LIST  OF  REFERENCES 


1,  Richmond,  J.H.,  "Radiation  and  Scattering  by 

Thin-Hire  Structures  in  the  Complex  Freguency 
Domain,"  Report  290?-10,  July,  1973,  The  Ohio  State 
Unversity  Electroscience  Labortory,  Deoartment  of 
Electrical  Engineering:  prepared  under' Grant  NGL 
36-008-138  for  National  Aeronautics  and  Space 
Administration,  Langley  Research  Center,  Hampton, 
Virginia  23365. 


2. (a)  Richmond,  J.H.,  "Computer  Program  for  Thin-Wire 

Structures  in  a Homogeneous  Conducting  Medium,"  NASA 
Contractor  Report  CR-2399,  June  1974,  for  sale  by 
the  National  Technical  Information  Service, 
Springfield,  Virginia,  22151,  Price  13.75. 


(b)  Richmond,  J.H.,  "Computer  Program  for  Thin-Wire 
Structures  in  a Homogeneous  Conducting  Medium," 
Report  2902-12,  August  1973,  The  Ohio  State 
University  ElectroScience  Laboratory,  Department  of 
Electrical  Engineering  prepared  under  Grant  NGL 
36-008-138  for  National  Aeronautics  and  Space 
Administration,  Langley  Research  Center,  Hampton, 
Virginia  23665. 


3.  Richmond,  J.H.  and  Geary.  N.H.,  "Mutal  Impedance  of 

Nonplanar-Skew  Sinusoidal  Dipoles,"  Report  2902-18, 
August  1974,  The  Ohio  State  university 
ElectroScience  Laboratory,  Department  of  Electrical 
Engineering. 


4.  Miller,  E. K. , Poggio.A.J..  Burle,  G.J.,  and  Selden, 

E.S.,  "Analysis  of  Wire  Antennas  in  the  Presence  of 
a Conducting  Half  Space:  Part  I.  The  Vertical 
Antenna  in  free  Space,"  Canadian  Journal  of  Physics, 
50,  pp  879-888. 


5.  Hiller,  E. K. , Poggio.A.J.,  Burle,  G.J.,  and  Selden, 

E.S.,  "Analysis  of  Wire  Antennas  in  the  Presence  of 
a Conducting  Half  Space:  Part  II,  The  Horizontal 
Antenna  in  free  Space,  Canadian  Journal  of  Physics, 
50  pp  2614-2627. 


w 


